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ABSTRACT 


Several  issues  assoeiated  with  grid  resolution  and  the  sub-grid  seale  (SGS) 
parameterizations  in  a  mesoscale  model  are  addressed  in  this  study.  Of  partieular 
eoneern  is  related  to  inereasingly  high-resolution  mesoseale  atmospheric  numerical 
models,  in  that  sub-grid  scale  parameterization  of  atmospheric  processes  becomes  unclear 
when  the  grid  resolution  becomes  comparable  to  the  length  scale  of  the  phenomenon. 
These  issues  are  studied  through  the  careful  analysis  of  in  situ  observations  and 
mesoscale  model  testing.  In  this  study,  we  perform  observational  analysis  to  better 
understand  the  scales  of  turbulence  in  various  environmental  conditions.  We  then  explore 
the  capability  of  the  U.S.  Navy’s  current  Coupled  Ocean-atmosphere  Mesoscale 
Prediction  System  (COAMPS™),  using  a  case  study  of  boundary  layer  roll  vortices  in 
cold  air  outbreak  conditions  over  the  Japan/East  Sea.  We  analyze  the  model’s  ability  to 
perform  accurately  within  the  resolvable  scales  and  to  accurately  represent  the  boundary 
layer  turbulence  mixing  when  the  resolution  is  high.  Based  on  results  from  the  above 
analysis,  we  modify  the  existing  parameterizations  towards  more  realistic  representations 
of  the  turbulent  processes  over  a  relatively  wide  range  of  grid  resolution,  and  test  this 
modification  within  COAMPS™. 
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I.  INTRODUCTION 


A.  MESOSCALE  MODELING  TODAY 

There  are  three  major  requirements  for  improved  numerical  weather  prediction 
(NWP):  better  model  physics,  better  observational  data,  and  better  methods  for  data 
assimilation.  These  improvements  are  very  computationally  intensive,  and  thus  advances 
in  computer  power,  coupled  with  the  trend  toward  local  modeling  efforts,  has  allowed  for 
concentrated  study  of  both  historical  and  local  real-time  mesoscale  structures  and 
dynamics,  resulting  in  extensive  evaluation,  optimization,  and  improvement  in  these  three 
key  areas  of  NWP  that  continues  today  (Kalnay  et  al.  1998). 

Kalnay  et  al.  (1998)  trace  the  history  and  improvements  of  operational  NWP  at 
the  National  Center  for  Environmental  Prediction  (NCEP),  a  review  which  is 
representative  of  improvements  at  all  major  NWP  operational  centers.  NWP  has  evolved 
from  the  381  km  resolution  of  the  National  Meteorology  Center  (NMC)  1  level 
baro tropic  model  of  the  1950’s,  running  on  a  then  state-of-the-art  IBM  704 
supercomputer,  to  NCEP’s  operational  Eta  mesoscale  model,  with  12  km  resolution  and 
45  vertical  levels,  running  on  today’s  state-of-the  art  massively  parallel  computer  system. 

While  NWP  efforts  and  advancements  in  mesoscale  modeling  continue  at  the 
national  level  (with  NCEP’s  Eta  mesoscale  model  and  the  U.S.  Navy’s  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System  (COAMPS™)  a  transition  is  occurring  where 
mesoscale  NWP  has  spread  to  the  regional  level  (Mass  and  Kuo  1998).  Three  major 
advancements,  namely,  the  advent  of  modestly  priced  single  and  multi-processor  work 
stations,  the  availability  of  mesoscale  models  that  run  efficiently  on  workstations,  and  the 
accessibility  of  initial  and  boundary  conditions  and  forecast  fields  from  operational 
centers  via  the  internet  (Mass  and  Kuo  1998),  have  given  meteorologists  the  ability  to 
locally  run  and  utilize  the  output  of  non-hydrostatic,  three-dimensional  atmospheric 
numerical  models  over  scales  of  thousands  of  kilometers  for  research,  teaching,  and 
operational  use.  As  of  early  2002,  over  fifty  state-of-the-art  mesoscale  and  regional  real- 
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time  NWP  efforts  exist  at  universities,  government  laboratories,  and  researeh  facilities  in 
the  United  States. 

Computational  improvements  and  the  vast  application  of  mesoscale  models  for 
various  meteorological  phenomena  have  allowed  for  development  of  more  detailed  model 
physics  and  parameterization  schemes  within  the  research  community  for  use  in 
mesoscale  models.  One  example  is  in  precipitation  predictions.  Precipitation  is  one  of 
the  most  difficult  aspects  of  NWP,  since  it  is  often  organized  at  the  mesoscales  between 
a  few  tens  to  a  few  hundred  kilometers.  One  example  of  work  to  improve  precipitation 
forecasts  is  the  use  of  complex  land-surface  schemes,  which  have  seen  limited  use  in 
mesoscale  modeling  due  to  limitations  in  computing  power.  While  they  have  been 
proven  important  in  climate  change  studies  using  general  circulation  models  (GCM), 
Wen  et  al.  (2000)  examine  the  role  that  land  surface  schemes  play  in  short-range 
precipitation  forecasts  within  mesoscale  models.  Two  different  land  surface  schemes  are 
compared,  providing  lower  boundary  conditions  to  the  Mesoscale  Compressibility 
Community  Model  (MC2);  the  force-restore  method,  a  simple  land  surface  scheme  used 
in  many  GCM  and  atmospheric  models  that  contains  a  two-layer  closure  condition  on  the 
heat  and  moisture  balances  at  the  air-soil  interface,  and  the  Canadian  Land  Surface 
Scheme  (CLASS),  a  more  complex  surface-vegetation-atmosphere  transfer  scheme, 
which  is  a  one-dimensional  column  model  designed  to  represent  the  average 
characteristics  of  a  horizontal  grid  square.  The  two  schemes  produced  different 
partitioning  of  the  fluxes  and  different  precipitation  distributions,  and  the  simulated 
precipitation  was  closer  to  that  observed  using  the  more  complex  CLASS  scheme.  The 
results  of  Wen  et  al.  (2000)  suggest  that  the  impact  of  land  surface  schemes  can  be 
significant  in  short-range,  high-resolution  precipitation  simulations  with  complicated 
vegetation  variations. 

Lynn  et  al.  (2001)  use  the  fifth-generation  Pennsylvania  State  University-National 
Center  for  Atmospheric  Research  Mesoscale  Model  (MM5),  at  nested  resolutions  of  15 
and  5  km  to  test  the  effects  of  a  sophisticated  land-surface  model  (the  Parameterization 
for  Land- Atmosphere  Convective  Exchange  (PLACE))  coupled  to  a  1.5  order  TKE 
closure  model  (Gayno  1994;  Shafran  et  al.  2000)  for  the  simulation  of  an  observed  case 
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of  widespread  summertime  moist  eonvection  over  the  Florida  Peninsula.  The  PLACE 
model  represents  soil  and  surface  hydrology  in  considerably  more  detail  compared  to 
force-restore,  allowing  soil  moisture  to  change  with  time  as  it  interacts  with  evaporation, 
rainfall,  runoff,  and  soil  water  drainage  to  bedrock.  The  simulations  using  the  1 .5  order 
TKE  closure  and  the  PEACE  scheme  produced  the  best  overall  skill  in  terms  of  biases  of 
the  surface  variables,  rainfall,  and  percent  cloud  cover,  and  correctly  simulated  an  early 
isolated  convective  storm  which  did  not  develop  in  the  control  run,  which  used  the 
standard  MM5  land-surface  scheme,  the  force-restore  method  (described  earlier),  coupled 
to  a  simpler  first-order  closure  BE  scheme.  Eynn  et  al.  (2001)  conclude  that  the  MM5 
model  required  the  combined,  synergistic  effect  of  the  more  detailed  land-surface  and 
turbulence  schemes  to  produce  the  most  realistic  simulation  of  the  convective  event. 

The  parameterization  of  surface  fluxes  of  heat,  momentum,  and  moisture,  such  as 
the  Louis  scheme  (Louis  1979),  has  been  used  extensively  in  global  models  as  a  simple, 
computationally  efficient  method,  and  is  also  commonly  used  within  today’s  mesoscale 
models.  The  approach  has  been  improved  in  recent  years,  most  notably  in  the 
incorporation  of  different  physical  processes  that  have  previously  not  been  included  in 
the  calculation  of  the  exchange  coefficients.  Eor  example,  the  COARE  algorithm  (Eairall 
et  al.  1996)  gives  special  considerations  for  light  wind  conditions,  corrections  for  sea 
surface  temperature,  effects  of  precipitation,  and  differences  in  roughness  lengths  for 
momentum,  sensible,  and  latent  heat  fluxes  in  the  calculation  of  the  exchange 
coefficients.  The  COARE  Algorithm,  along  with  the  Eouis  scheme,  is  an  additional 
option  for  the  surface  flux  parameterization  in  the  Navy’s  CO  AMPS™  model.  Andreas 
and  Emanuel  (2001)  review  recent  numerical  experiments  of  tropical  cyclones  where  the 
effect  of  sea  spray  was  included  within  the  surface  flux  parameterization  scheme  for  the 
high  wind  speed  regime.  Most  of  the  tests  concluded  that  adding  the  effects  of  sea  spray 
affect  the  structure  of  the  marine  BE  and  the  rate  of  tropical  cyclone  intensification. 
Andreas  and  Emanuel  (2001)  further  conclude  that  the  parameterization  of  sea  spray 
effects  must  also  consider  the  sensible  heat  flux  given  up  by  re-entrant  spray  droplets  to 
more  correctly  estimate  their  effect  on  the  exchange  coefficients. 
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New  and  more  eomplex  boundary  layer  parameterization  sehemes  have  also  been 
developed  in  reeent  years.  BL  parameterizations  used  in  mesoseale  models  have  evolved 
from  the  simpler  parameterizations  of  BL  turbulenee  using  first  order,  loeal  elosure 
methods  sueh  as  K-theory  (also  ealled  gradient  transport  theory).  First  order  closure 
based  on  K-theory,  and  local  closure  schemes  in  general,  are  described  and  evaluated  in 
some  detail  in  Stull  (1988)  and  Holstag  and  Boville  (1993). 

First-order  closure  gives  no  information  on  turbulence  intensity  or  temperature 
variance,  and  has  difficulty  with  well-mixed  layers  that  have  zero  gradients  of  mean 
variables.  Hence,  higher-order  closure  techniques  are  more  advantageous  for  estimating 
flow  within  the  BL.  In  recent  years,  BL  parameterization  schemes  using  higher-order,  as 
well  as  non-local,  closure  methods  have  been  developed  for  use  within  mesoseale 
models.  However,  with  the  increase  in  computational  requirements  and  cost  associated 
with  these  more  complex  schemes,  today’s  mesoseale  models  commonly  use  schemes 
that  combine  higher  order  (typically  1.5  order)  with  local  closure,  or  first  order  with  non¬ 
local  closure.  For  example,  the  scheme  employed  in  COAMPS™  is  a  1.5  order  closure 
model  based  on  the  Mellor-Yamada  level-2.5  model  (Mellor  and  Yamada  1982)  that 
includes  a  prognostic  equation  for  turbulent  kinetic  energy.  The  fluxes  are  derived  from 
a  local-K  approach,  but  unlike  the  formulations  for  the  eddy  diffusivity  coefficients  in  the 
earlier  schemes,  in  which  the  coefficients  are  functions  of  the  bulk  Richardson  number, 
the  eddy  diffusivity  in  the  COAMPS^"^  scheme  is  given  a  complex  algebraic  function 
involving  the  predicted  mean  and  turbulent  variables.  Another  more  complex  BL 
parameterization  scheme  is  the  Medium-range  Forecasting  (MRF)  PBL  scheme  (Hong 
and  Pan  1996),  which  is  a  first-order,  non-local  scheme  where  the  tendencies  are 
dependent  on  the  bulk  characteristics  of  the  PBL  and  include  counter-gradient  transports 
of  temperature  and  moisture  that  account  for  contributions  from  large  scale  eddies.  The 
eddy  diffusivity  coefficients  in  the  MRF  PBL  scheme  are  functions  of  the  friction 
velocity  and  the  PBL  height. 

Some  testing  of  the  sensitivity  of  these  different  boundary  layer  parameterization 
schemes  has  taken  place.  For  example  in  Braun  and  Tao  (2000)  in  the  MM5  simulation 
of  Hurricane  Bob  (1991)  at  high  resolution  (4  km),  the  model  was  able  to  produce,  to 
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varying  degrees,  the  traek  and  intensity  of  the  hurrieane,  but  results  exhibited  strong 
sensitivity  among  different  parameterization  sehemes.  Simulated  mean  sea-level 
pressures  and  maximum  winds  varied  by  about  16  mb  and  15  m  s'\  respeetively,  with  the 
Burk-Thompson  seheme  produeing  the  strongest  storms,  and  the  MRF  PBL  seheme 
produeing  the  weakest  storm.  Simulated  horizontal  preeipitation  struetures  varied 
substantially  between  the  different  sehemes  as  well,  suggesting  that  aeeurate  foreeasts  of 
preeipitation  in  hurrieanes  ean  be  just  as  sensitive  to  the  formulation  of  the  PBL 
parameterization  as  they  are  to  the  eloud  mierophysieal  parameterizations.  Bright  and 
Mullen  (2002)  evaluate  the  sensitivity  within  MM5  of  four  different  BL  sehemes  for  the 
ease  of  monsoon  eonveetion  over  Arizona,  and  eonelude  that  the  non-loeal  sehemes 
tested  more  eorreetly  prediet  the  development  of  deep,  monsoon  BL  and  assoeiated 
eonveetion. 

Cloud  parameterizations  have  also  beeome  more  eomplex,  partieularly  within  the 
framework  of  mesoseale  modeling.  A  traditional  approaeh  for  eumulus  parameterization, 
sueh  as  the  widely  used  Kain-Fritseh  seheme  (Kain  and  Fritseh  1990),  is  appropriate  for 
large-seale  models  at  eonveetively  unstable  grid  points  for  the  ereation  of  sub-grid  seale 
(SGS)  implieit  elouds  (Molinari  and  Dudek  1992).  However,  the  traditional  approaeh  has 
been  proven  to  be  less  well-posed  within  mesoseale  models,  or,  for  that  matter,  any 
model  using  grid  spaeing  near  50  km  and  below  (Cotton  and  Anthes  1989).  Studies 
addressing  this  issue  will  be  diseussed  in  a  later  seetion  eoneeming  model 
parameterizations  in  high-resolution  models.  In  addition  to  improvements  in  eloud 
parameterizations  based  on  resolution,  Deng  et  al.  (2003)  have  developed  a  shallow- 
eonveetion  parameterization  for  use  in  mesoseale  models.  The  seheme,  whieh  is  based 
on  the  Kain-Fritseh  seheme,  is  elosely  assoeiated  with  BL  turbulenee  proeesses,  and  ean 
transition  to  either  a  deep-eonveetion  seheme  in  eonditionally  unstable  environments,  or 
to  an  explieit  (resolved  seale)  seheme  in  moist  stable  environments.  Results  from  initial 
applieations  of  the  shallow-eonveetion  seheme  within  the  MM5  model,  in  both  marine 
and  eontinental  environments,  are  eonsistent  with  observed  mesoseale  thermodynamie 
struetures  and  loeal  eloud-field  parameters. 
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The  most  prominent  improvement  in  mesoscale  modeling  is  the  inereased  ability 
to  observe  the  atmosphere,  both  remotely  and  in-situ,  and  the  effective  incorporation  of 
observations  into  model  forecast  through  improved  data  assimilation  techniques. 
Geostationary  visible,  IR,  and  microwave  satellite  observations  are  now  routinely  taken 
at  1-10  km  spatial  resolution,  while  geostationary  and  polar  orbiting  satellite  systems  use 
high-resolution  spectral  techniques  to  obtain  remotely  sensed  high  vertical  resolution 
soundings.  Polar  orb  iter  high-resolution  (-'10  km)  microwave  scatterometers  measure 
sea  surface  wind  speed  and  direction.  Space-based  observations  are  enhanced  with  in- 
situ  observations  from  various  types  of  instruments,  including  radiosonde,  dropwinsonde, 
wind  profiler,  and  aircraft.  Improvement  of  observation  resolution  in  continental  regions 
to  the  sub- 10  km  range  comes  via  the  addition  of  radar  observations  (the  National 
Weather  Service’s  Next-Generation  Radar  system)  and  mesonet  observing  systems. 
Internet  capability  makes  all  of  this  data  readily  available  for  use  in  regional  mesoscale 
modeling  efforts.  Increased  computational  ability  has  also  led  to  implementation  and 
subsequent  improvements  within  numerical  models  of  three-dimensional  and,  most 
recently,  four-dimensional  variational  data  assimilation  systems,  to  replace  previously 
used  optimum  interpolation  (01)  schemes  (Parrish  and  Berber  1992,  Rabier  et  al.  2000, 
Lorenc  et  al.  2000).  Studies  have  shown  that  the  variational  data  assimilation  method 
creates  a  more  accurate  starting  point  for  NWP  integration,  and  is  more  flexible  than  the 
01  schemes,  making  it  easier  to  incorporate  new  data  sources  and  utilizing  existing  data 
better  (Barker  et  al.  2003,  Rogers  et  al.  2001). 

In  addition  to  improvements  discussed  above,  horizontal  grid  resolution  in  NWP 
has  also  increased  to  add  another  dimension  to  the  potential  of  improved  forecasting  in 
mesoscale  NWP.  The  Navy’s  COAMPS™  is  being  applied  operationally  at  9  km,  and, 
as  discussed  in  Mass  and  Kuo  (1998),  several  groups  in  the  research  community  are 
experimenting  with  high-resolution  NWP,  with  a  number  of  efforts  applying  resolutions 
of  15  km  or  less.  The  NCEP  Eta  model  will  go  into  single  digit  resolution  in  the  near 
future,  with  the  recent  acquisition  of  a  massively  parallel  computer  system,  and  its 
successor,  the  Weather  Research  and  Eorecasting  Model  (WRE)  is  being  designed  for 
resolutions  between  1  and  10  km  (Mass  et  al  2002).  NCEP  is  also  developing  a  Non- 
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hydrostatic  Mesoscale  Model  (NMM)  (Janjic  et  al.  2001),  formulated  speeifieally  for 
high-resolution  simulations,  to  better  represent  proeesses  for  whieh  the  hydrostatie 
approximation  ceases  to  be  valid.  In  addition,  private  seetor  eompanies  are 
experimenting  with  high-resolution  foreeast  systems  to  provide  foreeasts,  at  resolutions 
below  10  km,  for  major  urban  areas  of  the  United  States  (Mass  et  al.  2002). 

B.  THE  IMPORTANCE  OF  HIGH  RESOLUTION  MESOSCALE 

MODELING  FOR  NAVY  OPERATIONS 

Mesoseale  NWP  models  are  used  in  the  study  and  operational  foreeasting  of 
terrain-indueed  and  synoptieally-indueed  atmospherie  mesoseale  weather  systems  (Pielke 
2002).  Terrain  indueed  systems  are  defined  as  those  foreed  primarily  by  surfaee 
inhomogeneities,  sueh  as  sea  breeze  eireulations,  mountain-valley  winds,  foreed  flow 
over  and  around  terrain,  urban  eireulations,  and  lake  effeet  preeipitation.  Synoptieally 
indueed  mesoseale  systems  are  those  foreed  primarily  by  instabilities  within  moving 
larger  seale  disturbanees,  and  inelude  fronts  and  squall  lines,  eonveetive  bands  in 
stratiform  elouds,  mesoseale  eonveetive  systems,  and  tropieal  eyclones. 

The  eapability  to  prediet  all  of  these  types  of  mesoseale  atmospherie  phenomena 
is  required  to  maximize  the  effeetiveness  of  the  U.S.  Navy’s  battle  group  and  amphibious 
ready  group  operations  abroad.  Any  of  the  above  weather  seenarios  is  a  potential 
environment  for  Navy  surveillanee,  seareh,  and  strike  operations,  as  well  as  the  medium 
in  whieh  it  must  funetion  without  damage  to  platforms  and  personnel.  U.S.  Navy 
meteorologists  supporting  fleet  operations  in  speeifie  loealized  areas  of  interest  use  the 
Navy  COAMPS^'^  atmospherie  model  as  a  primary  tool  to  support  and  to  maximize  the 
effeetiveness  of  operations.  Aeeurate,  high-resolution  atmospherie  model  output  would 
be  extremely  useful  to  the  Navy  foreeaster  providing  guidanee  in  fleet  operations.  Being 
able  to  aeeurately  resolve  the  detail  of  flow  fields  around  eomplex  terrain  to  pinpoint 
regions  of  strong  or  turbulent  winds,  or  to  aeeurately  foreeast  the  loeation  of  strong 
eonveetive  aetivity  near  a  eold  front  or  mesoseale  eonveetive  system,  for  example,  would 
play  a  key  role  in  the  overall  sueeess  of  a  partieular  mission. 

One  speeifie  example  where  weather  played  a  signifieant  role  in  strike  operations 
was  during  the  air  eampaigns  of  the  desert  eonfliets  in  the  Middle  East.  Aireraft  without 
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target  radars  are  especially  dependent  on  clear  weather  conditions  for  effective 
operations.  For  example,  dust  storms  and  low  cloud  cover  were  weather  events  that 
significantly  affected  air  operations.  These  phenomena  are  usually  short-lived  and 
confined  to  specific  areas,  which  can  leave  one  airfield  or  target  area  in  the  dark  or  under 
cloud  cover  for  hours,  severely  hampering  visibility,  while  at  the  same  time  conditions 
are  clear  at  another  airfield  or  target  20  miles  away.  The  development  of  high-resolution 
mesoscale  models,  that  could  accurately  resolve  flow  fields  affected  by  varying  local 
terrain  structure  and  coastal  topography,  are  important  tools  in  predicting  these  localized 
weather  events  and  providing  short-term  forecasts  for  military  flight  operations. 

As  mesoscale  model  resolutions  increase  to  the  sub-km  scales,  it  will  be  possible 
to  resolve  features  at  the  turbulence  scale  such  as  BL  roll  vortices  and  clear-air  turbulent 
eddies.  While  sub-km  resolution  modeling  is  not  presently  possible  operationally, 
researchers  are  discovering  its  potential  today.  For  example,  Burk  and  Haack  (2000)  use 
COAMPS™  to  model  coastally-trapped  wave  clouds  off  the  central  California  coastline 
using  a  horizontal  resolution  of  1/3  km,  and  present  the  highest  real  data  forecast  using 
COAMPS™  to  date.  The  simulation  accurately  predicted  long,  lineal  variations  in  the 
wind,  temperature,  and  moisture  fields,  with  modeled  wavelengths  ~  4km,  consistent 
with  cloud  features  observed  in  satellite  imagery.  The  ability  to  predict  the  development 
of  these  turbulent  features,  which  have  vertical  updrafts-downdrafts  of  up  to  4  ms'^, 
would  offer  considerable  improvement  to  turbulence  forecasting  for  low-level  flight 
operations.  The  effort  to  improve  short  term  forecast  accuracy  as  operational  mesoscale 
models  attain  increasingly  higher  resolution  is  critical  to  the  Navy’s  mission. 


C.  PREVIOUS  STUDIES  IN  HIGH  RESOLUTION  MESOSCALE 
MODELING 

Many  previous  studies  have  examined  the  effects  of  horizontal  resolution  on  the 
accuracy  of  mesoscale  NWP  forecasts  (Mass  et.  al  2002).  Some  examples  of  these 
studies,  modeling  different  mesoscale  weather  phenomena  at  high  horizontal  grid 
resolution,  are  reviewed  below.  In  most  cases,  subjective  evaluation  shows  that 
increasing  the  horizontal  resolution  below  10-15  km  produces  more  realistic,  better- 
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defined  strueture  of  the  atmospherie  features  being  modeled.  On  the  other  hand,  few 
studies  have  demonstrated  that,  objeetively  measured  over  time,  foreeast  aeeuracy 
inereases  as  grid  spaeing  deereases  below  10-15  km  (Mass  et  al.  2002). 

Mesoscale  models  are  approaehing  the  resolution  neeessary  to  resolve  key 
features  in  regions  of  complex  terrain  that  have  a  direct  and  significant  impact  on 
precipitation.  Several  studies  (Bruintjes  et  al.  1995;  Colle  and  Mass  1996;  Gaudet  and 
Cotton  1998;  and  Westrick  and  Mass  2001)  have  shown  that,  when  run  at  sufficiently 
high  resolution  (down  to  '-lO  km),  mesoscale  models  can  realistically  simulate  observed 
precipitation  structures  over  complex  terrain.  Simulating  gap  flow  through  the  Strait  of 
Juan  de  Fuca  using  MM5,  Colle  and  Mass  (2000a)  found  that  the  sharp  transition  of  the 
Strait-exit  became  more  realistic  as  horizontal  grid  spacing  was  decreased  from  12  to 
1.33  km,  but  the  simulation  underestimated  the  winds  by  2-5  ms'^  in  the  gap  exit  region 
due  to  the  30%  underestimation  of  the  sea  level  pressure  gradient  down  the  strait.  Colle 
and  Mass  (2000b)  found  significant  improvement  in  precipitation  forecast  skill  as  grid 
spacing  decreased  from  36  to  4  km.  For  most  locations,  there  was  little  overall 
improvement  going  from  4  to  1.33  km  except  along  the  upper  windward  slopes  and  in  the 
immediate  lee  of  major  barriers,  where  precipitation  amounts  were  enhanced  with  higher 
resolution. 

Non-hydrostatic  models  employing  grid  spacings  of  4  km  or  less  have  shown 
substantial  success  in  duplicating  observed  structural  evolution  of  a  range  of  mesoscale 
convective  systems  (e.g.,  Weisman  et  al.  1988;  Skamarock  et  al.  1994;  Droegemeier  et  al. 
1994).  For  example,  Bemadet  et  al.  (2000),  simulating  four  high  plains  convective 
events  from  approximately  80-  to  2-km  grid  spacing  using  the  Colorado  State  Regional 
Atmospheric  Modeling  System  (RAMS)  model,  found  that  2  km  was  necessary  to 
capture  convection  explicitly.  In  the  simulation  of  Hurricane  Danny,  Kuo  et  al.  (2001) 
found  substantial  improvement  in  the  radius  of  maximum  wind  and  the  eye  wall  and  rain 
band  structures  as  horizontal  grid  spacing  was  reduced  from  81  to  1  km.  The  3  km 
resolution  simulation  produced  a  radius  of  maximum  wind  at  15  km  consistent  with  radar 
observations,  with  results  further  improved  using  1  km  resolution.  The  study  suggested 
that  a  grid  resolution  of  at  least  3  km  was  necessary  to  capture  the  wind  field  structure  of 
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this  small,  minimal  hurricane,  and  further  suggests  the  use  of  1  km  resolution  or  higher  to 
properly  simulate  the  evolution  of  mesoseale  eonveetive  systems  and  their  interaetion 
with  the  hurrieane  vortex. 

Mid-latitude  eyclones  and  frontal  zones  have  been  simulated  using  mesoseale 
NWP  models  at  high  resolution  in  several  reeent  studies.  Simulating  a  Paeific  cold  front 
with  MM5,  Chien  et  al.  (2001)  found  that  decreasing  the  grid  spaeing  from  45  to  15  km 
did  not  greatly  ehange  the  system’s  synoptie  structure;  however,  deereasing  the  grid 
spacing  to  1.67  km  was  necessary  to  produee  the  observed  narrow  frontal  structure, 
assoeiated  eonveetive  rain  band,  and  gravity-eurrent  head-like  strueture  near  the  front’s 
leading  edge.  In  an  MM5  simulation  of  a  polar  low  that  developed  over  the  Bering  Sea, 
Breseh  et  al.  (1997)  found  that  although  a  grid  spaeing  of  6.7  km  produeed  sharper 
mesoseale  features  than  a  20  km  resolution  grid,  the  eentral  pressure  and  synoptie  seale 
structures  were  nearly  the  same  at  both  resolutions.  Buckley  and  Leslie  (2000)  found  that 
deereasing  grid  spaeing  from  50  to  10  km  in  the  University  of  New  South  Wales  High 
Resolution  Limited-Area  Atmospherie  Predietion  (HIRES)  model  greatly  improved  the 
wind  and  central  pressure  foreeasts  for  the  Boxing  Day  storm  of  Deeember  1998.  In 
general,  these  and  other  subjeetive  studies  suggest  that  moving  to  high  resolution  (grid 
spaeing  less  than  15  km)  does  produee  better  definition  of  frontal  zones  and  other 
mesoseale  weather  features,  but  does  not  signifieantly  improve  the  synoptie  evolution  of 
fronts  and  eyelones  (Mass  et  al.  2002). 


D,  HIGH  RESOLUTION  AND  MODEL  PARAMETERIZATIONS 

As  stated  earlier.  Cotton  and  Anthes  (1989)  suggest  that  the  eoneeptual  basis  for 
eumulus  parameterization  beeomes  unelear  when  model  grid  spaeing  falls  below  about 
50  km.  Along  these  lines,  some  researehers  have  eliminated  eumulus  eonveetion  in  high- 
resolution  models,  and  explieitly  simulate  eumulus  eonveetion  at  the  grid  scale  (e.g., 
Yamasaki  1977;  Rosenthal  1978).  Other  researehers  have  designed  eumulus 
parameterizations  with  elosure  assumptions  more  appropriate  for  grid  resolutions  below 
50  km  (Fritseh  and  Chappell  1980;  Frank  and  Cohen  1987).  At  resolutions  typieal  of 

today’s  mesoseale  models,  on  the  order  of  tens  of  kilometers,  the  two  sehemes  must  be 
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used  concurrently:  an  implicit  parameterization  scheme  for  the  SGS  convection,  and  an 
explicit  scheme  for  modeling  for  larger-scale  clouds  resolved  on  the  grid  scale.  Belair  and 
Mailhot  (2001)  examine  the  relative  roles  of  implicit  and  explicit  convection  schemes  in 
the  numerical  simulation  of  a  MCS  and  squall  line  using  the  MC2  mesoscale  model  at  2, 
6,  18,  and  50  km  resolutions.  At  2  km  resolution,  the  explicit  scheme  is  used  alone,  and 
the  model  realistically  reproduces  both  the  mesoscale  and  small-scale  features  of  the 
system  without  the  use  of  any  parameterization.  In  the  6  km  simulation,  where  the 
implicit  scheme  is  employed  for  the  SGS  convection,  both  the  convective  and  anvil 
regions  of  the  MCS  are  still  resolved  by  the  model,  but  the  leading  convective  band  is 
less  realistic  (wider,  shallower,  and  less  intense)  than  in  the  explicit  2  km  simulation. 
Belair  and  Mailhot  (2001)  explain  that  at  6  km  resolution,  the  grid  spacing  is  so  close  to 
the  convective  scale  that  the  leading  convective  band  of  the  squall  system  is  barely 
resolved  by  the  model,  and  is  thus  the  implicit  and  explicit  schemes  coincidentally 
parameterize  part  of  the  convection.  As  previously  emphasized  in  Frank  (1983)  and 
Molinari  and  Dudek  (1992),  the  use  of  implicit  schemes  at  resolutions  between  3  and  10 
km  is  in  violation  of  the  parameterization’s  fundamental  assumptions  and  closure 
hypothesis.  Belair  and  Mailhot  (2001)  state  that  because  of  the  lack  of  separation  between 
the  convective  and  grid  scales,  results  obtained  from  models  at  these  intermediate 
resolutions  must  be  examined  cautiously.  Belair  and  Mailhot  (2001)  conclude  that  the 
necessity  to  parameterize  convection,  even  at  this  high  resolution,  should  encourage 
researchers  to  find  ways  to  represent  the  effects  of  unresolved  convection  in  the  gap 
between  the  scales  of  direct  simulation  of  clouds  (~100  m)  (Arakawa  and  Chen  1987)  and 
•-lO  km,  where  explicit  modeling  of  clouds  is  insufficient  to  represent  the  convection 
realistically. 

The  treatment  of  SGS  turbulence,  similar  to  that  of  convection,  becomes  less  well 
posed  as  the  horizontal  resolution  increases.  As  resolution  increases  into  the  turbulent 
scale,  the  appropriateness  of  the  turbulence  parameterization  scheme  itself  now  becomes 
questionable,  and  the  need  to  re-evaluate  the  appropriateness  of  turbulence 
parameterizations  arises.  Belair  et  al.  (1998)  addressed  the  issue  of  approaching  the 
turbulence  scale  in  high-resolution  numerical  modeling.  They  assessed  the  impact  of 
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increasing  the  horizontal  resolution  from  10  km  to  1  km  on  the  simulation  (using  the 
Freneh  MESO-NH  limited  area  model)  of  surfaee  and  turbulent  fluxes  in  a  field 
experiment  over  cultivated  land  and  forest  regions  in  southwestern  Franee.  Comparing 
the  model  output  with  aireraft  measurements,  they  found  that  the  inereased  resolution  to  1 
km  showed  mueh  smaller  seale  detail  in  the  parameterized  fluxes  and  turbulent  kinetie 
energy,  but  there  were  oseillations  in  the  parameterized  quantities  that  did  not  appear 
reasonable.  They  eonelude  that  sinee  a  non-negligible  portion  of  the  turbulent  fluxes  are 
teehnieally  resolvable  by  the  1  km  model,  this  portion  of  the  energy  should  not  be 
parameterized  by  the  turbulent  seheme,  which  is  designed  to  parameterize  the  vertieal 
mixing  based  on  the  total  ensemble  of  turbulent  energy.  Their  study  thus  indieates  that  at 
intermediate  resolutions  (between  200  m  and  a  few  kilometers)  the  interpretation  of 
model  output  relating  to  turbulent  quantities  is  eomplieated,  and  that  the  higher  resolution 
does  not  neeessarily  improve  the  simulation  of  the  turbulenee.  This  raises  the  questions: 
at  what  resolution  should  the  validity  of  SGS  parameterization  be  questioned,  and  how 
ean  we  properly  parameterize  the  SGS  proeesses  when  the  model  resolution  is  elose  to 
the  seale  of  turbulent  eddies?  These  questions  have  no  known  published  answers,  and 
are  key  issues  of  this  study. 


E.  OBJECTIVES  OF  RESEARCH 

While  the  Navy  operational  COAMPS™  is  eurrently  using  an  81-27-9  km  nest, 
the  inevitable  inerease  in  available  eomputing  power,  both  at  Navy  METOC  eenters  and 
afloat,  will  allow  even  higher  resolution  model  runs  to  beeome  operational.  With  this 
eventuality  and  the  discussion  in  the  previous  seetion  in  mind,  the  main  goal  of  this  study 
is  to  answer  the  basie  questions  and  better  understand  the  issues  relating  to  the  grid 
resolution  dependenee  of  turbulenee  parameterizations  in  mesoseale  models  where  the 
resolution  approaehes  the  seale  of  turbulent  eddies.  Understanding  the  seale  of 
turbulenee,  and  the  factors  affecting  it,  is  also  an  important  objective  in  determining  in 
the  appropriateness  of  ensemble  turbulence  parameterization  in  a  high-resolution 
numerieal  simulation  when  the  grid  resolution  approaehes  the  seale  of  turbulent  eddies. 
Through  analyses  of  a  numerieal  simulation  using  a  ease  study  of  eold  air  outbreak 
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conditions  over  the  Japan/East  Sea  (JES),  we  explore  the  capability  of  the  eurrent 
COAMPS™  turbulenee  parameterization  scheme,  at  high  resolution,  in  simulating  the 
boundary  layer  turbulenee  mixing.  Results  from  these  simulations  are  compared  with 
those  from  direet  measurements.  Based  on  results  from  the  above  analysis,  we  attempt  to 
modify  the  existing  COAMPS™  turbulenee  parameterization  towards  more  realistie 
representations  of  the  turbulent  proeesses.  The  modifieation  will  be  applied  within 
COAMPS™  for  the  same  ease  study,  and  results  of  the  eontrol  and  modified  simulations 
will  be  eompared  and  verified  against  observations  to  determine  if  improvement  was 
indeed  made. 

F.  SUMMARY  OF  DISSERTATION 

Chapter  II  eontains  baekground  information  on  turbulent  seales,  the  ensemble 
turbulenee  parameterization  seheme  used  within  COAMPS™,  and  mixing  length 
formulations  used  within  different  BE  turbulenee  parameterization  sehemes.  Chapter  III 
eontains  an  observational  analysis  of  the  seale  of  turbulenee  and  the  integral  length  seale. 
The  analysis,  whieh  utilizes  aircraft  observations  from  three  different  eonvective  BE  field 
experiments,  provides  the  coneeptual  basis  for  modifieations  made  to  the  COAMPS™ 
mixing  length  used  in  the  turbulenee  parameterization  scheme.  In  ehapter  IV,  an  in-depth 
analysis  of  the  JES  eold  air  outbreak  case  study  is  performed,  where  the  BE  strueture 
from  aircraft  observations  is  eompared  with  that  simulated  in  the  COAMPS™  eontrol 
run.  It  is  determined  that  the  model  poorly  simulates  the  observed  BE  turbulent  strueture, 
in  terms  of  the  spatial  variability  near  the  minimum  resolvable  seale.  Chapter  IV  also 
includes  a  baekground  discussion  of  horizontal  roll  vortices,  whieh  are  observed  in  the 
BE  for  this  ease  study.  Chapter  V  explains  modifieations  made  to  the  COAMPS™ 
mixing  length,  whieh  is  applied  within  the  model  for  the  JES  cold  air  outbreak  ease  study 
in  an  effort  to  improve  the  simulated  spatial  variability.  The  results  from  the  simulation 
with  modified  mixing  length  are  eompared  to  that  of  the  eontrol  run  to  determine  if 
improvement  in  the  simulated  BE  strueture  is  aehieved.  Chapter  VI  provides  a  projeet 
summary  and  eonelusions,  and  reeommendations  for  further  study  of  this  issue. 
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II.  BACKGROUND 


A.  PHYSICAL  PROCESSES  GENERATING  BOUNDARY  LAYER 

TURBULENT  KINETIC  ENERGY 

Modeling  of  BL  turbulence  in  COAMPS^”^  is  the  main  focus  of  this  research,  thus 
it  is  important  to  understand  the  various  physical  processes  that  act  to  generate  and 
maintain  turbulence  there.  The  BL  is  that  layer  of  the  atmosphere  directly  influenced  by 
the  earth’s  surface,  responding  to  surface  forcings  on  a  time  scale  of  1  hour  or  less,  and 
contains  a  significant  amount  of  turbulent  energy  (Stull  1988).  The  height  of  the  BL  top 
is  dependent  on  the  balance  between  the  strength  of  the  subsidence  in  the  free  atmosphere 
and  the  strength  of  turbulence  within  the  BL.  The  vertical  motions  generated  by 
turbulence  act  to  mix  heat,  moisture,  and  momentum  fluxes  generated  at  the  surface 
upward. 

Turbulent  kinetic  energy  (TKE)  is  a  measure  of  the  intensity  of  turbulence  in  the 
BL,  and  is  directly  related  to  the  momentum,  heat,  and  moisture  transport  through  the  BL. 
The  individual  terms  of  the  TKE  budget  equation  (time  tendency  of  TKE)  describe  the 
physical  processes  that  generate  turbulence.  Here  we  present  a  simplified  version  of  the 
TKE  budget  equation  for  the  purposes  of  describing  the  important  processes  which 
generate  TKE  in  the  BE.  The  simplified  TKE  budget  equation  for  a  field  with  horizontal 
homogeneity  (Stull  1988), 


de 

dt 


6^  dz  dz  p  dz 


(2.1) 


I  II  III  IV  V  VI 

assumes  a  coordinate  system  aligned  with  the  mean  wind,  no  horizontal  advection  of 
TKE,  and  neglects  subsidence.  Term  I  represents  the  local  storage  or  tendency  of  TKE. 
Turbulent  kinetic  energy  (TKE)  is  generated  through  buoyancy  production  (term  II)  and 
shear  production  (term  III),  mainly  occurring  in  the  lower  half  boundary  layer.  A 
positive  buoyancy  flux  generates  TKE  and  represents  the  effects  of  thermals  in  the 
convective  BE;  for  example,  a  warm,  sunny  day  over  land  or  the  advection  of  cold  air 
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over  a  relatively  warmer  surface.  Static  stability  tends  to  suppress  or  consume  TKE  as  it 
generates  negative  buoyancy  flux.  Such  conditions  are  present  in  the  stable  BL  at  night 
over  land,  or  when  the  surface  is  colder  than  the  overlying  air.  Shear  production  occurs 
in  the  presence  of  a  mean  wind  shear.  The  interaction  tends  to  generate  more  turbulence, 
since  the  two  terms  are  usually  of  opposite  sign,  and  thus  shear  production  is  a  positive 
term  in  the  TKE  budget  equation.  Term  IV  is  the  vertical  flux  of  TKE,  and  represents  the 
vertical  transport  of  TKE  by  turbulent  eddies.  This  term  is  not  a  production  or  loss  term 
(its  value  is  zero  when  integrated  of  over  the  depth  of  the  BE),  but  acts  to  redistribute 
TKE  within  the  BE.  Term  IV  is  generally  negative  in  the  lower  half  of  the  BE,  as 
turbulence  generated  there  by  shear  and  buoyancy  production  is  transported  to  the  upper 
half  of  the  BL,  where  term  IV  is  generally  positive.  Term  V  is  a  pressure  correlation 
term,  which  describes  how  TKE  is  redistributed  in  different  directions.  Little  is  known 
about  this  term  from  observations,  and  it  is  often  calculated  as  the  residual  required  to 
balance  the  other  terms  in  the  TKE  budget.  Term  VI  is  the  dissipation  of  TKE, 
representing  the  molecular  destruction  of  turbulent  motions.  Turbulent  energy  is 
dissipated  within  the  BL  through  the  energy  cascade  process,  discussed  in  the  next 
section,  and  the  magnitude  of  this  dissipation  rate  is  proportional  to  the  rate  of  TKE 
production.  Dissipation  occurs  primarily  within  the  smallest  turbulent  eddies,  at  the 
smallest  scales,  while  production  of  TKE  occurs  in  the  largest  turbulent  eddies.  Thus,  the 
production  and  loss  mechanisms  defined  in  the  TKE  budget  occur  on  completely 
different  scales.  Eigure  2.1  is  a  figure  from  Stull  (1988)  that  shows  the  variation  of  the 
terms  with  height  within  a  fair  weather  convective  BL.  The  relative  magnitudes  of  the 
terms  shown  in  the  figure  are  derived  from  several  data  and  models  from  several  different 
experiments. 

While  not  the  focus  of  this  study,  entrainment  is  an  important  BL  process  that 
responds  to  and  feeds  back  to  turbulence  in  the  BL,  and  is  a  mechanism  for  exchange 
between  the  BL  and  the  free  atmosphere.  Entrainment  of  above  BL  air  from  the  free 
atmosphere  into  the  BL  is  the  mechanism  for  BL  vertical  growth,  and  allows  the  free 
atmosphere  to  interact  with  the  BL  and  respond  to  changes  in  surface  forcing.  In  the 
clear  BL,  turbulent  eddies  rising  through  the  BL  become  negatively  buoyant  as  they  reach 
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Figure  2,1.  Terms  in  the  TKE  budget  equation,  normalized  by  BL  height.  Shaded 
areas  indieate  ranges  of  values.  From  Stull  (1988). 


the  inversion,  but  overshoot  into  the  free  atmosphere  before  sinking  baek  into  the  BL. 
During  the  overshoot  into  the  inversion,  some  free  atmospherie  air  is  pushed  into  the  BL 
and  is  rapidly  mixed  into  the  BL  beeause  of  the  strong  turbulenee  there,  and  results  in  a 
one-way  entrainment  proeess  allowing  the  BL  to  grow  in  thiekness.  If  elouds  are  present 
in  the  BL,  warm  air  from  above  the  BL  entrained  into  the  top  of  a  eloud  ean  eool  and  sink 
through  evaporation  of  eloud  droplets.  This,  along  with  IR  eooling  at  the  eloud  top,  form 
negatively  buoyant  downdrafts  generating  more  turbulenee  that  enhanee  BL  mixing  and 
entrainment.  Beeause  the  strength  of  the  entrainment  is  direetly  proportional  to  the 
strength  of  BL  turbulenee,  ehanges  to  the  free  atmosphere  are  dependent  upon  the  extent 
of  turbulent  flow  within  the  BL. 

B,  SCALES  OF  MOTION  IN  THE  ATMOSPHERE 

1,  Atmospheric  Energy  Spectrum 

Atmospherie  flow  as  a  whole  ean  be  deseribed  as  eddies  of  different  sizes,  eaeh 
earrying  a  eertain  amount  of  kinetie  energy.  This  is  best  shown  in  the  sketeh  of  the 
energy  density  speetrum  of  wind  speed  (Figure  2.2).  Here,  the  leftmost  peak  eorresponds 
to  kinetie  energy  assoeiated  with  synoptie  scale  phenomena.  The  spectral  minimum  at 
about  1  hour  allows  for  a  separation  of  scales  between  the  synoptic  and  mesoscale 
features,  or  mean  flow,  with  microscale  features,  or  turbulence,  in  the  BL.  In  this  gap 
region  is  where  several  special  mesoscale  types  of  flow  are  known  to  exist,  such  as  deep 
convection,  large  roll  vortices,  and  local  circulations  such  as  land/sea  breezes  (Smedman 
and  Hogstrom  1975).  Because  of  these  phenomena,  the  “gap”  is  sometimes  not  well- 
defined,  and  is  likely  to  contain  different  types  of  motion,  which  fortunately  are  distinctly 
varied  in  terms  of  spectral  signature.  The  remainder  of  the  spectral  energy  to  the  right  of 
the  gap  is  associated  with  turbulent  energy  within  the  atmospheric  BL.  Note  that  the 
relative  magnitudes  of  the  synoptic  scale  and  turbulent  scale  spectral  peaks  in  Figure  2.2 
are  for  illustrative  purposes  only  and  are  not  to  scale;  synoptic  scale  energy  is  in  reality  2 
or  more  orders  of  magnitude  larger  than  turbulence  scale  energy.  The  magnitude  of  TKE 
will  be  illustrated  clearly  in  Chapter  III. 
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2,  Turbulent  Energy  Spectrum 

Turbulence  is  the  gustiness  superimposed  on  the  mean  wind  (Stull  1988)  and  can 
be  visualized  as  irregular  swirls  of  motion,  known  as  turbulent  eddies.  The  rightmost 
peak  in  the  atmospheric  energy  spectrum  of  Figure  2.2  is  associated  with  the  energy- 
containing  turbulence  eddies  within  the  atmospheric  BL.  These  eddies  are  on  the  scale  of 
the  BL  depth,  roughly  100  m  to  3  km  in  diameter  and  are  produced  directly  by  the 
dynamic  instability  of  the  mean  flow  via  shear  or  buoyancy  effects.  The  smallest  eddies 
(on  the  order  of  a  few  millimeters  in  size)  are  called  the  dissipation  eddies  because  of  the 
dissipating  effects  of  molecular  viscosity.  Eddies  that  scale  on  the  order  of  centimeters  to 
100  m,  exist  on  a  scale  known  as  the  inertial  sub-range,  and  are  produced  by  the  energy 
cascade  process  from  the  large  eddies.  There  is  a  net  transfer  of  energy  from  the  largest 
to  the  smallest  eddies,  where  large  eddies  create  eddy-size  wind  shear  regions  that 
generate  smaller  eddies,  a  cycle  which  continues  until  the  smallest  eddies  dissipate  into 
heat  via  molecular  viscosity.  After  many  cascade  processes,  the  eddies  become  so  small 
and  contain  so  little  energy  that  they  are  no  longer  considered  a  mechanism  in  the  transfer 
of  momentum  and  heat,  and  are  much  less  dependent  on  the  flow  than  large  eddies.  The 
issue  of  understanding  the  scales  of  turbulence  is  important  to  improved  high  resolution 
simulations  of  BL  turbulent  motions,  and  will  be  addressed  in  great  detail  in  Chapter  III, 
through  the  analysis  of  turbulence  spectra  from  high  rate  aircraft  observations. 

3,  Scales  of  Atmospheric  Turbulence 

Several  length  scales  have  been  used  to  denot  the  size  of  turbulent  eddies.  Some 
of  the  more  commonly  used  length  scales  are  presented  in  Kaimal  (1973).  The  integral 
length  scale  (A)  is  the  classical  length  scale  of  fluid  mechanics,  and  is  obtained  from  the 
integral  time  scale  (x)  using  Taylor’s  “frozen  field”  hypothesis  (Taylor  1938)  using 

00 

K  =  Ut  =  U^  p{t)dt  ,  (2.2) 

0 

where  p(t)  is  the  autocorrelation  coefficient.  The  wavelength  corresponding  to  the  peak 
of  the  logarithmic  power  spectrum  (km),  also  derived  using  Taylor’s  hypothesis,  is 
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Figure  2,2,  Spectra  of  wind  speed  near  the  surface,  illustrating  the  scales  of 
motion  within  the  atmosphere.  From  Stull  (1988). 
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(2.3) 


K=Ulf^, 

where  fm  is  the  peak  frequeney.  This  length  seale  is  widely  used  in  the  interpretation  of 
atmospherie  speetra,  sinee  it  ean  be  derived  direetly  from  the  logarithmie  power 
speetrum,  and  will  be  used  in  the  work  presented  here  in  our  diseussion  of  turbulent 
length  seales.  A  and  km  are  related  empirieally  in  Kaimal  (1973)  by 

4  =  InK  .  (2.4) 

Several  studies  have  addressed  the  seales  of  turbulenee  by  analyzing  the  speetra 
of  veloeity  eomponents  and  temperature  using  observations  taken  within  the  surfaee  layer 
and  outer  layer  of  the  BL.  Studies  of  this  nature  also  provide  insight  into  the  physieal 
meehanisms  that  are  aetive  in  the  turbulent  transfer  proeesses  in  the  BL.  Some  of  these 
studies,  pertinent  to  the  work  presented  here,  are  reviewed  below.  The  diseussion  will 
foeus  mainly  on  the  vertieal  veloeity  speetra  and  length  scales,  and  will  be  referenced 
later  in  the  analysis  of  vertical  velocity  spectra  from  different  experiments  over  the  ocean. 

Much  of  the  pioneering  work  in  the  study  of  turbulent  energy  spectra  utilized 
similarity  theory  (Kolmogorov  1941;  Obukov  1941)  to  normalize  spectral  intensities 
according  to  scaling  variables  appropriate  to  the  flow.  This  approach  allowed  researchers 
to  organize  results  and  generalize  the  behavior  of  turbulence.  For  example,  turbulent 
eddies  in  the  inertial  sub-range  obtain  their  energy  inertially  from  larger  eddies  and  lose 
their  energy  in  the  same  way  to  smaller  eddies.  For  a  steady  state  turbulent  flow,  the 
cascade  rate  of  energy  down  the  spectrum  must  balance  the  dissipation  rate  at  the 
smallest  eddy  sizes.  Thus,  three  variables  govern  this  process:  spectral  intensity  (S),  the 
wavenumber  (k)  and  the  dissipation  rate  (s).  Dimensional  analysis  leads  to  an  equation 
for  the  spectral  intensity  in  the  inertial  sub-range, 

S^=aa^‘^K^‘\  (2.5) 

where  a  is  the  Kolomogrov  constant,  with  a  typical  value  around  1.5  (Gossard  et.  al 
1982).  This  relationship  has  allowed  for  the  determination  of  the  inertial  sub-range 
within  any  measured  spectrum,  since  the  plot  of  S  vs.  k  on  a  log-log  scale  would  present 
a  straight  line  portion  of  slope  -5/3,  indicating  the  location  of  the  inertial  sub-range. 
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Kaimal  et.  al  (1972)  studied  the  behavior  of  turbulenee  speetra  in  the  surfaee 
layer,  defined  as  the  lowest  -10%  of  the  BL  (Stull  1988),  within  the  framework  of 
similarity  theory.  The  study  utilized  data  obtained  within  a  surfaee  layer  over  fiat 
grasslands.  The  goal  of  their  study  was  to  define  the  general  behavior  of  veloeity  and 
temperature  speetra  for  different  heights  within  the  surfaee  layer  and  under  different 
eonditions  of  stability.  In  their  approaeh  to  this  goal,  the  speetra  were  eollapsed  into 
eoineidenee  in  the  inertial  sub-range,  by  normalizing  the  veloeity  and  temperature  speetra 
to  remove  the  dependenee  on  stability  there.  The  ability  to  eollapse  the  eurves  by 
removing  the  stability  dependenee  implies  that  the  smaller  eddies  in  the  inertial  sub-range 
reeeive  all  of  their  energy  via  the  easeade  proeess  from  larger  eddies,  with  no  direet 
interaetion  with  the  mean  flow  or  mean  stratifieation.  The  result  allows  the  observation 
of  the  variation  in  the  speetra  as  a  funetion  of  stability  above  the  inertial  sub-range,  at  the 
seale  of  energy-eontaining  eddies.  The  resulting  generalized  speetral  eurves,  derived 
from  the  normalized  w  speetra,  are  plotted  for  a  range  z/L  values  from  stable  (2.0)  to 
unstable  (-2.0)  for  vertieal  veloeity  (w)  in  Figure  2.3.  The  eurves  provide  some  important 
implieations  for  the  eharaeteristies  and  the  seales  of  turbulent  flow  in  the  surfaee  layer. 

One  notable  feature  of  the  speetrum  of  the  stable  surfaee  layer  is  its  shape,  whieh 
is  insensitive  to  ehanges  in  height  and  stability.  Also,  we  see  that  the  peak  speetral 
intensity  is  redueed  as  the  stability  inereases.  This  is  a  reasonable  result,  sinee  stability 
opposes  turbulent  motions.  We  also  observe  a  shift  toward  higher  frequeneies,  or  smaller 
seales,  of  the  peak  speetral  intensity  as  stability  inereases,  possibly  due  to  stronger 
damping  by  buoyaney  forees  at  the  lower  frequeneies  (Stull  1988).  The  stippled  region 
for  the  range  -0.3  >  z/L  >  -2.0  indieates  the  absenee  of  any  well-defined  trend  or 
eonformanee  to  a  single  shape  in  the  speetra  for  the  unstable  region,  where  the  systematie 
progression  seen  on  the  stable  side  breaks  down  as  z/L  ehanges  sign.  The  generalized 
results  from  Kaimal  et  al.  (1972)  were  eonsistent  with  surfaee  layer  speetra  reported  by 
other  investigators.  Kaimal  (1973)  further  extended  the  use  of  similarity  theory  to  bring 
the  speetra  in  the  stable  surfaee  layer  into  one  universal  eurve,  by  developing  a  modified 
frequeney  seale.  In  establishing  a  relationship  between  this  new  frequeney  seale  and 
stability,  Kaimal  determined  that  the  integral  length  seale  of  eddies  in  the  stable  surfaee 
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Figure  2,3.  Generalized  surfaee  layer  vertieal  veloeity  speetrum  for  z/L  values 
ranging  from  +2.0  to  -2.0.  From  Kaimal  et.  al  (1972). 
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layer  was  inversely  proportional  to  the  Richardson  number  (Ri).  Kaimal  also  generalized 
that  the  scale  of  eddies  in  the  stable  surface  layer  becomes  more  proportional  to  height  as 
stability  decreases.  These  studies  provided  a  comprehensive  picture  of  surface  layer 
turbulent  structure  and  the  size  of  surface  layer  turbulence,  but  also  indicated  a  need  for 
similar  work  at  higher  levels  in  the  BL,  since  it  appeared  that  surface  layer  turbulence 
was  strongly  influenced  by  scales  of  motion  large  enough  to  encompass  the  whole  BL. 

Kaimal  et  al.  (1976)  were  able  to  similarly  generalize  the  spectra  in  the  convective 
mixed  layer  (ML)  using  similarity  theory.  This  experiment  utilized  observations  from 
various  levels  in  the  surface  layer  and  ML  within  a  well-mixed,  daytime  convective  BL 
over  flat  land.  They  found  that,  within  the  ML,  the  energy  in  the  inertial  sub-range 
remained  essentially  constant  with  height,  and  the  spectral  peaks  were  invariant  in 
intensity  and  location  on  the  frequency  scale.  This  near  uniform  spectral  behavior  over 
much  of  the  BL  allowed  for  the  generalization  of  the  normalized  ML  spectra  expressed  in 
mixed  layer  similarity  coordinates.  Figure  2.4  shows  the  generalized  mixed  layer  spectra 
for  the  velocity  components.  At  lower  frequencies  the  velocity  spectra  separate  as  a 
function  of  z/zi,  (where  Zi  is  the  height  BL  top).  The  large  spread  with  height  between 
vertical  velocity  spectra  implies  that  the  scales  of  vertical  velocity  are  strongly  dependent 
on  the  z/zi  compared  to  the  horizontal  velocity  components,  which  only  have  two  closely 
separated  general  ranges  of  z/zi.  As  z/zi  increases  from  .01  to  1.0,  the  position  of  the 
spectral  peak  for  w  shifts  to  lower  frequencies  (larger  scales).  For  u  and  v,  there  is 
virtually  no  shift  in  the  peak  for  different  z/zi.  Kaimal  et.  al.  (1976)  also  investigate  km  as 
a  function  of  z/zi  (Figure  2.5).  We  see  from  the  figure  that  for  w,  a  near-linear 
relationship  exists  between  the  dimensionless  peak  wavelength  and  z/zi,  in  the  surface 
layer  (up  to  0.1  Zi)  where 

4/z,=5.9z/z,.  (2.6) 

Above  the  surface  layer,  km  is  seen  to  increase  more  gradually  up  to  about  0.5zi, 
and  approaches  a  constant  value  of  ~L5zi  in  the  range  0.5zi  <  z  <  Zi  .The  horizontal  wind 
components  u  and  v  show  virtually  no  variation  in  km  with  height.  For  temperature,  the 
km  increases  linearly  with  height  in  the  surface  layer,  while  for  z>0.1zi,  km  increases 
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Figure  2,5.  Dimensionless  peak  wavelength  (km)  for  the  veloeity  eomponents  and 
temperature  plotted  as  a  funetion  of  z/zi .  From  Kaimal  et.  al  (1976). 


slightly  with  height  up  to  0.5zi,  then  decreases  slightly  above  0.7zi.  In  general,  the 
results  show  that  in  the  mixed  layer,  the  maximum  wavelength  for  ah  components  is 
roughly  equal  to  1.5zi,  which  corresponds  to  the  length  scale  of  large  thermals  that 
dominate  the  convective  BL  circulation  over  land. 

While  the  previously  mentioned  studies  involved  the  analysis  of  spectra  using 
observations  over  land,  Grossman  (1982)  performed  spectral  analysis  using  aircraft 
observations  of  vertical  velocity  over  the  ocean,  within  the  trade  wind  region  near 
Barbados.  His  goal  was  to  test  the  hypothesis  that  a  hierarchy  of  physical  mechanisms 
(the  major  division  of  this  hierarchy  being  roll  vortex  motion  and  cellular  convection), 
responding  to  wind  shear  and  buoyancy  forcing,  is  active  in  the  generation  of  turbulent 
eddies  within  the  oceanic  sub-cloud  layer.  Using  cases  from  near-  neutral  and  near  free- 
convective  conditions  for  both  along  and  cross  wind  runs,  at  levels  of  18,  46,  and  152  m 
within  fair  weather  BLs  that  averaged  around  600-900  m  in  depth,  Grossman  (1982) 
compares  the  characteristics  of  the  vertical  velocity  spectra,  which,  regardless  of  stability 
or  wind  direction  orientation,  showed  a  shift  toward  larger  values  of  ^  with  increasing 
height.  The  over-ocean  spectra  at  ah  levels  also  reveal  a  double  peak  in  the  vertical 
velocity  spectra,  at  the  1-2  km  scale  and  at  the  lOO’s  of  m  scale.  The  remaining 
discussion  here  will  focus  on  the  cross-wind  vertical  velocity  spectra  at  different  levels 
for  the  case  in  near  free-convective  conditions  (shown  in  Figure  2.6),  which  is  most 
applicable  to  this  study.  The  km  for  the  1-2  km  peak  is  invariant  with  height,  while  the 
lOO’s  of  m  peak  shifts  to  a  longer  km  with  height.  The  variance  for  both  the  1-2  km  peak 
and  the  lOO’s  of  m  peak  at  the  km  increases  with  height.  Based  on  the  relatively  constant 
km  of  the  spectra  peak  at  1-2  km,  he  concludes  that  the  limiting  size  of  eddies  over  the 
ocean  is  roughly  2zi  (1.2  to  1.8  km  for  the  cases  studied),  which  appears  to  apply 
regardless  of  stability  or  orientation  to  the  horizontal  wind  direction.  These  findings  are 
in  close  agreement  with  over-land  studies  of  vertical  velocity  spectra  in  the  surface  layer 
and  ML  discussed  previously.  Grossman  (1982)  states  that  the  spectral  results  presented 
in  his  study  give  insight  to  the  nature  of  physical  mechanisms  occurring  within  these 
BLs.  For  instance,  the  spectral  signature  at  18  m  shows  eddies  on  the  scales  of  lOO’s  of 
m,  likely  generated  by  convective  cells  of  similar  scales,  are  co-existing  with 
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Figure  2.6.  Vertical  velocity  spectra  over  the  ocean  for  a  convective  BL.  Each 
line  represents  the  spectra  from  a  different  vertical  level  within  the  BL,  as 
indicated  (height  in  m).  From  Grossman  (1982). 
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eddies  in  the  km  seales,  likely  generated  by  larger  seale  linear  struetures  sueh  as  roll 
vortiees. 

C.  REPRESENTATION  OF  TURBULENCE  IN  NUMERICAL  MODELS 

1.  Direct  Numerical  Simulations 

In  an  idealized  sense,  a  direet  simulation  of  BL  turbulenee  would  mean  that  the 
turbulenee  is  fully  resolved.  No  turbulenee  elosure  would  be  used,  and  the  goal  is  to 
solve  the  full  equations  of  motions  with  explieit  moleeular  diffusion  terms.  For  this  type 
of  model  to  sueeeed,  the  grid  resolution  of  the  model  must  be  finer  than  the  smallest  seale 
of  turbulent  eddies  in  the  solution.  The  direet  applieation  of  the  model  equations  is 
limited  to  spatial  seales  on  the  order  of  about  1  em  and  time  seales  of  less  than  1  seeond 
(Pielke  2002),  and  to  represent  the  atmosphere  aeeurately,  the  equations  must  be 
evaluated  over  those  spaee  and  time  intervals.  Thus,  a  direet  simulation,  resolving  down 
the  finest  seales  of  turbulent  dissipation,  of  mesoseale  eireulations  (whieh  seale  on  the 
order  of  10  to  100  km  in  the  horizontal,  and  up  to  10  km  in  the  vertieal)  would  require 
solutions  at  lO'^  to  10^^  grid  loeations.  This  far  exeeeds  the  eapaeity  of  any  existing  or 
foreseeable  eomputer  system.  Thus,  there  will  likely  always  be  a  requirement  to  aeeount 
for  unresolved  motions  to  some  degree  in  mesoseale  NWP. 

2,  Large  Eddy  Simulations 

The  Large  Eddy  Simulation  (LES)  is  in  a  middle  ground  between  direet  and 
ensemble-average  simulations.  The  horizontal  grid  resolution  is  high,  typieally  less  than 
100  m  (within  the  inertial  sub-range  of  turbulenee),  with  simulation  time  steps  of  less 
than  10  seeonds.  The  numerieal  simulation  explieitly  deseribes  the  large-seale  turbulenee 
and  assoeiated  fluxes  of  heat  and  momentum,  and  a  turbulenee  elosure  seheme  represents 
the  unresolved  small  seale  motions,  mainly  aeeounting  for  the  dissipation  of  the  resolved 
turbulenee  energy,  at  seales  below  the  inertial  sub-range.  In  general,  EES  models 
perform  well  by  explieitly  resolving  the  large  eddy  energy  well.  The  SGS  energy,  due  to 
the  dissipation  eddies,  is  dependent  on  the  parameterization  seheme,  but  most  studies 
eonelude  that  EES  model  results  are  insensitive  to  the  small  seale  turbulenee  elosure, 
exeept  near  boundaries  and  regions  of  statie  stability  (Mason  1994).  This  insensitivity 
stems  from  the  faet  that  small  eddies  eontribute  mueh  less  heat  and  momentum  transport 
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than  the  large  eddies.  Beeause  these  types  of  simulations  explieitly  resolve  large-seale 
turbulent  eddies,  LES  modeling  has  been  very  sueeessful  in  researeh  as  a  means  of  both 
predieting  the  flow  properties  of  speeifie  turbulent  flows  and  providing  flow  details 
whieh  are  often  used  by  researehers  as  a  substitute  for  real  observations  for  testing  and 
refining  turbulenee  elosure  models.  However,  the  high  vertieal  and  horizontal  resolution 
required  for  LES  modeling  is  eomputationally  restrietive,  limiting  its  use  operationally. 

3,  Ensemble  Turbulence  Closure 

In  contrast  to  direct  simulations,  models  with  ensemble  turbulence  closure  only 
numerically  describe  the  mean  flow,  where  the  equations  of  motion  for  resolved  flow  are 
averaged  continuously  in  space  over  a  volume  AxAyAz.  Therefore,  a  calculated  variable 
centered  at  some  grid  point  represents  the  average  over  the  surrounding  grid  volume. 
The  ensemble  properties  of  all  fluctuations  in  the  flow  (i.e.  turbulent  eddies)  are  assumed 
to  take  place  on  scales  smaller  than  the  resolvable  scale  of  the  model  grid  resolution, 
called  the  sub-grid  scale  (SGS).  These  turbulent  eddies  are  not  modeled  directly,  but 
instead  must  be  represented  through  parameterizations  based  on  the  resolvable  mean 
quantities.  The  accurate  parameterization  of  surface  and  BE  turbulent  fluxes  in  numerical 
models  has  been  a  particularly  challenging  task,  and  the  schemes  used  are  far  from 
perfect,  but  ensemble-average  simulations  are  less  limited  by  computer  resources  than 
direct  or  LES  simulations  and  are  often  the  only  option  in  modeling  some  weather 
systems. 

When  operating  at  high  resolution,  most  mesoscale  models,  including 
COAMPS™,  use  grid  resolutions  that  fall  within  the  scale  of  turbulence.  In  this  case, 
some  of  the  larger-scale  turbulent  eddies  are  explicitly  resolved,  although  the  BE 
turbulent  mixing  processes  are  generally  represented  through  ensemble  parameterization 
schemes.  Some  of  the  different  BE  turbulence  parameterization  schemes  used  in 
mesoscale  NWP  models  were  briefly  discussed  in  chapter  I.  The  BE  turbulence 
parameterization  scheme  specifically  used  in  COAMPS™,  which  will  be  used  in  our 
study,  will  be  described  in  the  next  section. 
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D,  COAMPS™  ENSEMBLE  BL  PARAMETERIZATION  SCHEME 


In  this  section  we  discuss  the  parameterizations  used  to  represent  the  planetary 
boundary  layer  in  the  COAMPS^'^  mesoscale  model.  The  grid-volume  averaging  of 
simplified  prognostic  equations  for  the  mean  variables  within  mesoscale  models  results  in 
averaged  vertical  sub-grid  scale  fluxes  that  must  be  parameterized  since  they  cannot  be 
resolved  by  the  model.  Vertical  sub-grid  scale  fluxes  are  reasonably  well-understood  and 
can  be  parameterized  accurately  in  terms  of  dependent  variables,  i.e.,  the  fluxes  are 
specified  using  experimental  data  and  fundamental  concepts.  In  developing  sub-grid 
scale  averaged  quantities,  the  preferred  representation  is  an  ensemble  average  over  the 
grid  volume,  rather  than  a  simple  grid-volume  average  (Pielke  2002).  The  ensemble 
average  represents  the  most  likely  value  of  the  sub-grid  scale  quantity,  whereas  the  grid- 
volume  average  represents  just  one  realization.  Wyngaard  (1982,  1983)  and  Cotton  and 
Anthes  (1989)  discuss  ensemble  averaging  in  more  depth. 

As  mentioned  in  chapter  I,  COAMPS™  uses  a  one-and-a-half  (1.5)  order  closure 
ensemble  parameterization  scheme  based  on  the  Mellor-Yamada  level  2.5  model  (Mellor 
and  Yamada  1982,  hereafter  MY82)  for  BL  turbulence  parameterization.  The  1.5  order 
closure  retains  the  equations  for  the  vertical  variances  in  the  mean  equations,  with  a 
prognostic  equation  for  TKE  used  in  place  of  the  velocity  variance  equations.  In 
addition,  eddy  coefficients  are  based  on  K-theory  for  the  vertical  flux  terms,  which 
results  in  the  algebraic  simplification  of  the  variance  equations  in  terms  of  empirical 
constants  and  the  parameterized  fluxes.  The  resulting  COAMPS^'^  prognostic  TKE 
equation  is 
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where  e  is  the  TKE  and  1  is  the  vertical  mixing  length  scale  (discussed  in  the  next 
section).  The  second  term  on  the  left  side  represents  the  turbulent  diffusion  of  TKE 
(Mellor  1973). 
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is  the  shear  produetion  term, 


P,=/3gw'e\  (2.9) 

is  the  buoyaney  produetion  term,  with  g  the  aeeeleration  due  to  gravity  and  P  a  eonstant, 
and  0v  the  virtual  potential  temperature.  The  dissipation  rate,  s,  ean  be  parameterized  as 


(2.10) 


where  Ai  is  the  dissipation  length  seale.  The  vertieal  fluxes  are  parameterized  in  the 
prognostie  TKE  (and  in  the  mean  equations)  using  K-theory  as 


-Wu'  =  K—  (2.11) 
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where 
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and  Sm,  Sh,  and  Sg  are  funetions  of  the  gradient  Riehardson  number,  Ri,  where 


R  =A - (2.15). 


dz  J  [  dz 


Inserting  all  of  the  parameterizations  into  the  TKE  equation  gives  the  parameterized  TKE 
equation  used  in  COAMPS™: 


Dfe^]  d\.^d(e^V\  TSeY  ]  e' 
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E, 


FORMULATIONS  OF  VERTICAL  MIXING  LENGTH  SCALES  WITHIN 
BL  TURBULENCE  PARAMETERIZATION  SCHEMES 


Most  turbulence  parameterization  schemes  use,  as  a  closure  assumption,  a  vertical 
mixing  length  scale  (hereafter  referred  to  as  mixing  length),  which  represents  some 
measure  of  the  intensity  of  turbulence  mixing.  Stull  (1988)  details  the  development  and 
use  of  mixing  length  theory,  from  which  the  mixing  length  formulation  is  based  on. 
Mixing  length  formulations  are  often  developed  by  trial  and  error,  in  an  attempt  to  make 
the  simulated  flow  field  match  observed  field  or  laboratory  cases.  Turbulence  models 
have  been  shown  to  be  sensitive  to  the  choice  of  mixing  length  formulation  (Therry  and 
Lacarrere  1983),  and  thus  several  formulations  for  mixing  length  exist  in  the  literature, 
with  some  being  very  specific  to  different  BL  structures  and  characteristics.  To  illustrate 
this,  some  brief  examples  of  different  formulations  are  provided  below.  The  vertical 
profile  of  mixing  length  within  the  BL,  as  modeled  in  COAMPS™,  will  be  graphically 
presented  and  its  formulation  discussed  in  further  detail  in  chapter  IV. 

One  of  the  simplest  formulations  of  mixing  length,  normally  only  appropriate  for 
use  in  LES  modeling,  is  the  approach  first  proposed  by  Smagorinsky  (1963).  The 
Smagorinsky  model  (followed  by  Lilly  1967,  Deardorff  1972,  Klemp  and  Wilhelmson 
1978)  calculates  a  basic  mixing  length, 

/  =  (AvAyAz)‘'\  (2.17) 


where  the  mixing  length  is  dependent  only  on  the  model  horizontal  and  vertical  grid 
resolution  and  which  represents  the  distance  over  which  ensemble  averaging  occurs. 
Some  formulations  have  been  developed  for  very  specific  BL  conditions.  Deardorff 
(1980)  suggested  a  mixing  length  formulation  to  account  for  small  mixing  length  in 
stably  stratified  regions  of  the  form 
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where  e  is  the  TKE.  This  formulation  is  also  used  within  COAMPS^'^  for  the  same 

purpose.  A  variety  of  mixing  length  formulations  have  been  used  successfully  to  model 

the  stable  nocturnal  boundary  layer,  because  of  the  absence  of  large  eddies  (Estournel  and 
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Guedalia  1987;  Lacser  and  Arya  1986;  Yu  1976;  and  Delage  1974).  For  example,  the 
formulation  from  Delage  (1974), 
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where  G  is  the  geostrophie  wind  speed,  f  the  Coriolis  parameter,  p  a  eonstant,  and  L  the 
Obukov  length  seale,  is  an  asymptotie  formulation  that  allows  for  different  atmospherie 
dependeneies  to  be  faetored  into  the  ealeulation  (i.e.,  height,  stability,  wind,  turbulence). 
Equation  2.19  is  based  on  the  formula  developed  by  Blackadar  (1962), 
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where  k  is  the  Von  Karman  constant  and  z  the  height,  and  lo  is  some  maximum  value  of 
length  scale  prescribed  for  a  mixed  layer.  Equation  2.20  is  used  for  the  mixing  length 
formulation  is  many  mesoscale  and  EES  model  parameterization  schemes.  The 
formulation  asymptotically  approaches  two  limits,  where  l~kz  as  z  approaches  zero,  and 
l~lo  as  z  approaches  infinity.  Equation  2.20  is  widely  used  with  lo  constant  or  derived 
using  formulations  of  varying  complexity.  The  MY82  formulation  of  vertical  mixing 
length,  used  within  the  MY82  turbulence  parameterization  scheme  used  in  COAMPS^'^, 
follows  Blackadar  (1962),  with 


(2.21) 


which  represents  the  ratio  of  first  to  zero  order  moments  of  vertical  distribution  of  TKE. 
Physically,  this  means  that  the  length  scale  in  the  upper  portion  of  the  boundary  layer, 
which  is  related  to  the  scale  of  turbulent  eddies,  is  sensitive  to  the  vertically  integrated 
TKE.  The  added  dependence  on  TKE  creates  a  prognostic  equation  for  mixing  length  in 
connection  with  the  TKE  equation,  resulting  in  eddy  coefficients  that  vary  in  time  and 
space  as  the  BE  structure  changes. 

Using  higher-order  closure  techniques  allows  the  use  of  even  more  complex 
formulations  for  mixing  length.  Therry  and  Eacarrere  (1983)  attempt  to  improve  the 
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mixing  length  formulation  for  turbulence  modeling  of  the  convective  BL,  based  on 
results  from  a  third-order  closure  scheme  and  experimental  studies.  They  derive  a  rather 
detailed  formulation  of  both  a  dissipation  and  turbulent  mixing  length  scale,  dependent 
on  TKE  flux  and  which  also  account  for  various  thermal  conditions  and  for  the  larger 
differences  between  vertical  and  horizontal  eddy  motions  as  the  effects  of  buoyancy 
become  larger.  Bougeault  and  Andre  (1986)  develop  a  mixing  length  formulation  for  use 
in  a  third-order  closure  turbulence  model,  specific  for  the  case  of  a  stratocumulus-topped 
BL  where  a  statically  unstable  region  exists  at  the  cloud  top.  Their  formulation  is  based 
on  a  parcel  displacement  scheme,  where  the  upward  and  downward  paths  of  an  imaginary 
parcel  are  calculated,  based  on  the  TKE  at  a  given  height  in  the  BL  and  the  buoyancy 
force  the  parcel  experiences  at  that  height.  The  mixing  length  is  calculated  as  an  average 
between  the  upward  and  downward  paths,  weighted  towards  the  smaller  path  length. 

In  summary,  it  is  evident  by  the  large  number  of  mixing  length  formulations 
proposed  in  the  literature,  that  the  choice  of  appropriate  mixing  length  for  turbulence 
closure  based  on  gradient  transfer  or  K-theory  is  not  always  obvious.  The  use  of  higher- 
order  closure  techniques  come  with  the  expense  of  added  computational  complexity,  and 
so  schemes  such  as  the  MY82  scheme  (used  within  many  of  today’s  mesoscale  models), 
which  retains  the  simplicity  of  K-theory  parameterization,  but  improves  the  model  by 
making  the  mixing  length  a  prognostic  variable  (based  on  TKE),  are  more  desirable.  The 
remaining  chapters  of  this  study  focus  on  the  adequacy  of  the  MY82  mixing  length 
formulation  used  within  the  COAMPS™  BL  parameterization  scheme  in  representing  the 
turbulent  structure  of  the  convective  BL. 
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III.  DETERMINING  THE  SCALE  OF  TURBULENCE  FROM 

OBSERVATIONS 


A,  OBJECTIVE 

To  better  understand  how  the  resolution  of  mesoseale  models  relates  to  the  seales 
of  turbulenee  within  the  BL,  we  performed  speetral  analysis  of  turbulenee  measurements 
from  near  the  surfaee  and  within  the  conveetive  mixed  layer  (ML)  using  high-rate  aircraft 
measurements  from  three  different  experiments  over  the  ocean.  The  objective  is  to 
generalize  the  different  scales  of  turbulence  and  to  directly  evaluate  the  SGS 
representation  of  the  turbulent  fluxes  and  variances.  The  spectral  analysis  performed 
here  will  result  in  turbulent  variances  as  a  function  of  the  length  scale  of  spectral 
integration  for  different  cases  of  convective  BL’s  over  the  ocean. 

Chapter  II  summarized  previous  studies  on  scales  of  BL  turbulence,  the  results  of 
which  will  be  utilized  here  as  a  basis  for  generalizing  the  observed  turbulent  scales.  Our 
focus  here  will  be  on  cases  with  strong  convective  BL’s,  since  it  was  shown  in  Chapter  II 
that  the  convective  BL  contains  the  largest  eddies  which  would  most  likely  approach  the 
scale  of  a  high  resolution  model. 

In  this  chapter,  we  analyze  turbulent  energy  spectra  and  calculate  length  scales, 
using  aircraft  observations  from  the  three  field  experiments  described  below,  in  an 
attempt  to  better  understand  the  contribution  to  turbulent  statistics  from  different  scales  of 
turbulence.  This  analysis  will  provide  the  justification  for  modifications  to  the  COAMPS 
model  to  more  correctly  parameterize  SGS  turbulence  in  high  resolution  simulations 
where  turbulence  is  partly  resolved. 


B,  ANALYSIS  OF  TURBULENCE  SPECTRA  FROM  OBSERVATIONS 
1,  Aircraft  Observations 

a.  Japan/East  Sea  (JES)  Experiment 

The  JES  experiment  took  place  in  January  and  February  2000.  The  main 
objectives  of  the  atmospheric  portion  of  the  JES  experiment  were  to  characterize  the 
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lower  atmosphere  over  the  JES  and  examine  the  role  played  by  the  marine  boundary  layer 
in  determining  air-sea  heat  fluxes  in  the  JES.  This  study  utilizes  the  atmospherie 
measurements  eolleeted  by  the  Twin  Otter  researeh  aircraft  operated  by  the  Center  for 
Interdisciplinary  Remotely  Piloted  Aircraft  Studies  (CIRPAS).  Measurements  utilized,  in 
this  chapter  and  in  later  chapters,  include  SST,  flight  level  wind  (u,  v,  and  w), 
temperature,  and  water  vapor,  as  well  as  derived  flight  level  potential  temperature  and 
moisture  fluxes  and  wind  stress.  The  Twin  Otter  aircraft  speed  was  70  ms"'  and 
measurement  flight  patterns  were  made  up  of  low-level  transects  and  soundings  to  map 
the  surface  flux  and  mean  horizontal  and  vertical  structure  of  the  boundary  layer.  Data 
was  collected  at  a  sampling  rate  of  10  Hz.  A  more  detailed  description  of  the  JES 
experiment  measurements  is  given  in  Kheilf  et  al.  (2003).  While  the  JES  experiment 
covered  a  wide  range  of  wintertime  synoptic  conditions,  this  study  will  focus  on  the  high 
wind  cases  of  frontal  and  post-frontal  large  scale  forcing  in  cold  air  outbreak  conditions. 

The  JES  experiment  Twin  Otter  measurements  from  level  flight  legs  flown 
near  the  surface  are  used  in  this  chapter  for  the  spectral  analysis  of  BE  turbulent  structure. 
In  later  chapters,  the  level  leg  observations  as  well  as  slant  ascent/decent  soundings  from 
JES  are  used  for  inter-comparison  with  high-resolution  COAMPS  output  fields.  A  brief 
summary  of  the  Twin  Otter  flight  legs  used  in  this  study  is  provided  in  Table  3.1. 

b.  Tropical  Ocean  Global  Atmosphere  Coupled  Ocean-Atmosphere 
Response  Experiment  (TOGA  COARE) 

The  TOGA  COARE  intensive  operating  period  (lOP)  took  place  between 
1  November  1992  and  28  Eebruary  1993.  TOGA  COARE  was  a  large  international  field 
experiment  conducted  to  study  the  atmospheric  and  oceanic  exchange  processes  over  the 
region  of  the  western  Pacific  known  as  the  "warm  pool".  This  region,  which  is  near  the 
equator  between  about  140°  E  and  170°  E,  is  a  region  of  very  warm  ocean  temperatures 
(-'28°C)  with  active  convection  and  significant  precipitation  linked  to  the  El  Nino  climate 
variation.  TOGA  COARE’s  primary  objective  was  to  achieve  significantly  more 
accurate  and  complete  descriptions  of  the  surface  meteorology  and  air-sea  fluxes,  with 
the  ultimate  goal  of  producing  a  dataset  that  will  be  used  to  improve  air-sea  interaction 
and  BE  parameterizations  in  numerical  models  and  to  validate  coupled  models.  Weller 
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Date 

Flight 

Leg 

Measurement  Start/ 

End  Time  (UTC) 

Leg  Position  (lat  (degrees  N) 

ion  (degrees  E))  /  Length 

(km) 

Mean 

Attitude  (m) 

02/01/00 

ROl 

0039-0110 

41.0 

133.2-  132.0/90 

41 

02/01/00 

R03 

0124-0144 

41.3 

132.0-  133.0/70 

40 

02/01/00 

R04 

0200-0217 

41.7 

132.7-  132.1  /45 

43 

02/01/00 

R08 

0326-0351 

41.0 

132.1-  133.2  /94 

39 

02/17/00 

R02 

0052-0116 

41.0 

133.0-  132.1  /71 

46 

02/17/00 

R04 

0129-0146 

41.3 

132.1-  132.9/61 

42 

02/17/00 

R06 

0200-0222 

41.7 

132.9-  132.1  /61 
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Date 

Sounding 

Start/ 

End  Time  (UTC) 

Mean  Position  (lat 

(degrees  N)  Ion 

(degrees  E)) 

Sounding  Depth 

Coverage 

(m) 

02/01/00 

S2 

0144-0150 

41.4 

133.0 

36-803 

02/01/00 

S3 

0150-0157 

41.5 

132.8 

803  -  27 

02/17/00 

SI 

0147-0153 

41.4 

132.9 

40  -813 

02/17/00 

S2 

0155-0159 

41.6 

132.9 

808  -40 

02/17/00 

S7 

0409-0414 

41.1 

132.9 

41  -878 

02/17/00 

S8 

0414-0420 

41.2 

132.8 

917  -60 

02/17/00 

S9 

0420-0426 

41.3 

132.7 

58  -782 

02/17/00 

SIO 

0426-0432 

41.4 

132.6 

784  -64 

02/17/00 

Sll 

0432-0436 

41.5 

132.5 

63  -806 

02/17/00 

S12 

0436-0441 

41.6 

132.4 

805  -90 

02/17/00 

S13 

0441-0447 

41.7 

132.3 

88  -  775 

02/17/00 

S14 

0447-0453 

41.8 

132.2 

779  -35 

Table  3.1.  JES  low  level  flight  legs  and  sounding  legs  analyzed  in  this  study. 
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and  Anderson  (1996),  describe  the  wind  conditions  during  the  lOP  as  dominated  by 
periods  of  low  wind  (mean  wind  speed  ~2  ms'^),  with  some  short-lived  events  and 
westerly  wind  burst  events  where  the  wind  speed  reached  ~9-12  ms"'  for  a  short  duration. 
These  wind  bursts  excite  the  local  air-sea  heat  and  moisture  fluxes,  enhancing  the 
convection  and  BL  turbulence  even  more. 

Measurements  from  TOGA  CO  ARE  utilized  in  this  study  include 
horizontal  and  vertical  wind  velocity,  temperature,  and  moisture  fields  collected  by  the 
National  Center  for  Atmospheric  Research  (NCAR)  Electra  aircraft,  which  flew  32 
flights  to  study  boundary  layer  fluxes  and  convective  systems.  The  Electra  aircraft  made 
turbulence,  radiation,  and  cloud  microphysics  measurements  in  the  BE.  The 
measurements  used  in  this  study  consisted  of  horizontal  legs  at  different  altitudes  in  and 
above  the  marine  BE,  and  vertical  soundings  throughout  the  BE.  The  aircraft  recorded 
measurements  at  a  constant  speed  of  100  m  s' ^  Horizontal  flight  legs  consisted  of  20- 
minute  E-shaped  legs  at  various  vertical  levels  within  and  above  the  BE.  Observations 
captured  cloud  structure  and  measured  the  turbulent  flux  profdes  at  various  levels  in  the 
BE.  Weller  and  Anderson  (1996)  describes  the  Electra  atmospheric  data  collection 
instrumentation  used  in  TOGA  COARE.  The  data  was  sampled  at  a  50  Hz  sampling  rate 
and  recorded  at  20  Hz  for  data  processing.  These  measurements  were  obtained  from  the 
data  archives  at  the  National  Center  for  Atmospheric  Research  (NCAR).  Descriptions  of 
the  synoptic  and  BE  meteorological  conditions  of  the  TOGA  COARE  flights  analyzed 
are  based  on  information  from  the  TOGA  COARE  lOP  Operations  Summary  (1993). 
These  descriptions  will  be  used  to  explain  the  characteristics  of  the  BE  turbulence 
structure  along  these  flight  tracks.  Table  3.2  lists  the  flight  date,  time,  average  BE  depth 
(as  determined  from  sounding  measurements),  and  average  level  leg  length  and  wind 
speed  for  the  boundary  layer  flight  measurement  legs  that  were  used  in  this  study.  In 
cases  where  convective  clouds  were  present,  the  lifting  condensation  level  (ECE)  was 
calculated  from  sounding  measurements  and  taken  as  the  BE  height  (Stull  1988). 
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Fit 

# 

Date 

Measurement 

Start/ 

End  Time  (UTC) 

Avg  BL  Depth  or 

LCL  (m) 

Average  BL 

level  leg  length 

(km) 

Avg.  BL  wind 

speed 

(m/s) 

10 

12/9/92 

2230-0300 

300 

119 

2.5 

12 

12/10/92 

1830-2200 

450 

52 

3.0 

19 

1/13-14/93 

2230-0200 

450 

59 

2.0 

20 

1/14-15/93 

2230-0200 

500 

60 

1.9 

22 

1/17-18/93 

0030-0330 

400 

119 

4.9 

23 

1/19/93 

0045-0330 

910 

122 

7.4 

24 

1/26-27/93 

2230-0050 

250 

49 

4.8 

26 

1/28-29/93 

2245-0200 

570 

61 

6.5 

Table  3.2.  TOGA  COARE  level  flight  legs  analyzed  in  this  study. 
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c.  FIRE  Experiment 


During  June  and  July  1987,  the  First  International  Satellite  Cloud  and  Climate  Projeet 
(ISCCP)  Regional  Experiment  (FIRE)  Marine  Stratoeumulus  Intensive  Field 
Observations  (IFO)  was  eondueted  off  the  eoast  of  southern  California.  FIRE  obtained  a 
eomprehensive  set  of  boundary  layer  and  above  boundary  layer  measurements  in  12 
eoordinated  multi-missions.  Albreeht  et  al.  (1988)  and  Cox  et  al.  (1987)  give  a  detailed 
deseription  of  projeet  objeetives.  In  addition,  Kloesel  et  al.  (1988)  gives  a  summary  of 
experiment  operations  and  general  meteorologieal  eonditions  during  the  period.  In  this 
ehapter,  we  utilize  horizontal  turbulenee  leg  data  obtained  by  the  NCAR  Eleetra  aireraft 
(as  deseribed  above  for  the  TOGA  COARE  experiment)  during  FIRE  on  July  3,  1987 
(flight  3),  at  various  levels  in  the  boundary  layer. 

2,  Statistical  Calculations  for  Spectral  Estimates 

In  this  study,  we  used  MATLAB,  the  data  aequisition,  analysis,  and  visualization 
software  applieation  eommonly  used  in  the  seientifie  eommunity,  as  an  integral  tool  for 
the  detailed  analysis  of  both  model  output  and  the  real  data.  For  the  graphieal 
representation  of  the  speetra  analysis,  the  MATLAB  funetion  PSD  was  used  to  estimate 
the  power  speetral  density  of  a  diserete  time  series  using  Weleh’s  averaging  method.  The 
time  series  is  linearly  detrended,  whieh  reduees  the  eontamination  of  the  speetra  by 
wavelengths  longer  than  the  data  series.  The  power  speetral  density  as  a  funetion  of 
wave  number  is  then  estimated  to  a  95%  eonfidenee  interval.  The  speetra  are  ealeulated 
between  the  smallest  resolvable  seale  (or  largest  resolvable  (Nyquist)  frequeney),  defined 
by  the  data  frequeney  (for  observations),  and  the  largest  resolvable  seale,  defined  by  the 
length  of  the  observation  leg. 

Using  the  energy  density  speetrum,  we  also  ealeulate  the  Ogives,  defined  as  the 
eumulative  integral  of  the  speetrum  from  the  smallest  wavelength,  (where  the 
eontribution  to  the  energy  is  generally  small)  to  some  speeified  wavelength  of  speetral 
integration  ()io)  (Rogers  et  al.  1995).  Theoretieally,  Ogives  ean  be  expressed  as 
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where  the  denominator  is  the  total  integrated  energy  over  all  wavelengths,  or  variance, 
which  normalizes  the  calculation.  In  discrete  form,  we  set  the  lower  limit  in  integrations 
to  the  wavelength  corresponding  to  the  Nyquist  frequency  (smallest  resolvable 
wavelength)  for  both  integrals  in  Equation  3.1,  and  the  upper  limit  to  3  km  (generalized 
estimate  of  the  largest  turbulent  eddy  size  for  the  cases  analyzed  here),  in  the  variance 
calculation  in  order  to  only  consider  spectral  energy  associated  with  turbulence.  The 
resulting  plot  of  the  Ogives  versus  ko  allows  us  to  analyze  the  fraction  of  total  turbulent 
energy  from  scales  smaller  than  a  given  scale  limit. 

3,  Spectral  Analysis 

a.  Method  and  Objectives 

Here  we  present  examples  of  the  energy  density  spectrum  of  vertical 
velocity  as  a  function  of  wave  number,  and  the  Ogives  vs.  wavelength  of  spectral 
integration.  In  these  examples,  wavelength  (k)  and  wave  number  (k)  are  related  by 

k  =  27i/k,  (3.2) 

where  k  is  calculated  from  the  sampling  frequency  and  the  aircraft  speed.  For  reference, 
a  wave  number  of  1  km'^  corresponds  to  a  wavelength  of  about  6  km. 

The  spectra  are  presented  in  log-linear  fashion,  so  that  the  area  under  the 
energy  curve  is  proportional  to  the  variance,  and  a  peak  in  the  spectrum  can  be 
interpreted  as  the  wavelength  that  contributes  most  to  the  variance.  As  the  scale 
increases,  the  statistical  significance  and  thus  the  smoothness  of  the  spectral  estimate 
decreases.  However,  the  length  of  all  observation  legs  in  these  examples,  which  varies 
from  70  to  130  km,  is  considered  adequate  for  the  spectral  representation  of  the  entire 
turbulence  scale,  which  is  considered  here  to  be  at  wave  numbers  2  km'^  and  higher 
(turbulent  eddies  3  km  and  smaller).  For  better  presentation  of  the  turbulent  energy 
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spectra,  the  signal  from  scales  larger  than  6  km  was  filtered  out  of  the  observed  vertical 
velocity,  to  remove  signals  from  the  less  statistically  significant  waves. 


The  spectra  analysis  will  focus  primarily  on  two  aspects:  1)  determination 
of  the  peak  wavelength,  which  will  indicate  the  size  of  turbulent  eddies  dominating  the 
fluctuations  of  turbulent  variables  within  the  BL,  and  2)  the  Ogives,  or  fraction  of  total 
turbulent  energy  from  scales  smaller  than  a  given  scale  limit.  Ogives  allows  us  to 
determine  the  fraction  of  total  turbulent  energy  that  could  potentially  be  resolved  by  a 
high-resolution  mesoscale  model,  and  thus  the  appropriateness  of  an  ensemble  turbulence 
parameterization. 

While  TKE  depends  on  the  variance  of  both  horizontal  and  vertical  wind, 
we  choose  to  present  only  the  spectra  of  vertical  velocity.  The  horizontal  wind  spectra 
were  analyzed  as  well,  and  the  shape  and  turbulent  scales  are  quite  similar  to  that  of  the 
vertical  velocity  spectra.  The  main  difference  is  the  significant  mesoscale  contribution  to 
the  spectral  energies  in  the  horizontal  wind  spectra,  which  are  not  present  in  the  vertical 
velocity  spectra.  While  vertical  motion  is  limited  by  the  depth  of  the  BL,  there  are  no 
boundaries  on  the  horizontal  flow  to  restrict  the  scales  of  motion.  These  mesoscale 
contributions  to  the  horizontal  wind  spectra  are  also  noted  in  Nucciarone  and  Young 
(1991)  for  the  FIRE  experiment. 

b.  Surface  Layer  Spectra 

Figures  3.1  through  3.3  show  three  examples  of  turbulent  characteristics  of 
the  surface  layer  of  the  convective  BL  using  near-surface  high-rate  vertical  velocity 
measured  by  aircraft.  The  surface  layer  is  defined  as  a  height  (z)  within  the  BL  that  is  less 
than  or  equal  to  1/10  of  the  BL  height  (zi).  Using  a  normalized  height  scale  z*  equal  to 
z/zi,  the  surface  layer  is  then  defined  for  a  z*  of  0.1  or  less.  Figure  3.1  is 
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JES  Aircraft  Observations  -  Mean  Spectra 


Distance  along  track  (km) 


Wavelength  (km) 


Figure  3,1.  Vertical  Velocity  (w)  along  the  flight  leg  (a),  energy  density  spectrum  of  w 
vs.  wave  number  (b)  and  Ogives  of  w  (c)  vs.  wavelength,  using  near-surface  observations 
from  the  JES  experiment. 
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TOGA  COARE  Aircraft  Observations  -  Leg  FI  23b 


0  10  20  30  40  50  60 


Distance  along  track  (km) 


Figure  3.2,  Same  as  Figure  3.1,  exeept  using  observations  from  a  near-surfaee  flight  leg 
during  flight  23  of  the  TOGA  COARE  experiment. 
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FIRE  Aircraft  Observations  -  Leg  HI  3a 
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Figure  3.3.  Same  as  Figure  3.1,  exeept  using  aireraft  observations  from  a  near-surface 
flight  leg  during  flight  3  of  the  FIRE  experiment. 
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an  example  from  the  JES  experiment,  on  January  31,  2000.  The  surfaee  layer  was 
unstable  (z/L  =  -0.35),  with  an  air-sea  temperature  difference  of  nearly  10  degrees.  The 
horizontal  wind  speed  was,  on  average,  greater  than  12  ms‘\  For  this  flight,  the  BE  depth 
averaged  around  700  m  as  determined  from  sounding  profdes.  The  spectrum  and  Ogives 
in  Figure  3.1  are  averages  of  four  different  low-level  legs  from  this  flight,  flown  at  40  m 
above  the  sea  surface  (z*  =  0.07).  Figure  3.2  shows  an  example  from  flight  23  of  the 
TOGA  CO  ARE  experiment.  For  this  flight,  the  surface  layer  was  also  unstable  (z/F  =  - 
0.86),  with  an  air-sea  temperature  difference  of  3  degrees  and  a  mean  BE  horizontal  wind 
speed  of  about  8  ms'^  The  BE  height  (taken  as  the  lifting  condensation  level)  averaged 
around  400  m.  The  observations  used  are  from  a  110  km  long  flight  leg  flown  at  31  m 
(z*=  0.08)  above  the  sea  surface.  Figure  3.3  shows  a  typical  vertical  velocity  spectra  and 
Ogives  from  flight  3  of  the  FIRE  experiment.  During  this  night  mission,  the  surface  layer 
was  unstable  (z/F  =  -0.38)  with  an  air-sea  temperature  difference  of  0.8  degrees,  and  an 
average  BE  horizontal  wind  that  exceeded  10  ms‘\  The  BE  was  very  well  mixed  and 
turbulent,  with  an  average  depth  of  about  750  m.  The  observations  used  are  from  a  120 
km  long  flight  leg  flown  at  56  m  above  the  sea  surface  (z*=  0.08). 

We  can  see  distinct  features  of  the  near-surface  turbulence  for  the 
observations  in  these  three  cases  of  unstable  surface  layers.  Comparing  the  top  panels  in 
each  of  the  Figures  3.1  through  3.3,  the  JES  example  clearly  has  the  largest  vertical 
velocity  variation.  This  fact  is  also  evident  in  the  larger  magnitude  of  the  spectra 
compared  to  the  TOGA  COARE  and  FIRE  examples.  For  all  three  examples,  the  spectral 
signature  of  turbulent  energy  appears  clearly  at  wave  numbers  (wavelengths)  to  the  right 
of  about  2  km‘'  (3  km)  on  the  x-axis,  and  the  peak  scale  of  turbulence  eddies  (km) 
generally  exists  within  the  wavelength  range  of  200  to  300  m,  representing  the  dominant 
eddy  size  contributing  to  the  total  TKE  at  the  measurement  level  near  the  surface.  The 
inconsistency  in  the  shape  of  the  spectra  at  larger  scales  between  the  three  cases  with 
varying  values  of  z/E  is  consistent  with  the  findings  of  Kaimal  et  al.  (1972)  for  vertical 
velocity  spectra  in  the  unstable  surface  layer.  When  plotted  on  a  log-log  scale  (not  shown 
here),  the  inertial  sub-range  in  the  near-surface  energy  spectra  (as  well  as  the  spectra 
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from  higher  levels  in  the  ML,  shown  later)  follow  the  -5/3  law  of  the  inertial  sub-range, 
and  the  seale  of  the  inertial  sub-range  (whieh  appears  to  begin  around  60  km'^  (100  m 
wavelength)  in  all  three  cases),  does  not  vary  among  the  different  cases  (Kaimal  et  al. 
1972).  In  the  JES  case,  we  observe  the  existence  of  a  double-peak  in  spectral  intensities, 
as  was  observed  in  Grossman  (1982)  (see  chapter  II).  The  peak  near  I  km,  while  not 
dominant  near  the  surface,  signifies  the  existence  of  observed  horizontal  roll  vortices 
(discussed  chapter  IV)  produced  in  a  post-frontal  cold  air  outbreak  over  the  relatively 
warmer  ocean.  The  much  larger  magnitude  spectra  in  the  JES  ME  case,  compared  to  the 
examples  from  TOGA  COARE  and  EIRE,  is  due  to  the  relatively  more  intense  turbulence 
forced  by  large  surface  fluxes  associated  with  the  unstable  surface  layer. 

The  Ogives  plots  in  Eigures  3.1  through  3.3  indicate  that,  near  the  surface, 
50%  of  the  total  turbulent  variance  is  contained  in  eddies  smaller  than  about  200  m  in 
scale.  The  largest  eddies,  greater  than  1  km  in  scale,  make  up  about  10%  of  the 
turbulence  variance  near  the  surface.  This  implies  that  a  mesoscale  model  operating  at  a 
horizontal  resolution  of  500  m  could  theoretically  resolve  a  1  km  (2Ax)  wave,  and  thus 
10%  of  the  turbulent  energy  near  the  surface. 

c.  Mixed  Layer  Spectra 

Vertical  velocity  measurements  from  the  TOGA  COARE  and  EIRE 
experiments  taken  from  within  the  ML  of  the  convective  BL  are  also  analyzed  to 
determine  the  scale  of  turbulence  at  these  levels.  The  vertical  velocity  along  the  leg, 
spectra,  and  Ogives  plots  in  Eigure  3.4  are  from  flight  23  of  the  TOGA  COARE 
experiment.  The  observations  used  are  from  a  flight  leg  flown  at  221  m  above  the  sea 
surface,  in  the  middle  of  the  ME  (z*=  0.55).  Eigure  3.5  illustrates  these  same  parameters 
for  flight  3  of  the  EIRE  experiment,  from  a  leg  flown  at  481  m  above  the  sea  surface,  also 
near  the  middle  of  the  ML  (z*=  0.64). 

Compared  to  the  near  the  surface,  we  observe  distinct  and  consistent 
differences  in  the  ML  turbulence  scales,  fairly  consistent  with  the  findings  of  Grossman 
(1982)  for  over  the  ocean  MLs.  Eor  comparison  purposes,  the  near-surface  and  ML 
spectra  are  plotted  on  the  same  axis  in  Eigures  3.6  (for  EIRE)  and  3.7  (for  TOGA 


49 


TOGA  COARE  Aircraft  Observations  -  Leg  H3  3a 
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Figure  3.4.  Same  as  Figure  3.2  except  using  aircraft  measurements  from  a  flight  leg 
flown  at  221  m  above  the  sea  surface. 
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FIRE  Aircraft  Observations  -  Leg  H2  3b 
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Figure  3.5.  Same  as  Figure  3.3,  exeept  using  aireraft  measurements  from  a  flight  leg 
flown  at  481  m  above  the  sea  surfaee. 
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Figure  3.6.  Vertical  velocity  spectra  from  FIRE  aircraft  observations,  from  figures  3.3 
and  3.5,  plotted  together  for  comparison. 
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COARE).  In  general,  the  magnitude  of  the  spectral  peak  at  the  km  scale  is  much  larger 
than  near  the  surface,  and  the  eddies  at  the  km  scale  are  now  dominant.  The  km  has  not 
changed  significantly  in  magnitude  nor  has  it  shifted  much  compared  to  the  near-surface 
spectra.  The  energy  drop  at  6  km  (i.e.,  the  filter  scale)  is  now  more  abrupt  than  it  was 
near  the  surface,  since  there  is  now  a  significant  amount  of  energy  at  larger  end  of  the 
turbulent  scale  near  the  spectral  gap.  These  differences  between  the  surface  and  ML 
spectra  are  consistent  in  both  the  FIRE  and  TOGA  COARE  spectra.  We  also  observed 
complexity  in  the  large  scales  of  the  individual  spectra,  where  there  are  several  peaks 
with  significant  energy,  and  no  single  dominant  km.  These  variations  are  expected  since 
the  calculations  are  from  only  one  flight  leg. 

For  the  TOGA  COARE  case,  km  has  increased  from  about  400  m  near  the 
surface,  to  between  600  m  and  1 .5  km  in  the  ML,  and  the  spectral  magnitude  near  the  1 
km  scale  has  nearly  doubled  (evident  from  the  larger  variation  of  the  vertical  velocity 
along  the  flight  leg  (top  panel  in  the  figures),  compared  to  near  the  surface.  With  a  BE 
height  of  400  m,  the  energetic  large-scale  eddies  have  a  size  of  about  1.5  Zi,  consistent 
with  the  findings  from  Kaimal  (1976).  More  dramatic  differences  are  seen  between  the 
surface  and  ML  spectra  from  FIRE.  Due  to  the  more  energetic  large  scale  eddies,  the  ML 
spectrum  is  nearly  four  times  larger  in  magnitude  in  the  larger  scales,  near  1  km, 
compared  to  near  the  surface.  The  more  drastic  increase  in  magnitude  between  the  near¬ 
surface  and  ML  spectra  compared  to  that  of  the  TOGA  COARE  example  is  most  likely 
due  to  the  higher  BE  and  more  turbulent  ML  of  the  FIRE  case  (see  Wang  and  Albrecht 
1994).  Unfortunately,  no  ML  flight  legs  were  flown  during  the  JES  experiment  to 
include  in  this  comparison.  The  empirical  relationship  in  Kaimal  (1976)  appears  to  be 
valid  since  it  is  supported  by  the  examples  shown  here  from  FIRE  and  TOGA  COARE. 

Consistent  with  the  shift  in  turbulent  energy  to  larger  scales,  the  Ogives 
calculations  for  the  ML  examples  (lower  panel.  Figures  3.4  and  3.5)  show  a  shift  in  the 
fraction  of  turbulent  energy  at  a  given  scale.  We  see  that  the  ML  Ogives  indicate  only 
about  15  to  20%  of  the  total  turbulence  variance  at  wavelengths  below  600  m,  over  30% 

of  the  total  turbulent  eddy  energy  scales  above  1  km.  A  mesoscale  model  operating  at 
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500  m  resolution  could  theoretically  resolve  four  times  more  turbulent  energy  in  the  ML 
(30%),  compared  to  near  the  surface  (10%),  implying  that  the  resolvability  issue  is  even 
more  important  within  the  ML. 

d.  Conclusions 

The  purpose  of  the  spectral  analysis  was  to  determine  the  scales  at  which 
turbulence  in  the  atmospheric  BL  occurs,  and  to  determine  the  extent  to  which  a  high- 
resolution  mesoscale  model  using  ensemble  parameterization  accounts  for  that 
turbulence.  The  vertical  velocity  spectra  from  over-the-ocean  convective  BLs  near  the 
surface  and  within  the  ML  were  analyzed.  No  overland  cases  were  shown  here,  but  past 
studies  by  Kaimal  (1976)  and  Young  (1987)  show  that  the  vertical  velocity  spectra  over 
land  are  similar  in  terms  of  shape  and  height  dependence  to  the  spectra  over-the-ocean, 
but  have  a  somewhat  larger  peak  wavelength.  While  it  is  difficult  to  make  general 
conclusions  from  a  limited  number  of  cases,  our  results  were  consistent  and  in  agreement 
with  previous  studies.  In  addition,  we  discovered  that  while  the  majority  of  turbulence  in 
the  convective  ML  indeed  occurs  at  scales  near  1  km,  a  significant  fraction  of  the  total 
turbulence  can  extend  to  larger  scales  near  the  spectral  gap,  blurring  the  convenient 
separation  between  the  mesoscale  and  turbulent  scales. 

The  Ogives  plots  suggest  that  nearly  30%  of  the  turbulent  eddy  energy  in 
the  ML  exists  at  scales  above  1  km.  This  suggests  that  as  the  resolution  of  a  mesoscale 
model  approaches  the  sub-kilometer  scales,  the  model  would  partially  resolve  a 
significant  portion  of  the  total  BL  turbulent  energy.  Theoretically,  an  ensemble 
turbulence  parameterization  represents  the  turbulence  from  all  scales,  so  if  indeed  the 
larger  scale  turbulence  is  theoretically  resolved  by  a  high  resolution  mesoscale  model,  the 
SGS  parameterization  should  only  account  for  part  of  the  total  turbulence  variance.  The 
spectral  analysis  performed  here  thus  provides  an  indication  of  the  amount  of  turbulent 
energy  we  can  expect  a  high-resolution  mesoscale  model  BL  ensemble  turbulence 
parameterization  to  account  for  at  different  levels  within  the  convective  BL.  The  results 
suggest  that  it  is  no  longer  appropriate  to  use  ensemble  turbulence  closure  in  a  mesoscale 
model  when  the  grid  resolution  goes  below  1  km. 
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C.  TURBULENCE  INTEGRAL  LENGTH  SCALES  FROM  OBSERVATIONS 
1.  Importance  of  Integral  Length  Scale  Analysis 

To  further  understand  the  seale  of  BL  turbulenee,  the  same  aireraft  observations 
used  in  the  spectral  analysis  are  used  here  to  obtain  the  integral  length  scale  (defined  in 
chapter  II).  The  purpose  here  is  to  quantify  the  length  scale  and  its  variation  with  height 
in  the  BL.  This  analysis  will  allow  us  to  determine  the  appropriateness  of  the  mixing 
length  used  in  ensemble  turbulence  parameterizations  when  horizontal  resolution  is 
smaller  than  the  scale  of  the  largest  turbulent  eddies  in  the  BL,  and  will  be  the  basis  for 
adjustments  made  to  the  ensemble  turbulence  parameterization  within  COAMPS,  via  the 
mixing  length  calculation,  to  account  for  the  partial  resolution  of  the  turbulent  eddies. 

2,  Integral  Length  Scale  Calculation 

a.  Formulation 

For  each  horizontal  flight  leg  used  in  the  subsequent  analysis,  the  integral 
length  scale  of  vertical  velocity  is  calculated  for  FIRE  flight  3  and  the  TOGA  COARE 
flights  using  Equation  2.2.  The  autocorrelation  coefficient  (p(t))  is  calculated  by  finding 
the  covariance  of  the  vertical  velocity  and  the  time  lag  at  each  point  along  the  flight  leg 
(the  co-variance  is  normalized  so  that  the  co-variances  at  zero  lag  are  identically  1.0). 
p(t)  is  the  time  lag  corresponding  to  the  point  along  the  leg  where  the  normalized 
covariance  first  falls  below  e"'.  U  is  the  aircraft  speed  relative  to  the  wind  speed  in  the 
direction  of  the  flight  path. 

b.  Calculation  of  Partial  Integral  Length  Scales 

Using  a  3’^‘^  order  Butterworth  filter  in  MATEAB,  a  high  pass  filter  was 
designed  to  remove  low  frequency  perturbations  above  a  specified  filter  scale  E.  Eigure 
3.8  shows  an  example  of  the  application  of  the  filter  function,  where  the  energy  from 
eddies  larger  than  500  m  in  scale  is  filtered  out.  Note  that  the  filtering  can  be  observed  in 
the  time  series  of  the  filtered  vertical  velocity,  where  the  low  frequency  oscillations  due 
to  the  larger  eddies,  seen  in  the  unfiltered  time  series,  are  not  present. 
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Figure  3.8.  Example  of  time  series  and  spectra  of  vertical  velocity  after  filtering  out 
energy  from  eddies  scaling  larger  than  500  m.  Observations  are  taken  from  a  horizontal 
flight  leg  during  FIRE  Flight  3. 
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The  filter  is  applied  to  the  FIRE  flight  3  and  TOGA  COARE  data  for  the 
calculation  of  p(t).  The  result  is  the  calculation  of  partial  integral  length  scale,  defined  as 
the  integral  scale  of  turbulence  below  a  certain  filter  scale  E.  The  ensemble  integral  scale 
representing  all  scales  of  turbulence  was  calculated  by  filtering  out  the  signal  from  waves 
larger  than  3  km  in  the  vertical  velocity  time  series,  ensuring  that  the  energy  from  all 
turbulent  eddies  are  considered  in  the  calculation.  As  implied  by  our  results  in  the 
spectral  analysis  section,  a  3  km  filter  would  capture  the  energy  from  all  scales  of 
turbulence  for  the  FIRE  and  TOGA  COARE  flight  data,  and  is  likely  appropriate  for 
most  convective  BLs. 

This  filter  procedure  was  used  to  calculate  the  integral  length  scale  for  a 
range  of  E  from  3  km  down  to  200  m.  The  partial  integral  scales  from  along  and  cross 
wind  legs  (flown  in  succession  at  the  same  level)  were  averaged  together,  based  on 
Eenschow  and  Stankov  (1986),  where  it  was  determined  that  there  is  only  a  slight 
difference  between  integral  scales  of  the  wind  components  calculated  from  along  and 
cross  wind  flight  legs. 

3,  Results 

Following  Eenschow  and  Stankov  (1986),  we  plot  the  partial  integral  length  scale, 
normalized  by  the  BE  height,  for  the  different  filter  scales  (Figure  3.9).  The  different 
symbols  in  Figure  3.9  represent  the  partial  integral  length  scale  for  different  L,  as 
calculated  from  the  observations.  Using  the  normalized  height  scale  z*,  we  find  that  for 
both  the  FIRE  and  TOGA  COARE  observations,  the  ensemble  integral  length  scale  (i.e., 
at  the  3  km  filter  scale)  for  vertical  velocity  (Aens)  is  well  represented  at  most  z*  in  the 
BE  by 

A^^^/z.  =0.43zy"  .  (3.2) 

This  relationship  also  followed  for  the  integral  length  scale  for  vertical  velocity  of 
Eenschow  and  Stankov  (1986),  using  aircraft  measurements  taken  in  convective  BEs  over 
both  ocean  and  land  surfaces.  The  relationship  does  not  fit  the  observations  as  well  near 
the  top  of  the  BE,  where  the  integral  scale  depends  more  on  the  stability  of  the  overlying 
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Int.  Scale  vs.  z/z.  for  various  filter  scales  (L)  -  TOGA  COARE  and  FIRE 


Figure  3,9,  Integral  length  scale,  normalized  by  the  BL  height  (zi),  for  different  filter 
scales  (upper  limit  on  the  scale  of  motion)  (L).  Calculations  used  observations  from 
several  horizontal  flight  legs  at  various  levels  in  the  BL  from  FIRE  flight  3  and  TOGA 
COARE.  Different  symbols  represent  the  partial  integral  scale  calculated  along  the  flight 
leg  for  different  E,  as  shown  in  the  legend.  The  line  that  fits  to  the  observations  for  E  =  3 
km  is  calculated  from  Equation  3.2.  The  line  that  best  fits  the  observations  for  L  =  200  m 
is  a  vertical  line.  Eines  that  fit  to  the  observations  for  L  between  200  m  and  3  km  are 
calculated  from  Equation  3.3. 
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layer  and  less  on  the  scale  of  turbulent  BL  eddies  (Carruthers  and  Hunt  1986). 

As  L  decreases,  such  that  only  smaller  and  smaller  scales  of  motion  are 
considered  in  the  calculation,  the  integral  length  scale  expectedly  decreases  at  all  levels  in 
the  BL,  as  shown  in  Figure  3.9.  The  difference  between  integral  scales  for  different  L  is 
largest  in  the  middle  and  upper  BL,  where  the  majority  of  turbulent  eddies  are  large,  and 
smallest  near  the  surface,  where  the  size  of  most  eddies  is  small.  As  L  decreases,  the 
integral  length  scale  also  becomes  less  variant  with  height  in  the  BL,  and  we  see  that  for 
L  =  200  m,  representing  eddies  within  the  inertial  sub-range,  there  is  practically  no 
variation  in  the  integral  scale  with  height. 

Attempts  were  made  to  represent  the  partial  integral  length  scale  using  the 
ensemble  length  scale  (Aens)  and  that  from  the  200  m  cutoff  filter  (A2oom)-  The  latter  is 
considered  to  represent  the  SGS  length  scale  for  a  LES  simulation,  although  the  cutoff 
limit  for  LES  is  generally  smaller  than  200  m.  Equation  3.3  is  selected  among  many 
other  alternatives  for  several  reasons.  Eirst,  this  relationship  should  be  applicable  for  the 
ensemble  parameterizations  currently  used  in  COAMPS.  Second,  it  should  generate  a 
length  scale  that  approaches  what  is  normally  used  in  EES  as  the  cutoff  length  falls  into 
the  range  of  the  inertial  sub-range.  The  above  requirements  ensure  the  asymptotic 
property  of  the  length  scale  representation  at  large  and  small  E.  Third,  it  should  fit  the 
observed  partial  integral  length  scale  shown  in  Eigure  3.9  for  intermediate  E,  where  E 
relates  to  a  horizontal  grid  resolution  with  E  =  2  Ax.  Therefore,  in  a  high-resolution 
mesoscale  model,  the  partial  integral  length  scale  using  Equation  (3.3)  is  essentially  a 
weighted  average  between  the  ensemble  and  the  EES  length  scales,  where  the  weighting 
varies  with  Ax: 
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(3.3) 


The  parameter  R  in  Equation  (3.3)  should  be  obtained  empirically.  Eor  the  EIRE  and 
TOGA  COARE  convective  BE  observations,  R  is  determined  to  be  3.8  km'\  At 

Ax  =  500  m,  is  about  0.85. 
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The  physical  processes  that  may  affect  R  is  worth  discussion.  Since 

provides  a  weighting  between  the  ensemble  and  the  minimum  length  scales, 
conceivably,  it  is  largely  affected  by  the  shape  of  the  energy  density  spectrum.  As  a 
result,  factors  affecting  the  property  of  the  energy  density  spectrum  should  eventually 
affect  the  magnitude  of  R.  The  results  here  were  obtained  from  measurements  within 
convective  BLs  over  the  ocean.  In  general,  R  may  differ  as  the  BL  stability  varies,  which 
is  a  main  factor  affecting  the  boundary  layer  turbulent  spectra.  However,  as  revealed  by 
Kaimal  et  al.  (1976),  the  stability  dependence  of  the  turbulence  spectra  in  the  convective 
boundary  layer  is  much  weaker  compared  to  that  of  the  stable  boundary  layer. 
Consequently,  we  do  not  expect  significant  variation  of  R  for  other  cases. 

Equation  (3.3)  will  not  approach  the  LES  length  scale  presently  used  in  most  of 
the  EES  models  at  Ax  smaller  than  100  m,  the  minimum  filter  scale  used  in  our  length 
scale  calculation.  Ideally,  instead  of  using  a  constant  A2oom  ,  one  should  use  the  EES 
length  scale  which  is  also  a  function  of  the  filter  scale  Ax.  This  is  left  for  future  research. 

This  experiment  illustrates  two  major  points.  Eirst,  we  can  make  the  analogy  that, 
for  a  given  model  resolution  Ax,  the  turbulent  length  scale  (mixing  length)  required  to 
appropriately  represent  SGS  turbulence  is  represented  by  the  observations  and  associated 
curve  for  an  E  of  2Ax  as  shown  in  Eigure  3.9.  It  is  clear  that  the  mixing  length  profile 
that  accounts  for  all  scales  of  turbulent  eddies  in  the  BE  (i.e.,  Aens)  is  inappropriate  for 
the  COAMPS  BE  parameterization  scheme  when  the  resolution  is  less  than  around  1.5 
km.  Thus,  the  parameterization  scheme  overestimates  the  SGS  turbulent  fluxes  and  the 
TKE,  based  on  Equations  2.7  through  2.16  in  chapter  II.  The  formulation  we  derive  from 
observations  is  the  basis  for  an  appropriate  scale-dependent  mixing  length  for  model 
turbulence  parameterization  schemes  when  the  model  resolution  is  smaller  than  the  size 
of  large  turbulent  eddies.  This  formulation  also  provides  a  smooth  transition  from  an 
ensemble  mixing  length,  when  Ax  is  large,  to  a  typical  EES  turbulent  parameterization 
mixing  length,  when  Ax  is  small.  In  essence,  the  appropriate  mixing  length  profile  for  a 
given  Ax  in  Equation  3.3  is  a  weighted  average  between  a  typical  mixing  length  profile 
appropriate  for  an  EES  model  (appropriate  for  eddies  scaling  in  the  inertial  sub-range  and 
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below)  at  the  small  end,  and  a  mixing  length  profile  of  the  mesoseale  model  (appropriate 
for  all  scales  of  turbulence)  at  the  large  end. 

4,  Conclusions 

As  the  resolution  of  mesoseale  NWP  models  approaches  the  turbulence  scales,  the 
appropriateness  the  mixing  length  used  in  ensemble  parameterization  becomes 
questionable.  The  filter  scale  applied  here  to  the  calculation  of  the  integral  length  scale  is 
analogous  to  twice  the  horizontal  grid  resolution  of  NWP  models.  For  example,  a  filter 
scale  of  1  km  gave  an  appropriate  integral  length  scale  representing  eddies  scaling  1  km 
or  less.  A  high-resolution  NWP  model  using  a  horizontal  grid  resolution  of  500  m  should 
parameterize  all  turbulence  on  the  SGS,  which  is  technically  1  km  for  this  resolution. 
Thus,  as  discussed  in  chapter  II,  the  mixing  length  used  in  the  parameterization  of  the 
SGS  eddies  should  also  represent  the  1  km  and  smaller  eddies,  and  should  not  use  a 
mixing  length  representing  the  turbulent  eddies  from  all  scales.  Examination  of  Figure 
3.9  provides  the  justification  for  this  argument,  as  well  as  providing  guidance  in  adjusting 
the  length  scale  within  the  COAMPS  model  to  appropriately  consider  the  scales  of 
turbulence  that  are  being  parameterized  when  the  horizontal  resolution  in  and  below  the 
range  of  scale  of  the  largest  turbulent  eddies.  The  formulation  derived  (Equation  3.3) 
will  be  needed  to  adjust  to  the  current  COAMPS  MY82  mixing  length  formulation  and 
will  be  tested  in  COAMPS  in  chapter  V. 
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IV.  CASE  STUDY  USING  HIGH-RESOLUTION  COAMPS™ 

A.  INTRODUCTION 

Using  COAMPS^'^,  we  simulated  a  cold  air  outbreak  (CAO)  case  during  the  JES 
experiment  in  which  extensive  roll  vortices  were  observed.  In  this  chapter,  we  will 
introduce  this  case  and  present  the  simulated  COAMPS™  results.  The  objective  of  the 
initial  COAMPS^'^  simulation  of  this  case  study  is  to  analyze  the  model’s  performance  at 
high-resolution  (500  m).  We  will  compare  the  high-resolution  model  output  fields  with 
satellite  imagery,  aircraft  observations,  and  results  of  past  studies  on  similar  weather 
phenomena,  to  determine  the  model’s  ability  to  predict  the  state  of  the  atmosphere,  on 
both  the  mesoscale  and  turbulent  scales.  This  chapter  also  discusses  the  effects  of  model 
horizontal  diffusion  and  advection  on  the  model  results.  The  results  of  the  initial  case 
study  will  serve  as  the  baseline  simulation  (or  control  run)  with  which  we  can  compare 
the  results  from  a  modified  COAMPS™  simulation. 

B,  COAMPS™ 

The  mesoscale  model  used  in  this  research  is  in  the  Navy's  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System  (COAMPS™)  developed  at  the  Naval 
Research  Laboratory  (NRL),  Monterey,  CA.  COAMPS™  is  described  in  detail  in  Hodur 
(1997)  and  Hodur  and  Doyle  (1998).  The  model  includes  a  data  assimilation  system,  non¬ 
hydrostatic  atmosphere  model  components,  a  hydrostatic  ocean  circulation  model,  and  an 
ocean  wave  model.  The  atmospheric  and  ocean  models  can  be  used  separately,  or  in  a 
fully  coupled  mode  (still  in  development)  which  allows  surface  fluxes  of  heat, 
momentum,  and  moisture  to  be  exchanged  across  the  air-ocean  interface  every  time  step. 

The  COAMPS™  atmospheric  model  is  currently  used  operationally  by  U.S.  Navy 
forecasters  for  short-term  regional  forecasting  (up  to  72  hours),  and  is  designed  for  both 
idealized  and  real  data  simulations.  For  real  data  simulations  (performed  in  this  study), 
COAMPS™  uses  a  complete  data  assimilation  system  comprised  of  data  quality  control, 
analysis,  initialization,  and  forecast  model  components.  The  model  grid  projection  is 
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specified  along  with  the  latitude  and  longitude  of  any  one  point  in  the  grid,  making 
COAMPS™  globally  locatable.  The  model  utilizes  nested  grids,  allowing  any  number  of 
meshes,  with  a  3:1  ratio  reduction  in  grid  spacing  between  the  grids.  Inner  grids  can  be 
specified  arbitrarily  within  the  confines  of  the  next  coarser  grid,  allowing  the  highest 
resolution  to  be  focused  over  a  specific  region  of  interest. 

The  COAMPS™  data  analysis  routine  produces  a  set  of  initial  conditions  for  the 
forecast  model.  The  analysis  specifies  3D  atmospheric  fields  and  lateral  boundary 
conditions  from  the  interpolation  of  global  atmospheric  model  analysis  or  forecast  fields 
generated  by  the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  to 
the  COAMPS™  grid(s).  The  interpolated  COAMPS™  analysis  fields  for  the  u  and  v 
wind  components  and  geopotential  height  are  blended  with  real-time  observations  (such 
as  rawinsonde,  SSMI,  satellite-derived  data,  and  aircraft  observations),  if  available,  using 
3 -dimensional  multi-variate  optimum  interpolation  (MVOI)  analysis.  These  model  and 
observation  fields  are  available  from  the  Fleet  Numerical  Meteorology  and 
Oceanography  Center,  Monterey,  CA.  In  addition  to  the  initial  atmospheric  fields,  several 
surface  conditions  are  set  with  various  fields.  Surface  terrain  height  is  obtained  from  the 
Defense  Mapping  Agency  (DMA)  Digital  Terrain  Elevation  Data  (DTED)  level  1  data  (1 
km  resolution),  and  is  bilinearly  interpolated  to  the  COAMPS™  grid(s).  Surface  albedo, 
surface  roughness,  ground  temperature,  and  ground  wetness  are  also  specified  initially. 
Monthly  climatological  global  fields  of  albedo  and  surface  roughness  are  bilinearly 
interpolated  to  the  COAMPS™  grid(s).  Albedo  is  set  to  0.09  over  the  water.  Ground 
wetness  and  ground  temperature  are  specified  from  a  previous  COAMPS™  forecast  if 
available.  SST  is  set  by  bilinearly  interpolating  the  Elect  Numerical  Meteorology  and 
Oceanography  Center  (ENMOC)  global  SST  analysis  to  the  COAMPS™  grid(s). 

The  COAMPS™  analysis  is  performed  on  the  Arakawa-Eamb  scheme  A  grid  (no 
staggering),  while  both  the  horizontal  and  vertical  grids  in  the  forecast  model  are 
staggered.  Scalar  forecast  variables  are  defined  at  the  center  of  the  model  grid  box,  with 
the  forecast  wind  fields  defined  at  half  the  grid  distance  from  the  center  of  the  model  grid 
box.  The  analysis  fields  are  interpolated  to  the  initial  forecast  fields  using  a  bicubic 
spline  to  the  Arakawa  grid  C  grid  configuration  (staggering)  in  the  horizontal,  then 
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interpolated  in  the  vertical  to  the  model  vertical  coordinate.  The  vertical  coordinate  in 
COAMPS™  is  a  terrain-following  sigma-z  system  (Gal-Chen  and  Sommerville  1975). 

The  equations  used  in  the  COAMPS™  atmospheric  model  are  comprised  of  a 
non-hydrostatic,  compressible  formulation  of  the  primitive  equations,  and 
parameterizations  for  the  sub-grid  scale  mixing,  surface  fluxes,  explicit  moist  physics, 
cumulus  convective  processes,  and  radiation.  The  equations  are  solved  in  three 
dimensions  with  the  terrain-following  sigma-z  vertical  coordinate.  The  complete  list  of 
prognostic  equations  and  their  origination  is  found  in  Hodur  (1997).  The  numerical 
scheme  used  is  the  leapfrog  (centered-in-time,  centered-in-space)  second-order 
differencing.  The  PBL  and  free  atmospheric  turbulent  mixing  and  diffusion  are  modeled 
after  MY82  as  described  in  chapter  II. 

C.  JES  CASE  STUDY  -  DESCRIPTION 

1,  Overview 

In  this  section  we  describe  the  weather  event  that  will  be  simulated  in  our  study. 
The  entire  modeling  region  (Figure  4.1)  expands  to  the  west  over  North  and  South  Korea 
and  eastern  China,  to  the  north  over  part  of  southern  Russia,  and  to  the  east  and  south 
over  the  larger  part  of  the  JES.  The  small  region  of  interest  that  will  enclose  the  high- 
resolution  domain  of  our  simulation  (outlined  in  Figure  4.1)  is  centered  over  the  western 
JES  south  of  Vladivostok,  located  at  approximately  41  to  42  degrees  north  latitude  and 
132  to  133  degrees  east  longitude,  enclosing  about  1000  km^.  The  coastal  mountain 
ranges  west  of  the  JES,  the  Hamgyong  and  Taeback  Mountains  in  Korea  and  the  Sikhote- 
Alin  Range  in  Russia,  are  also  seen  clearly  in  Figure  4.1.  A  summary  of  the  synoptic 
situation  and  resulting  mesoscale  weather  features  for  the  case,  using  descriptions  from 
Dorman  et  al.  (2002),  is  given  below.  A  discussion  of  convective  roll  vortices,  the  main 
weather  phenomenon  of  this  case,  is  also  presented. 

2.  Synoptic  Conditions  for  JES0216 

On  February  16  and  17,  2000,  a  Very  Cold  Siberian  Air  Outbreak  (VCSAO) 
occurred  over  the  JES.  As  described  in  Dorman  et  al.  (2002),  these  events  occur  just  a 
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Figure  4.1,  Region  of  model  simulation,  centered  over  the  JES.  The  high-resolution 
domain  of  the  simulation,  the  region  of  interest,  is  outlined  in  the  figure.  The  coastal 
mountain  ranges  are  seen  as  the  shaded  red  regions. 
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few  days  a  month  during  the  winter  over  the  JES.  A  VCSAO  oceurs  when  the  Siberian 
High  moves  farther  south  of  its  mean  position,  generating  northerly  winds  with  -30°  to  - 
20°  C  air  over  the  Russian  Coast  and  JES,  initiating  a  strong  CAO.  The  The  Navy 
Operational  Global  Atmospherie  Predietion  System  (NOGAPS™)  analysis  fields  of 
pressure  and  1000  to  500  mb  thiekness  for  16  Eeb  2000  1200Z  and  17  Eeb  2000  OOOOZ 
(Eigure  4.2)  and  the  vertical  p-velocity  (omega)  at  850  mb  for  17  Eeb  2000  OOOOZ 
(Eigure  4.3b)  show  that  the  Siberian  high  is  indeed  located  farther  south  into  China, 
directly  west  of  the  JES,  creating  subsidence  (indicated  by  the  positive  omega  values) 
over  our  small  region  of  interest.  NOGAPS^*^  analysis  fields  also  show  that  the  cold  front 
associated  with  the  low  pressure  center  northeast  of  the  JES,  and  strongest  cold  air 
advection,  has  already  moved  east  of  the  JES,  indicated  by  the  strong  gradient  of 
thickness  contours  in  this  region.  Even  though  the  thickness  gradient  is  now  weaker  over 
our  region  of  interest,  very  cold  air,  indicated  by  the  low  thickness  values,  continues  to  be 
advected  eastward  over  the  JES.  The  NOGAPS™  analyses  in  Eigure  4.2  also  show  that 
the  synoptic  conditions  have  not  changed  much  in  the  12  hours  between  these  analyses. 

Ten-day  composite  SST  imagery  (1  to  10  Eebruary  2000)  from  the  National 
Ocean  and  Atmospheric  Administration’s  polar  orbiting  satellite  number  15  (NOAA-15) 
Advanced  Very  High-resolution  Radiometer  (AVHRR)  and  the  NOGAPS™  SST 
analysis  (Eigure  4.4)  show  the  structure  and  gradient  of  the  SST  field  for  this  case.  The 
NOGAPS™  SST  analysis  was  uniformly  reduced  2  °C  to  better  match  the  magnitude  of 
observed  SST  in  our  region  of  interest  in  the  western  JES.  This  corrected  SST  analysis 
will  be  used  in  the  initialization  of  our  COAMPS  simulations.  The  satellite  imagery  and 
NOGAPS  analysis  show  near  0  °C  SST  over  our  small  area  of  interest,  while  a  stronger 
NW  to  SE  SST  gradient,  from  colder  to  warmer  temperatures,  exists  over  the 
southeastern  JES.  The  equivalent  potential  temperature  at  1013  mb  (Eigure  4.3a) 
indicates  that  the  near-surface  air  temperature,  between  258  and  264  K  in  our  region  of 
interest,  is  8  to  10  °C  cooler  than  the  SST.  The  large  air-sea  temperature  difference  is  the 
driving  force  in  the  down-wind  development  of  the  turbulent  BE  over  the  JES  for  this 
case,  as  the  BE  air  warms  as  it  moves  over  the  relatively  warmer  ocean,  generating  large 
surface  fluxes  of  heat  and  moisture. 
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Figure  4,2,  NOGAPS  Global  Analysis  Fields,  Valid  (a)  16  Feb  2000,  1200Z  and  (b)  17 
Feb  2000,  OOOOZ,  showing  mean  sea  level  pressure  in  millibars  (solid)  and  geopotential 
thickness  in  meters  (dashed).  The  region  of  interest  for  the  high-resolution  domain  is 
outlined. 
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Figure  4,3,  Same  as  in  Figure  4.2,  exeept  for  (a)  near-surfaee  equivalent  potential 
temperature  (K)  and  (b)  omega  (mbs"').  The  region  of  interest  for  the  high-resolution 
domain  is  outlined. 
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Figure  4,4,  SST  (in  degrees  C)  variation  over  the  JES.  (a)  NOGAPS  SST  Field  for  16 
February  2000,  1200  Z.  (b)  10-day  eomposite  (1  to  10  February  2000)  SST  derived 
from  NOAA  AVHRR  Satellite  Imagery.  The  region  of  interest  for  the  high-resolution 
domain  is  outlined. 
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Dorman  et  al.  (2002)  discusses  the  differences  between  typical  CAO  conditions  and 
VCSAO  conditions  over  the  JES.  In  typical  CAO  conditions,  a  strong  surface  inversion 
and  weaker  northerly  winds  trap  the  eoldest  inland  air  to  the  west  of  the  Russian  eoastal 
mountains.  Reduced  drag  over  the  JES  causes  surface  divergence,  which  pulls  air  down 
the  lee  of  the  mountains  and  adiabatically  warms  the  coastal  air  such  that  the  air-sea 
temperature  difference  is  only  around  5°C.  In  VCSAO  conditions,  the  much  colder 
surface  air,  which  reduces  stability  and  eliminates  the  surface  inversion,  coupled  with 
stronger  winds,  allows  much  colder  air  to  advect  over  the  mountains  and  out  over  the  JES 
compared  to  the  typical  CAO,  creating  much  larger  air  sea  temperature  differences  and 
surface  heat  fluxes.  These  differenees  between  typieal  CAO  and  VCSAO  eonditions 
result  in  the  larger  magnitude  of  turbulence  and  more  extensive  and  detailed  formation  of 
heat-flux  driven  horizontal  roll  vortices  and  cloud  streets  in  the  VCSAO.  The  VCSAO 
case  of  17  Eebruary  2000,  our  case  study,  will  be  referred  to  from  here  on  as  JES0216. 

3.  BL  Roll  Vortices 

a.  Satellite  Imagery 

We  can  see  the  down-wind  development  of  the  BE  over  the  JES  on  17 
Eebruary  2000  from  the  cloud  structure.  Visible  satellite  imagery  from  the  Japanese 
Meteorologieal  Ageney  (JMA)  geostationary  satellite  (Eigure  4.5)  shows  80%  of  the 
Japan/East  Sea  covered  in  convective  roll  elouds  and  stratus.  The  elouds  appear  to  be 
more  roll-like  in  the  western  JES,  in  the  region  of  stronger  large-seale  subsidence  shown 
in  the  omega  field  (Eigure  4.4),  and  more  cellular  in  the  eastern  JES  near  Japan,  where 
omega  at  850  mb  is  negative  and  the  large-scale  subsidence  is  absent  in  the  presence  of 
the  upper-level  trough  (not  shown).  There  is  a  line  of  thin  clouds  streaming  off  the 
Russian  eoast  and  thorough  the  eenter  of  our  region  of  interest  (outlined  in  Eigure  4.5), 
while  the  thickest  clouds  form  farther  off  shore  where  there  is  more  BE  turbulence  to 
generate  large  eddies  and  moisten  the  originally  dry  eontinental  air  to  the  point  where 
low-level  convective  clouds  develop  at  the  top  of  the  BE  below  the  inversion.  Roll 
vortices  were  observed  over  the  JES  on  several  days  during  the  JES  experiment  (Dorman 
et  al.  2002)  in  addition  to  the  JES0216  case. 
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Figure  4,5,  JMA  visible  satellite  imagery  over  the  JES  on  February  17,  2000,  0200  Z. 
Thin  cloud  streets  are  present  within  the  region  of  interest  (outlined). 
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b.  General  Structure  and  Characteristics  of  Roll  Vortices 

Roll  vortices  are  frequently  observed  in  the  marine  BL,  especially  in  CAO 
conditions.  Etling  and  Brown  (1993)  and  Young  et  al.  (2002)  provide  reviews  of  the 
convective  roll  vortices.  They  loosely  define  roll  vortices  as  quasi  two-dimensional  (i.e., 
nearly  linear)  organized  large  eddies,  extending  through  the  entire  depth  of  the  PBL, 
with  properties  different  from  those  of  small-scale  turbulence.  Oecurring  usually  in  high 
wind  conditions,  these  features  exhibit  coherent  perturbations  in  the  wind,  temperature, 
and  humidity  fields  (e.g.,  LeMone  and  Pennell  1976).  The  rolls  eontribute  significantly  to 
the  vertieal  fluxes  of  heat,  moisture,  and  humidity  from  the  surfaee  through  the  entire 
depth  of  the  PBL,  and  thus  play  a  signifieant  role  BL  in  mixing,  transport,  and 
entrainment  processes. 

A  common  feature  assoeiated  with  roll  vortices,  but  not  always  present,  is 
the  manifestation  of  so-ealled  cloud  streets  that  are  easily  seen  in  satellite  imagery 
(Ligure  4.5).  If  adequate  moisture  is  present  for  their  formation,  elouds  streets  can  form 
in  the  updraft  regions  of  the  roll  circulations,  while  eloud-free  areas  exist  in  between  the 
streets  in  downdraft  regions.  The  rolls  and  associated  cloud  streets  are  usually  oriented 
more  or  less  parallel  to  the  mean  PBL  wind  direetion.  Cloud  streets  are  most 
impressively  formed  during  eold-air  outbreaks  over  the  oceans,  but  ean  form  over  land  as 
well.  Ligure  4.6  shows  a  simple  schematic  of  horizontal  roll  vortices  structure  and  the 
resulting  cloud  streets  in  the  PBL. 

Etling  and  Brown  (1993)  cite  many  studies  that  have  been  condueted,  for 
both  over  land  and  over  oeean  cases,  to  investigate  the  vertieal  and  horizontal  structure  of 
roll  vortiees.  The  results  can  be  summarized  as  follows: 

Vertical  Extent  (H)  1-2  km 

Wavelength/horizontal  scale  (k)  2-20  km 

Aspect  ratio  (k/H)  2-15 

Down-wind  extent  (E)  10-1 000  km 
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Orientation  of  roll  axis  to  the  mean  wind  (s)  -20  to  +30  deg 

Life  time  (x)  1-72  h 

The  rolls  have  a  vertieal  seale  roughly  equal  to  the  BL  depth,  and  a  wide  range  of 
horizontal  scales,  leading  to  a  wide  range  of  aspect  ratios.  Observational  and  LES 
modeling  results  summarized  by  Young  et  al.  (2002)  give  the  aspect  ratio  as  a  function  of 
BL  depth.  The  results  are  somewhat  varied,  but  for  CAO  cases  with  BL  depths  near  800 
m  (similar  to  JES0216),  aspect  ratios  between  4  and  7  are  most  common.  High- 
resolution  observations  indicate  the  presence  of  multiple  thermals  within  the  roll  updraft 
regions  (Brummer  1999).  Typically,  two  or  three  stronger  updraft  “cores”,  rather  than  a 
single  updraft,  are  located  within  a  region  of  overall  weak  roll  updrafts,  leading  to  a 
rather  frothy  appearance  in  satellite  imagery  for  any  associated  cloud  streets,  as  is 
apparent  in  Ligure  4.5  for  JES0216. 

Conditions  for  the  formation  of  these  rolls  has  been  the  subject  of  many 
observational  studies  (see  Brown  1980).  The  cause  of  the  rolls  has  been  attributed  to 
both  shear  instability  of  the  Ekman  BL  (inflection  point  instability)  and  convective 
instability,  but  typically  a  combination  of  these  two  instability  mechanisms  generates  the 
roll  structure.  However,  it  has  been  shown  by  Brummer  and  Latif  (1985)  and  by  Etling 
and  Raasch  (1987)  that  inflection  point  instability  is  suppressed  when  the  height  of  the 
inflection  in  the  BL  wind  is  greater  than  6zi.  In  general,  theory  suggests  that  roll  vortices 
occur  for  all  convection  regimes  in  the  presence  of  little  shear,  and  for  neutral  and  even 
slightly  stable  stratification  in  the  rotating  frame  of  reference.  Nearly  all  field  studies  of 
roll  vortices  have  been  conducted  during  CAO’s  over  the  ocean,  or  in  the  daytime  PBL 
over  land,  which  implies  that  the  rolls  are  related  to  unstable  stratification  and  convective 
instability.  Thus,  these  large  eddies  routinely  form  in  an  unstably  stratified  PBL  heated  at 
the  surface. 


The  parameter, 

L  0  ul 


(4.1) 
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where  L  is  the  Monin-Obukov  length  scale  for  the  surface  layer,  is  often  used  as  a 
measure  of  PBL  stratification,  and  has  become  the  criterion  measurement  for  the  onset  of 
roll  convection.  Numerical  and  observational  studies  of  convective  BLs  (see  Etling  and 
Brown  1993)  related  types  of  thermal  convection  to  Zi/L  as  a  stratification  parameter  with 
a  negative  z/L  indicating  thermal  instability.  Studies  conclude  in  general  that  for  -5< 
Z./L  <  -1.5,  convection  with  roll  vortices  may  exist,  but  for  Zi/L  <  -25,  three-dimensional 
convection  is  more  likely.  Hence  the  roll-regime  exists  in  the  unstably  stratified  regime,  - 
25  <  Zi/L  <  -1.5,  but  Etling  and  Brown  (1993)  state  that  roll  vortices  and  cloud  streets  will 
most  likely  be  observed  for  moderately  unstable  stratification,  -15  <  Zi/L  <  -5,  with  Zi/L 
usually  less  (more)  negative  for  moderate  (low)  winds  and  moderate  (strong)  surface 
heating. 

Because  of  the  quasi  two-dimensional  structure  of  rolls,  many  two- 
dimensional  modeling  studies  of  these  features  have  been  employed.  The  first  modern 
numerical  attempt  at  simulating  roll  development  was  due  to  Mason  and  Sykes  (1980) 
who  used  a  finite  difference  model  to  simulate  the  development  of  roll  vortices  for  the 
neutrally  stratified  BE.  These  authors  extended  their  work  to  investigate  various  aspects 
of  two-dimensional  roll  circulations  in  the  unstably-stratified  BE  (Mason  and  Sykes 
1982;  Mason  1983;  Mason  1985).  Three-dimensional  modeling  of  roll  vortices  has 
mostly  been  limited  to  EES  modeling  of  the  PBL.  Numerical  studies  by  Chlond  (1992) 
and  others  cited  in  Etling  and  Brown  (1993)  indicate  that  large-eddy  simulations  of  cases 
with  convective  instability  produce  roll  vortices  in  the  PBL.  However,  in  contrast  to  2-D 
numerical  results  and  observations  of  PBL  rolls,  they  do  not  show  up  clearly  as 
permanent  features  (showing  up  at  early  integration  times,  then  disappearing  later  on). 
The  lack  of  permanency  may  be  attributed  to  the  limited  domain  sizes  used  in  these  EES 
experiments.  Modeling  of  roll  vortices  using  high-resolution  mesoscale  models  has  not 
yet  been  documented. 

Complicating  a  simple  picture  of  these  two-dimensional  rolls  are  other 
smaller  scale  PBL  eddies,  whose  energy  exchanges  with  the  rolls  may  be  important.  As 
was  illustrated  in  chapter  III,  partitioning  of  PBL  vertical  fluxes  due  to  different  scales  of 
motion,  i.e.  small  scale  turbulence  and  larger  eddies,  may  be  done  using  spectral  analysis 
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of  observations  or  model  results.  Spectra  of  aircraft  measurements  during  the  GALE 
experiment  analysis  performed  by  Chou  and  Zimmerman  (1989)  illustrate  that  the  mean 
vertical  fluxes  of  heat,  momentum,  and  moisture  can  be  attributed  to  vertical  motions  on 
different  scales:  small  scale  turbulence  <  200  m),  thermal  plumes  (200  m  <  Z<  2  km), 
and  roll  vortices  (k  >  2  km).  The  energy  exchange  between  large  and  small  scale 
turbulent  eddies  implies  the  importance  of  correctly  estimating  the  turbulent  scale  energy 
in  a  mesoscale  NWP  model  in  order  to  correctly  predict  the  roll  structure  occurring  on  the 
mesoscale. 

c.  Aircraft  Observations 

Aircraft  observations,  from  level  legs  near  the  surface  and  BE  soundings, 
were  taken  within  our  small  region  of  interest  during  the  JES  experiment  on  17  Eebruary 
2000  OOOOZ.  Three  level  flight  legs,  each  flown  at  a  height  of  about  40  meters,  and  ten 
soundings  through  the  depth  of  the  BE  were  flown  within  this  region.  Specific  times  and 
exact  locations  of  the  flight  legs  and  soundings  can  be  found  in  Table  3.1.  The  general 
locations  of  the  flight  legs  and  soundings  within  the  region  are  shown  in  Eigure  4.7.  The 
three  level  flight  legs  were  flown  in  the  East-West  direction,  located  in  northern,  central, 
and  southern  positions  within  the  region.  These  three  legs  were  flown  within  one  hour, 
between  0100  Z  and  0200  Z.  The  high-rate  (10  Hz)  aircraft  observations  from  the  near¬ 
surface  legs  are  filtered  here  to  remove  variations  greater  than  0.1  Hz  for  visual 
presentation.  The  BE  soundings  were  made  between  heights  of35mto917m  above  sea 
level.  Eight  of  the  sounding  legs  were  flown  in  a  zig-zag  pattern  in  the  vertical,  forming 
a  vertical  cross-section  in  the  along-wind  direction  (see  Eigure  4.10).  These  observations 
reveal  the  structure  of  the  BE  on  this  day  in  more  detail. 

Eigure  4.8  shows  the  time  series  of  vertical  velocity,  absolute  humidity, 
and  potential  temperature  from  the  near-surface  flight  legs.  With  the  level  flight  legs 
oriented  at  about  a  45  degree  angle  to  the  mean  wind  direction,  the  surface  layer 
observations  reveal  both  a  signal  of  the  down-wind  development  of  the  BE  while  at  the 
same  time  traversing  across  the  updraft  and  downdraft  regions  of  the  turbulent  rolls.  As 
expected,  there  is  a  general  moistening  and  warming  of  the  air  mass,  both  from  north  to 
south  and  east  to  west.  There  is  a  well-defined  oscillation  in  the  temperature  and  water 
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Figure  4,7.  (a)  General  loeation  of  aircraft  observation  legs  and  soundings  within  the 
region  of  interest  over  the  western  JES  (outlined  region  in  figures  4.1  through  4.5). 
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Figure  4,8.  From  left  to  right:  Observed  vertical  velocity  (m  s"'),  water  vapor  (g  kg  -1) 
and  potential  temperature  (K)  from  near-surface  (40  m  above  sea  level)  flight  legs  from 
north  (R06)  to  south  (R02)  in  the  region  of  interest  in  the  western  JES.  Observations 
taken  17  Feb  2000,  OlOOZ  -  0200Z.  Horizontal  axes  are  distance  (D,  in  km)  from  the 
beginning  of  the  leg.  See  Figure  4.7  for  the  location  of  the  flight  legs. 
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vapor  traces,  which  shows  the  presence  of  the  roll  vortices.  The  peaks  of  higher 
temperature  and  water  vapor  are  collocated,  as  expected,  with  the  peaks  in  vertical 
velocity,  as  warm,  moist  air  is  transported  upward  from  the  ocean  surface  within  updrafts. 
The  vertical  velocity  is  also  seen  to  become  more  variable  and  larger  in  magnitude  from 
north  to  south  as  the  BL  develops  and  becomes  more  turbulent,  with  stronger  and  more 
clearly  defined  updraft  and  downdraft  regions.  We  also  observe  the  stronger  updraft 
cores  within  regions  of  weaker  upward  vertical  velocity,  as  is  commonly  observed  in 
cases  of  roll  convection  (Young  et  al.  2002).  The  length  scale  of  the  rolls  can  be  seen 
most  clearly  in  the  western  half  of  the  central  and  southernmost  legs,  where  the 
horizontal  distance  between  successive  updraft  regions  and  peaks  in  moisture  is  5-7  km. 
Accounting  for  the  aircraft  flight  path  with  respect  to  the  mean  horizontal  wind  direction 
and  roll  orientation,  the  actual  scale  of  the  rolls  is  likely  closer  to  3-5  km.  Assuming  the 
BL  depth  is  around  850  m,  based  on  the  maximum  height  of  the  BL  soundings,  the  roll 
aspect  ratio  is  4  to  6.  These  values  of  roll  scale  and  aspect  ratio  are  in  agreement  with 
cases  of  roll  convection  as  described  earlier. 

Figure  4.9  reveals  more  typical  BL  conditions  associated  with  a 
convective  roll  vortex  regime.  The  Figure  presents  the  time  series  of  surface  layer  wind, 
air-sea  temperature  difference,  and  combined  sensible  and  latent  heat  fluxes  from  the 
three  near-surface  level  flight  legs.  The  flight  level  winds  range  from  13  to  16  m  s'^  and 
air-sea  temperature  differences  in  the  western  JES  are  around  8.5  degrees,  resulting  in 
surface  heat  fluxes  (sensible  and  latent  heat)  of  up  to  250  W  m'  .  We  see  that  the 
combination  of  higher  wind  and  large  air-sea  temperature  difference  creates  the  largest 
heat  fluxes  along  the  northernmost  leg.  The  decrease  in  heat  flux  from  north  to  south 
appears  to  be  dominated  by  the  corresponding  decrease  in  wind  speed.  The  vertical 
cross-sections  in  Figure  4.10  are  derived  from  aircraft  soundings  in  the  intensive 
measurement  area.  The  cross-sections  show  clearly  that,  in  the  down-wind  direction  (left 
to  right  in  the  cross-section),  the  BL  is  indeed  undergoing  slight  warming  (Figure  4.10  b) 
and  moistening  (Figure  4.10  a),  as  was  also  shown  in  the  near-surface  observations.  It  is 
also  evident  that  the  BL  is  becoming  moister  at  higher  levels  down-wind.  As  a  result  of 
turbulent  mixing  in  the  boundary  layer,  the  initially  dry  air  originating  from  over  land 
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Figure  4,9.  Same  as  Figure  4.8,  exeept  for  (from  left  to  right)  horizontal  wind  speed  (m 
s'^),  air-sea  temperature  difference  (K),  and  Total  (latent  and  sensible)  heat  flux. 
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Figure  4,10,  BL  cross-sections  in  the  up-wind  direction  (along  the  flight  track)  of  (a) 
water  vapor  (g  kg“')  and  (b)  potential  temperature  (K)  derived  from  eight  sounding  legs. 
The  sounding  flight  track  is  shown  within  the  contour  plots.  See  Figure  4.7  for  the 
loeations  of  the  soundings. 
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gradually  moistens  reaches  to  higher  and  higher  levels  with  down-wind  distance  within 
the  BL,  where  the  roll  clouds  eventually  form  in  the  updraft  regions.  In  addition,  Figures 
4.10  a  and  b  reveal  that  the  aircraft  soundings  did  not  reach  the  inversion  layer  normally 
above  the  boundary  layer  since  there  is  no  inversion  present  from  the  soundings.  Thus, 
the  boundary  layer  top  is  above  900  m.  The  inversion  is  probably  only  slightly  higher 
than  900  m,  since  the  experiment  soundings  are  intended  to  cover  the  entire  depth  of  the 
BL  (Kheilf  et  al.  2003).  Assuming  a  BL  depth  of  950  m,  and  taking  the  mean  potential 
temperature  of  the  BL  from  the  JES  soundings,  the  surface  layer  sensible  heat  flux,  and 
u*  calculated  from  the  surface  layer  flux  legs,  the  average  Zi/L  in  the  BL  for  our  region  of 
interest  is  -  5.1  (Equation  4.1).  This  value,  predicting  a  moderately  unstable  regime,  falls 
within  the  range  of  z/E  given  in  the  literature  for  convective  roll  development. 

D,  DESCRIPTION  OF  COAMPS™  SIMULATIONS  OVER  THE  JES 

We  performed  a  high-resolution  COAMPS™  simulation  over  the  region  of 
interest  in  the  JES  and  for  the  JES0216  case  described  in  the  previous  section.  The 
modeling  region  consisted  of  4  nested  domains.  The  horizontal  resolutions,  sizes  for  each 
domain,  and  simulation  times  for  the  simulation  are  described  in  Table  4.1,  and  the 
locations  of  the  domains  are  shown  in  Eigure  4.11.  The  geographic  locations  of  the  outer 
domain  (domain  1)  was  chosen  such  that  the  domain  extended  westward  over  the  coastal 
mountains  and  farther  to  the  west  over  land,  such  that  topographic  effects  upwind  would 
be  simulated  and  modeled  within  the  inner  nests.  The  high-resolution  inner  domain,  at  a 
resolution  of  500  m,  was  limited  in  size  by  computational  constraints,  and  is  the  primary 
focus  of  this  study.  Unless  otherwise  specified,  all  simulation  results  presented  will  be 
from  the  JES0216  high-resolution  innermost  domain.  The  aircraft  flight  paths  and 
locations  of  both  surface  layer  flux  legs  and  BL  soundings  within  the  inner  0.5  km 
resolution  nest  shown  earlier  (Eigure  4.7)  will  be  used  in  comparison  to  model  output. 
The  simulation  start  time,  12  to  16  hours  prior  to  the  time  of  the  observations,  was  chosen 
to  allow  for  adequate  spin-up  time  for  the  model.  Eor  this  simulation,  the  model 
contained  30  vertical  levels,  of  which  12  were  located  below  1  km  for  better  resolution  of 
the  marine  BE.  The  simulation  involved  calculations  at  30  second  time  steps  over  16 
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Table  4.1.  Domain  Set  Up  for  COAMPS™  JES0216 


Domain  1  (outer) 
Resolution  (km)  /  size  (km^) 

Domain  2 

Resolution  (km)/  size  (km^) 

Domain  3 

Resolution  (km)/  size  (km^) 

Domain  4 
(innermost) 
Resolution 
(km)/size  (km^) 

13.5/  1107x  1107 

4.5  /585  X  585 

1.5  /272X  272 

0.5/ 100  X  100 

Figure  4.11.  Four  nested  COAMPS™  domains  over  the  JES  region.  Observation 
legs  used  for  model-observation  inter-comparison  (shown  in  Figure  4.7)  are  within  the 
innermost  domain. 
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hours  for  2.88  x  10^  grid  points,  which  required  a  wall  elock  simulation  time  of  over  120 
hours  running  on  32  processors  on  a  CRAY  SVlex  supereomputer  at  the  Arctie  Regional 
Supereomputing  Center  in  Anehorage,  Alaska. 

The  persistence  of  the  synoptie  eonditions  from  16  February  1200Z  and  17  February 
OOOOZ  (Figure  4.2)  imply,  assuming  adequate  spin-up  time  (typically  around  6  to  9  hours 
for  COAMPS™  (Hodur  1997),  that  a  model  foreeast  initialized  from  the  16  February 
1200Z  analysis  fields  of  this  region  will  be  a  valid  representation  of  the  development  of 
the  BL  turbulent  structure.  It  was  determined  by  inspection  of  the  model  forecast  fields 
throughout  the  forecast  period  that  the  12-hour  forecast  best  represents  the  turbulent  BL 
structure.  Thus  we  will  present  the  simulated  fields  from  the  12-hour  foreeast  in  our 
discussion  and  in  comparison  to  aircraft  observations.  The  observational  analysis  will  be 
based  on  the  limited  set  of  observations  available  for  comparison  to  the  model  results,  so 
we  will  also  rely  on  past  modeling  and  observational  studies  to  supplement  our  analysis 
of  our  case  study. 


E.  ANALYSIS  OF  COAMPS™  SIMULATIONS 

1.  Effects  of  Model  Numerics 

Before  we  begin  the  analysis  of  the  model  output  from  the  JES0216  case,  a 
diseussion  on  the  issue  of  model  resolvability  is  in  order.  As  discussed  in  Pielke  (2002) 
and  Avissar  et  al.  (1990),  the  energy  easeade  of  turbulence  to  smaller  and  smaller 
wavelengths  (described  in  chapter  II)  cannot  occur  in  mesoscale  numerical  models. 
When  wave  interactions  occur  in  the  real  world,  a  speetra  of  wavelengths  result. 
Eventually  the  smallest  eddies  are  removed  by  moleeular  dissipation.  In  a  numerical 
model,  however,  the  existence  of  a  diserete  grid  results  in  aliasing  in  the  resolvable  scale 
phenomena  due  to  perturbations  in  the  SGS  fields  (also  referred  to  as  short  waves). 
These  erroneously  produeed  larger  seale  waves  again  interact  to  produce  larger  and 
smaller  waves.  Because  the  proper  cascade  of  energy  to  smaller  and  smaller  scales  is 
disrupted,  a  fictitious  energy  buildup  oceurs  in  the  model.  Pielke  (2002)  states  that  sueh 
short  waves  are  badly  handled  and  poorly  represented  in  amplitude  and/or  phase  by 
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various  numerical  schemes  on  a  computational  grid,  even  for  linear  scale  advection. 
Thus,  the  solution  adopted  to  control  the  shortest  wavelengths  is  to  damp  or  remove  these 
waves  using  numerical  filters.  In  COAMPS™  and  other  mesoscale  NWP  models,  an 
explicit  diffusion  equation  is  used  to  parameterize  the  sub-grid  scale  correlation  terms  so 
that  energy  is  extracted  from  the  averaged  equations  in  a  manner  consistent  with  reality. 
In  COAMPS™,  a  background  second-order  (for  points  on  the  grid  boundaries)  and 
fourth-order  (for  all  other  grid  points)  diffusion  of  all  prognostic  variables  (except  for 
pressure  perturbation)  is  carried  out  on  sigma  surfaces.  Hodur  (1997)  describes  the  form 
of  the  diffusion  operator  and  criteria  for  choice  of  diffusion  coefficients.  Unfortunately, 
as  discussed  in  chapter  II,  only  the  vertical  sub-grid  scale  correlation  terms  are  reasonably 
well  known  and  can  be  parameterized  accurately  in  terms  of  dependent  variables.  In 
contrast,  horizontal  sub-grid  scale  mixing  is  crudely  estimated,  since  the  horizontal 
averaging  scale  is  much  larger  than  the  vertical  scale,  and  due  to  a  lack  of  understanding 
of  the  structure  of  horizontal  mixing  over  heterogenous  surfaces.  Therefore,  the  form  of 
the  horizontal  sub-grid  scale  correlation  terms  is  chosen  to  control  non-linear  aliasing, 
and  not  to  represent  the  actual  physical  processes  accurately. 

Avissar  et  al.  (1990)  states  that  because  of  added  numerical  diffusion  schemes, 
motions  on  scales  no  less  than  4Ax  are  fully  resolved  in  numerical  models.  However, 
more  recent  studies  based  on  spectral  analysis  of  model  simulated  fields  yield  slightly 
different  conclusions.  Bryan  et  al.  (2003)  conducted  a  series  of  large  eddy  simulations  of 
deep  moist  convection  using  different  6**'  order  explicit  diffusion  coefficients.  Analysis 
of  the  energy  spectra  revealed  a  scale  at  which  there  was  a  discrete  separation  between 
the  physical  and  computational  solutions.  They  conclude  that  only  information  at  scales 
greater  than  6Ax  in  the  energy  spectra  represent  a  physical  solution  and  should  be 
considered  useful  in  most  applications.  This  conclusion  is  supported  in  Harris  et  al. 
(2001),  who  compare  an  explicit  simulation  at  3  km  resolution  of  a  squall  line  to  radar 
data  of  the  same  case.  Power  spectrum  analysis  gave  excellent  comparisons  from  large 
scales  to  about  5Ax,  supporting  their  conclusion  that  only  information  larger  than  this 
scale  can  be  considered  physically  real.  Both  Bryan  et  al.  (2003)  and  Harris  et  al.  (2001) 
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Figure  4.12,  COAMPS  surface  layer  vertieal  veloeity  speetra  from  simulations  of 
JES0216  using  4*''  order  diffusion  eoeffieients  of -0.002  (solid)  and  -.004  (dashed). 
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conclude  that  numerical  diffusion,  both  implicit  and  explicit,  is  the  responsible  for  the  fall 
off  in  variability  at  small  scales. 

Evidently,  the  numerieal  diffusion  seheme  will  affeet  the  physieal  representation 
of  perturbations  on  the  smallest  resolvable  seales.  To  test  this  diffusion  and  resolvability 
issue  for  the  JES0216  ease,  we  performed  two  COAMPS™  simulations  for  ease  JES0216 
with  different  4**^  order  diffusion  eoeffieients  (0.002  and  0.004).  A  speetral  analysis  was 
then  done  for  the  simulated  vertieal  velocity  field.  The  energy  density  speetra  along 
observation  leg  R02  (see  Eigure  4.7)  for  the  two  simulations  are  shown  in  Eigure  4.12. 
The  speetra  show  the  separation  between  physical  and  computational  solutions  in  the 
vertieal  veloeity  speetra  referred  to  in  Bryan  et  al.  (2003).  We  find  that  the  speetra  for 
both  runs  are  eomparable  to  eaeh  other  at  seales  between  3  km  and  15  km.  Below  3  km 
(6Ax),  the  magnitude  of  the  speetra  deereases  sharply  from  the  observations.  We 
eonelude  that,  for  the  JES0216  model  simulation,  we  eannot  expeet  to  fully  resolve 
perturbations  less  than  3  km  within  the  high-resolution  (Ax  =  500  m)  inner  grid.  Thus, 
while  the  model  can  theoretically  resolve  a  1  km  (2Ax)  wave  in  the  inner  grid,  all 
perturbations  at  or  below  the  3  km  seale  will  be  damped  to  some  degree,  if  not 
eompletely,  by  the  diffusion  seheme  and  cannot  be  well  represented  in  our  results.  In  the 
COAMPS™  simulated  results  presented  in  the  remainder  of  this  ehapter,  a  4**^  order 
diffusion  eoeffieient  of  0.002  was  used.  This  simulation  will  be  referred  to  as  the  eontrol 
run. 

Another  numerieal  issue  worthy  of  discussion  is  the  adveetion  effect  from  the 
outer  to  the  inner  domains.  Sinee  the  mean  wind  direction  is  from  the  northwest  in 
JES0216,  the  effeet  of  the  model  physies,  within  the  innermost,  highest  resolution 
domain  (i.e.,  smaller  explieit  and  numerieal  diffusion),  are  not  eompletely  seen  in  the 
model  output  in  the  northwestern  region  of  this  small  domain  (10,000  km  in  size), 
espeeially  near  the  upwind  boundaries.  As  a  result,  the  model  fields  for  the  northwestern 
region  of  the  domain  show  spatial  variability  similar  to  the  eorresponding  geographie 
region  in  domain  3.  As  features  advect  into  the  relatively  small  inner  domain  from  its 
parent  domain,  the  model  physies  at  higher  resolution  will  have  more  notieeable  effeet  as 
the  down-wind  distanee  inereases  and  the  spatial  variation  of  the  fields  gets  eloser  to  the 
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Figure  4,13.  Energy  density  speetra  COAMPS™  vertieal  veloeity  at  40  m  above  sea 
level  for  the  parent  domain  (solid)  and  the  innermost  domain  (dashed)  from  the 
simulation  of  JES0216  in  both  (a)  northern  (along  flight  leg  R06)  and  (b)  southern 
(along  flight  leg  R02)  geographie  regions  within  the  innermost  domain.  See  Eigure  4.7 
for  the  loeations  of  the  flight  legs. 
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true  resolvability  of  the  higher  grid  resolution.  To  illustrate  this  eoncept  for  the  mean 
fields,  a  eomparison  between  the  innermost  domain  and  its  parent  domain  (domain  3  in 
Table  4.1)  of  the  spectra  in  the  BL  (at  400  m  above  sea  level)  for  the  w  wind  component 
is  shown  in  Figure  4.13.  The  spectra  for  each  of  the  domains  are  calculated  from  model 
output  along  west  to  east  lines  for  R02  (southernmost  flight  leg)  and  R06  (northernmost 
flight  leg,  see  Figure  4.7)  It  is  clear  from  the  figure  that  the  spectral  energy  in  the 
resolvable  scales  (scale  (wave  number)  greater  than  3  km  (2  km'^)  for  the  inner-most 
domain  and  10  km  (0.6  km"')  for  the  parent  domain)  is  similar  for  both  domains  in  the 
northern  region  of  domain  4,  and  increases  with  distance  north  to  south  for  each  domain. 
However,  the  increase  in  spectral  energy  for  the  innermost  domain  from  north  to  south, 
nearly  one  order  of  magnitude,  is  clearly  larger  than  for  the  parent  domain.  This 
illustration  of  the  down-wind  effect  of  the  higher  resolution  model  physics  on  the  spatial 
variability  within  a  given  domain  will  be  important  in  chapter  V,  when  comparisons  are 
made  between  two  different  simulations  with  different  SGS  mixing  lengths. 

2,  Analysis  of  COAMPS™  Results  from  the  Control  Run 

a.  BL  Development 

We  will  now  determine  the  extent  to  which  COAMPS™  represents  the  BL 
structure  and  features  typical  of  CAO  conditions.  Figure  4.14  shows  the  COAMPS™ 
near-surface  wind  and  temperature  fields  within  the  13.5  km  resolution  outer  domain, 
depicting  conditions  that  exist  within  the  innermost  domain.  The  surface  wind  and 
surface  air  temperatures  over  the  geographic  region  of  the  innermost  domain  (outlined  in 
Figure  4.14)  are  from  the  northwest  with  magnitudes  up  to  30  kts  (15  m  s"')  and  between 
261  and  263  K,  respectively.  The  temperature  gradient  over  most  of  the  JES  is  relatively 
perpendicular  to  the  wind  direction.  The  effects  of  topography  on  the  wind  magnitude 
are  evident,  where  a  region  of  weaker  wind  within  the  innermost  domain  is  present 
downstream  of  the  coastal  mountains. 

Horizontal  and  vertical  cross-sections  of  wind  speed,  vertical  velocity, 
potential  temperature,  and  moisture  are  shown  in  Figures  4.15  and  4.16.  The  horizontal 
cross-sections  for  the  entire  innermost  domain  region  are  shown  for  a  height  of  390  m 
(near  the  middle  of  the  BL).  The  vertical  cross  section  is  taken  along  the  line  shown  in 
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Figure  4,14,  COAMPS  12.5  km  resolution  domain  forecast  fields,  Valid  17  Feb  2000, 
OOOOZ.  The  figure  shows  the  wind  speed  in  kts  (barbs)  and  potential  temperature  in  K 
at  40  m  above  sea  level. 
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Figure  4.15.  COAMPS™  JES0216  simulated  fields,  Valid  17  Feb  2000,  OOOOZ. 
Horizontal  cross-sections  at  390  m  above  sea  level  and  vertical  cross-sections  along  the 
line  shown  in  the  horizontal  cross-sections  for  potential  temperature  (a  and  b)  and  water 
vapor  (c  and  d).  Contour  intervals  are  0.2  K  for  potential  temperature  and  0.05  g  kg'^ 
for  water  vapor. 
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the  horizontal  cross-sections,  in  the  down-wind  direction  (viewing  the  plot  from  left  to 
right),  and  from  the  surface  to  above  the  BL.  To  avoid  erroneous  data  near  the  domain 
boundaries,  the  plot  begins  and  ends  10  grid  points  (5  km)  in  from  the  boundaries.  Figure 
4.15a  show  small  spatial  variation  in  potential  temperature  that  resembles  that  of  the  SST 
variation  in  the  region,  indicating  the  strong  surface  forcing  in  this  case.  The  winter 
boundary  layer  over  the  simulated  region  is  generally  dry,  with  the  specific  humidity 
varying  around  1  g  kg‘\  However,  downstream  moistening  is  evident  from  both  the 
horizontal  cross-section  (Figure  4.15c)  and  the  downwind  vertical  cross-section  (Figure 
4.15d).  From  the  vertical  cross-section  of  potential  temperature  (Figure  4.15b),  we  can 
also  identify  the  downwind  development  of  the  boundary  layer  as  the  air  column  moved 
over  the  JES.  Here,  the  boundary  layer  growth  was  significant  in  the  first  40  km  of  the 
vertical  cross-section.  There  was  also  rapid  warming  as  the  BL  grew.  The  BL  height 
stabilized  at  around  700  m  above  sea  level  beyond  the  first  40  km  where  the  BL  became 
better  mixed  vertically  in  both  potential  temperature  and  water  vapor.  This  simulated  BL 
height  is  lower  than  the  observation  from  the  aircraft  soundings  (above  900  m),  likely  the 
result  of  an  under-predicted  BL  air  temperature.  Compared  to  the  aircraft  observations, 
the  BL  potential  temperature  was  in  general  under-predicted  by  COAMPS^"^  by  1  to  1 .5 
K.  COAMPS™  generally  has  a  cold  bias  recognized  by  NRL  scientists  (Dr.  James 
Doyle,  NRL  Monterey,  personal  communication).  Unfortunately,  the  source  of  this  cold 
bias  is  not  known  yet.  The  simulated  BL  water  vapor  is  consistent  with  the  observed 
amount.  We  also  note  that  while  some  thin  clouds  are  seen  in  the  satellite  imagery  of 
Ligure  4.5  for  this  case,  COAMPS™  did  not  produce  any  clouds  in  the  control  run. 

The  simulated  wind  speed  in  the  down-wind  direction  at  390  m  above  sea 
level  (Ligure  4.16  a)  shows  a  region  of  weaker  wind  with  little  spatial  variability  in  the 
center  of  the  grid,  along  the  line  of  the  vertical  cross-section,  and  stronger  winds  with 
slightly  more  spatial  variation  in  other  regions.  As  shown  in  the  outer  domain  wind  field 
(Ligure  4.14),  the  difference  in  wind  speed  magnitudes  across  the  domain  are  likely  due 
to  the  topographic  effects  of  the  nearby  coastal  mountains  and  valleys  to  the  west.  The 
simulated  mean  BL  wind  speed  magnitude  along  the  center/down-wind  region  of  the 
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Figure  4,16,  Same  as  in  Figure  4.15,  exeept  for  wind  speed  (a  and  b)  and  vertieal 
veloeity  (e  and  d).  Contour  intervals  are  1  ms'^  for  wind  speed  and  0.1  ms"'  for  vertieal 
veloeity. 
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domain  is  around  14  m  s'^,  comparable  to  the  near-surfaee  aireraft  observations  (see 
Figure  4.9).  However,  the  inerease  in  wind  speed  magnitudes  in  the  southwest  and 
northeast  eomers  of  the  domain  is  not  as  large  in  the  observations  and  thus  are  over- 
predieted  by  the  model.  In  addition,  more  horizontal  spatial  variability  in  the  simulated 
wind  speed  is  expeeted  in  the  resolvable  seales,  sinee  the  wind  speed  variations  in  the 
aireraft  observations  in  Figure  4.9  are  on  the  seale  of  4  to  10  km.  In  the  vertieal  eross- 
seetion  (Figure  4.16  b)  there  is  indeed  little  variation  along  the  down-wind  direetion 
through  the  depth  of  the  BL,  but  signifieant  wind  shear  is  seen  above  the  inversion,  with 
an  evident  infleetion  point  at  around  1  km  above  sea  level.  Unfortunately,  the  wind 
speed  inerease  above  the  BL  eannot  be  verified  sinee  there  were  no  observations  above 
the  BL. 

The  simulated  vertieal  veloeity  eross-seetions  (Figure  4.16  e  and  d)  should 
indieate  the  BL  updraft  and  downdraft  regions  on  the  resolvable  seale.  While  the 
observations  showed  near-surfaee  updrafts  and  downdrafts  of  up  to  1  m  s’'  in  magnitude, 
the  model  prediets  magnitudes  along  the  eross-seetion  within  the  surfaee  layer,  and  most 
of  the  BL,  of  only  0.2  m  s'*.  One  noteworthy  roll-like  turbulent  eddy  strueture  is  seen  on 
the  horizontal  eross-seetion  at  390  m  above  sea  level  in  the  southwest  eorner  of  the  grid, 
with  updraft/downdraft  magnitudes  eomparable  to  those  expeeted  near  the  middle  of  a 
eonveetive  BL,  and  a  width  of  around  5  km.  Another  more  eellular-like  eddy  appears  in 
the  NE  eorner  of  the  horizontal  eross-seetion.  Other  than  these  isolated  features,  the 
simulated  vertieal  veloeity  magnitudes  overall  are  too  weak,  and  we  see  again  that  the 
spatial  variability  of  eddies  is  in  general  mueh  less  than  in  the  observations. 

Another  notable  feature  in  the  vertieal  eross-seetion  of  vertieal  veloeity  is 
the  presenee  of  eddies  above  the  BL  inversion.  These  features  are  likely  internal  gravity 
waves,  whieh,  as  diseussed  in  Clark  et  al.  (1986),  Mason  and  Sykes  (1982)  and  others, 
are  features  that  exist  in  nature  as  well  as  being  simulated  in  numerieal  models.  These 
features  are  eommonly  seen  in  eases  with  a  eonveetive  BL  and  strong  vertieal  wind  shear 
above  the  inversion.  The  waves  are  produeed  as  thermally-foreed  eddies  formed  in  the 
BL  push  up  on  the  inversion  and  into  the  stable  layer.  These  ripples  then  produee 
vertieally-propagating  gravity  waves.  The  wind  shear  above  the  inversion  also 
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contributes  to  the  gravity  wave  produetion  and  propagation,  as  the  differential  air  motions 
exist  relative  to  the  thermals.  The  shear  is  obstructed  (known  as  “form  drag”)  by  the  BL 
eddies,  and  as  the  wind  flows  up  and  over  eddies  to  uphold  the  continuity  of  mass, 
vertically  propagating  waves  form.  Thermal  foreing  is  a  pre-requisite  for  the  produetion 
of  internal  gravity  waves,  while  the  shear  faeilitates  the  wave  produetion  and 
propagation.  Mason  and  Sykes  (1982)  obtained  numerical  solutions  without  clouds  where 
BL  eddies  were  produced  and  coupled  with  internal  gravity  waves  of  distinet  but  small 
magnitude  (0.4  m  s"').  Clark  et  al.  (1986)  simulated  idealized  cases  of  conveetive  BLs 
where  internal  gravity  waves  were  produced  in  the  eases  of  moderate  to  strong  shear,  and 
describe  the  development  of  the  gravity  waves.  Once  the  eddies  produeed  in  the 
simulation  fill  the  BL,  they  become  obstaeles  to  the  flow,  causing  the  formation  of  the 
gravity  waves,  whieh  are  launehed  vertieally  up  and  over  the  BL  thermals  into  the  stable 
layer.  The  eddies  propagate  both  horizontally  and  vertically  away  from  their  origin, 
inereasing  in  both  horizontal  and  vertical  scale.  While  this  is  all  very  interesting,  we  will 
eonelude  by  aeknowledging  the  realistic  existence  of  these  features  as  eited  and  diseussed 
in  the  literature  and  move  on.  Again,  the  existence  of  the  simulated  internal  gravity 
waves  is  however  not  verified  for  our  case,  since  we  have  no  observations  above  the  BL. 

b.  Surface  Fluxes 

While  we  will  save  more  detailed  deseriptions  of  the  parameterized 
turbulent  fields  (fluxes  and  TKE)  for  next  ehapter,  we  will  briefly  diseuss  the  model 
surface  fluxes  (Figure  4.17)  in  comparison  to  observation.  The  heat  and  momentum 
fluxes  show  a  region  of  relatively  weaker  fluxes  in  the  eenter  of  the  domain,  along  the 
diagonal  in  the  down-wind  direetion,  which  corresponds  to  the  region  of  relatively 
weaker  (but  eomparable  to  observations)  and  relatively  homogenous  simulated  horizontal 
winds.  The  combined  sensible  and  latent  heat  fluxes  are  eomparable  to  observations  in 
the  eenter/down-wind  region,  and  clearly  overestimated  in  the  regions  of  over-predicted 
winds  in  the  northeast  and  southwest.  Also,  beeause  the  comparison  here  is  between  bulk 
flux  method  parameterized  surface  fluxes  and  the  observed  fluxes  at  40  m  above  sea 
level,  we  do  expeet  some  differences  between  the  observed  and  simulated  surfaee  fluxes. 
There  is  better  horizontal  spatial  variability  in  the  fluxes  in  the  southeastern  region  of  the 
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Figure  4,17,  Horizontal  cross  sections  from  JES0216  control  run,  valid  17  Feb  2000, 

OOOOZ  for  (a)  surface  sensible  and  (b)  surfaee  latent  heat  fluxes  (W  m'  ),  and  (c)  surface 
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wind  stress  (N  m'  ).  Contour  intervals  are  5  W  m'  for  sensible  and  latent  and  0.02  N 
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m‘  for  wind  stress. 
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domain  (outlined  in  each  of  cross-sections  in  Figure  4.14),  which  was  not  as  evident  in 
the  cross-sections  of  the  simulated  mean  fields.  This  is  due  to  the  affect  of  advection 
from  the  1.5  km  resolution  parent  domain  into  the  500  m  resolution  innermost  domain, 
with  the  model  physics  improving  the  spatial  resolvability  down-wind,  as  discussed 
earlier.  More  issues  will  be  raised  and  closer  model-observations  comparisons  will  be 
made  using  the  SGS  BL  turbulent  fluxes  in  the  next  chapter. 

We  also  calculate  the  zi/L  from  the  simulated  surface  sensible  heat  flux, 
u*,  BL  height  and  potential  temperature.  The  mean  zi/L  for  the  domain  4  region  is 
calculated  as  -5,  which  is  comparable  to  the  calculation  made  from  the  observations  and 
thus  also  predicts  a  regime  where  convective  rolls  would  develop. 

F,  SUMMARY  AND  CONCLUSIONS 

In  this  chapter,  we  described  the  JES0216  case  study.  Global  model  analysis 
fields  showed  a  persistent  synoptic  situation  in  the  region  of  interest  in  the  western  JES 
(post-cold  frontal,  CAO  conditions  with  strong  subsidence)  between  16  Eebruary  1200  Z 
and  17  Eebruary  0000  Z.  Satellite  imagery  and  aircraft  observations  illustrated  the 
presence  of  convective  roll  vortices  in  the  location  of  the  high-resolution  inner  domain  of 
our  control  run  simulation.  The  value  of  Zi/E  calculated  from  the  aircraft  observations 
was  -  4.6,  which  falls  within  the  Zi/E  range  given  in  the  literature  for  convective  roll 
development.  The  aircraft  observations  also  indicated  large  air-sea  temperature 
differences  coupled  with  relatively  high  near-surface  winds.  The  combined  effects 
generated  large  surface  heat  fluxes,  which  drove  the  development  of  the  rolls.  The  scale 
of  the  rolls,  based  on  mean  field  perturbations  in  the  near-surface  aircraft  observations 
was  around  3  to  5  km. 

A  COAMPS™  control  run  simulation  for  the  JES0216  case  was  performed.  In- 
depth  analysis  of  the  model  output  fields  showed  that  that  the  model  is  adequately 
representing  the  mean  state  of  the  BE  for  this  post-frontal,  CAO  case.  The  control  run 
showed  good  comparison  to  observations  in  terms  of  magnitudes  of  the  mean  fields  and 
surface  heat  fluxes.  The  control  run  realistically  simulated  development  of  the  BE  in  the 
down-wind  direction,  becoming  warmer  and  moistening  with  distance  off-shore. 
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However,  the  control  run  failed  to  simulate  the  convective  roll  vortices.  Simulated 
vertical  velocities  were  too  weak,  and  there  was  an  overall  lack  of  horizontal  spatial 
variability  in  the  simulated  mean  fields  compared  to  the  observations.  In  light  of  the 
discussion  on  the  explicit  and  numerical  effects  of  diffusion  and  advection  on  the 
simulation,  the  results  indicate  that  the  BL  structure  on  the  true  resolvable  scale  (>  3  km) 
is  not  well  represented  by  the  model.  Therefore,  adjustments  will  be  made  within 
COAMPS™  in  an  attempt  to  improve  the  simulation  of  the  resolvable  scale  features. 
The  method,  objective  and  rationale  for  this  adjustment,  and  the  results  of  the  simulation 
using  the  modified  COAMPS™  code,  will  be  discussed  in  chapter  V. 
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V.  COAMPS™  CASE  STUDY  SIMULATION  USING  A 

RESOLUTION-DEPENDENT  TURBULENCE 
PARAMETERIZATION 

A,  MODIFICATION  TO  COAMPS™  MIXING  LENGTH 

1.  Rationale  and  Motivation 

It  is  apparent  from  the  results  in  chapter  IV  that,  based  on  the  known  limits  of 
resolvability  in  the  COAMPS™  model,  the  highest  resolution  domain  should  be 
resolving  a  greater  degree  of  variability.  Our  proposed  theory  to  explain  the  lack  of 
variation  in  the  model  results  for  the  JES0216  control  run  is  that  the  model  mixing  length, 
used  in  the  turbulent  parameterization  calculations,  is  too  large.  Based  on  the  use  of 
mixing  length  in  the  SGS  parameterization  of  the  vertical  fluxes,  with  the  turbulent 
mixing  length  too  large,  the  vertical  mixing  that  occurs  with  each  time  step  in  the  model 
would  generate  too  much  mixing  of  properties  in  the  BL,  and  thus  likely  underestimates 
the  spatial  variability  on  the  true  resolvable  scale  (6Ax). 

Recall  from  the  length  scale  analysis  in  chapter  III  that,  for  a  horizontal  grid 
spacing  that  falls  below  the  scale  of  large  turbulent  eddies  (as  is  in  the  JES0216 
simulation  with  a  Ax  of  500  m)  the  integral  length  scale  profile  should  be  smaller  in 
magnitude  than  that  calculated  by  an  ensemble  turbulence  parameterization  scheme  (see 
Eigure  3.9).  The  smaller  length  scale  would  more  correctly  represent  the  turbulent  eddies 
that  scale  at  and  below  the  technical  resolvability  (2Ax)  of  the  domain.  It  can  be 
hypothesized,  based  on  the  reference  by  Chou  and  Zimmerman  (1989)  to  the  dependence 
of  the  mesoscale  roll  vortices  on  the  large  scale  turbulent  eddies,  that  improved 
simulation  of  the  large-scale  turbulence  could  also  improve  the  simulation  of  the 
mesoscale  convective  rolls. 

This  chapter  will  discuss  the  changes  made  to  the  mixing  length  calculation  in 
COAMPS™  in  an  effort  to  improve  the  horizontal  spatial  variability  of  the  high- 
resolution  output  to  be  more  realistic,  and  compare  the  modified  mixing  length  simulated 
output  to  the  JES0216  control  run  simulation  and  to  observations  to  see  if  improvement  is 
indeed  made.  Based  on  the  results  of  our  control  simulation  of  the  JES0216  case,  we  will 
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apply  an  attempted  improvement  to  the  COAMPS™  SGS  parameterization  in  the  high- 
resolution  innermost  domain,  in  an  effort  to  better  simulate  the  true  spatial  variability  of 
the  atmosphere  in  high-resolution  output. 

2.  Modified  Mixing  Length  Formulation  for  COAMPS™ 

Based  on  our  study  in  previous  chapters,  we  introduce  a  resolution-dependent 
mixing  length  (Iax)  formulation  into  COAMPS^'^  with  modifications  suitable  for  the 
model.  This  new  formulation  (applied  only  within  the  high  resolution  nest  of  our 
simulation)  adjusts  the  original  COAMPS^*^  ensemble  turbulence  parameterization 
mixing  length  calculated  using  the  MY82  vertically  integrated  TKE  formulation,  to  a 
magnitude  and  profile  appropriate  for  the  resolution  of  the  simulation.  The  equation  used 
to  adjust  the  mixing  length  is  based  on  equation  3.3,  but  adjusted  for  COAMPS™.  Here, 
the  integral  scale  accounting  for  all  scales  of  turbulence  (Aens  in  equation  3.3)  is  replaced 
with  the  mixing  length  as  originally  calculated  in  COAMPS™  using  the  MY82 
formulation  (1my82),  and  the  integral  scale  accounting  for  the  smallest  scales  considered 
(A2oom  in  equation  3.3)  is  taken  as  a  constant  value,  Iles,  of  37.5  m.  Therefore,  equation 
3.3  becomes,  for  our  numerical  adjustments. 
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with  R  equal  to  3.8  when  Ax  is  in  km,  as  in  equation  3.3.  As  RAx  becomes  large, 
equation  5.1  approaches  1my82,  where  all  turbulence  scales  below  2Ax  and  should  be 
entirely  accounted  for  in  the  parameterization.  As  Ax  approaches  LES  resolution  (i.e..  Ax 
~  0(100  m)),  equation  5.1  approaches  Iles,  which  would  be  appropriate  for 
parameterization  of  eddies  in  the  inertial  sub-range  and  below.  Thus,  this  modified 
equation  allows  for  a  smooth  transition  from  a  mixing  length  appropriate  for  total 
ensemble  parameterization  at  larger  Ax,  to  one  appropriate  for  an  EES-type 
parameterization  at  small  Ax. 

Equation  5.1  was  applied  in  a  simulation  of  JES0216,  where  the  simulation  is 
identical  to  that  presented  in  chapter  IV,  except  for  the  modification  to  the  mixing  length. 
We  will  refer  to  the  modified  simulation  as  RESDEP,  for  the  resolution-dependent 
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mixing  length  modification.  The  modified  mixing  length  at  the  12-hour  forecast  time  is 
shown  in  Figures  5.1  and  5.2.  Here,  Figure  5.1  shows  the  domain- averaged  mixing 
length  profile.  The  magnitude  of  the  resolution-dependent  mixing  length  is  seen  to 
decrease  in  magnitude  from  that  of  the  control  run  by  about  one-half.  Based  on  the 
results  of  the  observational  study  of  integral  scales  in  chapter  111,  the  magnitude  decrease 
and  shape  of  the  profile  is  appropriate  for  parameterization  of  1  km  eddies  and  smaller. 
The  horizontal  cross-section  of  mixing  length  at  390  m  above  sea  level  is  shown  in 
Figure  5.2.  This  level  is  near  the  middle  of  the  BL  where  the  magnitude  of  mixing  length 
is  generally  largest,  as  seen  from  Figure  5.1.  From  Figure  5.2,  we  see  the  general 
decrease  in  the  magnitude  shown  in  Figure  5.1.  The  steeper  slope  of  the  NOMOD 
mixing  length  profile  (Imysi)  near  390  m  generates  more  variation  of  the  mixing  length  in 
the  horizontal  cross-section  across  the  domain  at  this  level.  Because  the  mixing  length 
formulations  are  sensitive  to  the  vertically  integrated  TKE  (Equation  2.21),  the  horizontal 
distribution  of  the  mixing  length  is  also  dependent  on  the  distribution  of  TKE  in  the  BE, 
which  is  driven  by  the  surface  forcing  (total  heat  flux  for  JES0216). 

Next  we  will  examine  the  extent  to  which  this  new  resolution-dependent  mixing 
length  formulation  affects  the  simulation  of  the  BL  structure  for  JES0216.  Recall  from 
chapter  IV  that  the  Zi/L  of  -5,  calculated  from  both  observations  and  the  NOMOD 
simulation,  suggested  the  possible  formation  of  horizontal  roll  vortices  within  the  BL  for 
JES0216.  Indeed  the  roll  vortices  were  observed  at  5  km  spacing.  Erom  our  test 
simulations  on  model  numerics  effects,  COAMPS™,  at  500  m  resolution,  should  be  able 
to  resolve  these  roll  structures.  Unfortunately,  the  rolls  were  not  simulated  well  in  the 
NOMOD  simulation.  It  is  hoped  that  the  modification  to  the  mixing  length  in  the 
RESDEP  simulation  will  improve  the  simulated  structure  of  the  convective  rolls.  We 
will  analyze  and  compare  the  control  and  modified  simulations  from  two  aspects:  the 
differences  in  the  resolvable  scale  structure  of  the  BL,  in  terms  of  both  magnitude  and 
spatial  variability  of  the  mean  quantities,  and  the  differences  in  the  magnitude  and 
structure  of  the  SGS  parameterized  fluxes  and  TKE. 
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Figure  5,1,  COAMPS™  domain-averaged  profiles  of  mixing  length  (m)  from 
NOMOD  (solid)  and  from  the  resolution-dependent  formulation  from  RESDEP 
(dashed)  for  JES0216.  Data  shown  is  from  the  12-hour  forecast  fields,  valid  17  Eeb 
2000,  OOOOZ. 
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Figure  5,2,  Same  as  in  Figure  5.1,  exeept  horizontal  eross-sections  for  the  entire 
domain  at  390  m  above  sea  level,  for  NOMOD  (a)  and  RESDEP  (b)  simulations. 
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B,  EFFECTS  OF  MODIFIED  MIXING  LENGTH  ON  SIMULATING  THE 

ROLL  STRUCTURE 

1.  Boundary  Layer  Mean  Structure 

In  the  following  comparisons  of  the  mean  quantities  (horizontal  wind, 
temperature,  and  water  vapor),  we  will  focus  on  the  southeastern  region  of  the  domain, 
where  the  spatial  variability  is  considered  to  be  representative  of  the  500  m  resolution 
(due  to  advection  affect  from  the  parent  nest,  as  discussed  in  chapter  IV).  This  region  is 
defined,  from  inspection  of  the  spatial  variability  of  the  SGS  fluxes  from  the  control  run, 
by  the  region  outlined  within  the  high-resolution  domain  in  Figure  4.17.  Unless 
otherwise  noted,  the  horizontal  boundaries  of  the  horizontal  and  vertical  cross-sections 
presented  in  this  chapter  are  defined  by  this  region  of  the  domain.  In  addition,  the 
vertical  cross-sections  are  now  taken  along  a  West-to-East  line  at  the  southern  end  of  this 
region  instead  of  a  line  in  the  downwind  direction  as  was  done  for  the  vertical  cross- 
sections  presented  in  chapter  IV. 

In  the  comparison  between  the  horizontal  and  vertical  cross-sections  of  the  mean 
quantities  between  NOMOD  and  RESDEP  in  Eigures  5.3  to  5.6,  it  is  immediately 
apparent  that  the  horizontal  spatial  variability  is  increased  in  RESDEP.  In  the  horizontal 
cross-sections,  a  clear  signature  of  the  observed  roll  structure  is  now  simulated  in  all  of 
the  mean  quantities,  as  represented  by  the  longitudinal  alignment  of  small-scale  features 
in  the  horizontal  wind,  temperature,  and  moisture  with  the  mean  wind  direction  (from 
NW  to  SE),  which  were  not  present  in  the  NOMOD  simulation.  The  signature  of  the 
rolls  is  clearest  in  the  water  vapor  and  vertical  velocity  fields.  Erom  the  horizontal  cross- 
sections  of  water  vapor  and  vertical  velocity  from  RESDEP  (Eigures  5.5  c  and  5.6  c),  the 
scale  of  the  rolls  (the  measured  distance  across  adjacent  moist  updraft  and  drier 
downdraft  features)  vary  between  3  and  7  km,  which  is  fairly  consistent  with  the  scale  of 
the  rolls  determined  from  the  aircraft  observations.  The  BE  height  and  structure  of  the 
inversion  in  RESDEP  is  very  similar  to  that  of  NOMOD,  clearly  seen  in  the  vertical 
cross-sections  of  potential  temperature  (Eigure  5.4  b  and  d).  Using  a 
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Figure  5.3,  Comparisons  of  simulated  wind  fields  between  NOMOD  and  RESDEP. 
Horizontal  cross-sections  and  W-E  vertical  cross-sections  of  the  horizontal  wind  (m  s"') 
for  NOMOD  (a  and  b)  and  RESDEP  (c  and  d).  Area  shown  in  horizontal  cross-sections 
(a  and  c)  is  the  southeastern  region  (-60x60  km  )  of  the  domain.  Vertical  cross-sections 
(b  and  d)  are  along  the  line  shown  in  the  horizontal  cross-sections.  Contour  interval  is  1 
m  s"'  for  all  plots. 
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Figure  5,4,  Same  as  in  Figure  5.3,  except  for  potential  temperature  (K).  Contour 
interval  is  0.25  for  all  plots. 
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Figure  5,5,  Same  as  Figure  5.3,  except  for  water  vapor  (g  kg'').  Contour  interval  is  0.1 
for  all  plots. 
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Figure  5,6,  Same  as  in  Figure  5.3,  except  for  vertical  velocity  (m  s"').  Contour  interval 
is  0.5  for  all  plots. 


no 


mean  BL  height  of  this  eross-seetion  (around  600  m),  the  aspeet  ratio  of  the  rolls  is  5  to 

1 1 .6,  whieh  is  somewhat  larger  (due  to  the  lower  than  observed  simulated  BL  depth)  than 
the  4  to  6  range  of  aspeet  ratio  caleulated  from  the  observations  in  chapter  IV,  but  is 
within  the  range  for  a  BL  containing  roll  vortices  (Etling  and  Brown  1993).  In  Figure 

5.6,  the  differences  in  eddy  structure  modeled  by  the  two  simulations  are  evident.  In  the 
NOMOD  run,  the  only  roll  updraft  and  downdraft  events  were  simulated  on  the  western 
end  of  the  cross-section,  with  weak  resolvable  scale  vertical  motion  for  the  rest  of  the 
cross-section.  In  contrast,  the  RESDEP  simulation  shows  considerable 
updrafts/downdrafts  throughout  the  cross-section,  which  clearly  defines  the  role  structure 
of  the  BE.  Furthermore,  the  RESDEP  results  show  narrower  updraft  cores  and  broader 
downdraft  regions  being  realistically  simulated.  The  vertical  cross-sections  for  the 
different  mean  quantities  are  cutting  through  the  rolls  at  about  a  45  degree  angle,  and 
together  they  show  the  simulated  BE  roll  structure  clearly,  with  the  warm,  moist  eddies 
corresponding  to  updraft  cores,  and  cooler,  drier  eddies  corresponding  to  the  downdraft 
cores.  The  simulated  magnitude  of  roll  updrafts/downdrafts  in  RESDEP  is  also  larger 
and  more  realistic,  with  the  maximum  updraft  core  reaching  1.3  m  s'\  twice  that  of  the 
NOMOD  simulation. 

To  further  determine  if  the  spatial  variability  and  magnitudes  of  the  mean  fields 
simulated  in  RESDEP  is  realistic,  we  directly  compare  the  simulated  fields  to  aircraft 
observations.  Due  to  the  relative  lack  of  spatial  coverage  of  the  observations  for 
JES0216,  our  model-observation  validation  is  limited  only  to  near  the  surface  and  to 
profiles  within  the  BE.  The  simulated  mean  quantities  are  compared  to  the  JES  sounding 
observations  in  Figures  5.7  and  5.8.  The  sounding  used  for  comparison  in  this  figure  is 
S7,  whose  location  within  the  domain  is  shown  in  Figure  4.7.  Because  the  aircraft  travels 
8  to  10  km  from  SE  to  NW  during  the  sounding,  there  is  significant  variability  seen  in  the 
observed  profiles  of  the  horizontal  wind  components  (Figure  5.7),  and  some,  but  less, 
variability  in  the  observed  water  vapor  and  potential  temperature  profiles  (Figure  5.8)  as 
the  aircraft  passes  through  locally  different  air  masses  in  the  horizontal.  This  variability  is 
not  present  in  the  simulated  profile,  which  is  made  up  of  data  from  the  corresponding 
vertical  levels  at  the  same  horizontal  grid  point,  located  in  the  region  of  the  S7 
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Figure  5,7,  COAMPS  vertical  profiles  of  (a)  wind  speed,  (b)  wind  direction,  (c)  U 
wind  component,  and  (d)  V  wind  component  (units  shown)  for  NOMOD  (solid), 
RESDEP  (dashed),  and  aircraft  sounding  observations  taken  near  the  location  of  the 
model  vertical  profiles  (dotted). 
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Figure  5,8,  Same  as  Figure  5.7,  exeept  for  potential  temperature  (a)  and  water  vapor 
(b). 
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observations.  We  observe  only  small  differences  between  the  simulated  profiles  at  this 
particular  location  in  the  domain,  and  in  comparison  to  observations,  there  is  fairly  good 
agreement  in  both  the  trend  with  height  and  magnitudes  for  all  quantities  except  the 
potential  temperature,  which  is  simulated  between  1  and  2  degrees  too  cool,  as  already 
discussed  in  chapter  IV.  This  is  at  least  part  of  the  reason  for  the  lower  than  observed  BL 
depth  and  smaller  than  observed  roll  aspect  ratio.  Significant  wind  shear  in  both 
magnitude  and  direction  are  seen  above  the  BL  inversion.  The  simulated  profiles  for 
RESDEP  clearly  show  significant  wind  shear  near  the  middle  of  the  BE,  with  an 
inflection  point  at  a  height  of  about  530  m,  and  significant  veering  of  the  wind  within  the 
BE  as  well.  This  BE  wind  shear  depicts  the  inflection  point  instability  that,  in  conjunction 
with  convective  instability,  enhances  the  development  of  convective  rolls,  as  discussed  in 
Etling  and  Brown  (1993).  The  signature  of  the  inflection  point  instability  is  not  present 
in  the  NOMOD  wind  profiles.  The  wind  shear  above  the  inversion  is  present  in  both 
simulations,  and  enhances  the  internal  gravity  waves  in  the  free  atmosphere  above  the  BE 
(discussed  in  chapter  IV),  visible  in  the  vertical  velocity  vertical  cross-sections  (Eigure 
5.6  b  and  d)  for  both  simulations.  The  observed  and  simulated  BE  wind  profiles,  which 
show  very  little  wind  shear  within  the  BE,  also  indicate  that  the  mechanism  for  the  roll 
development  is  mainly  convective  instability  vice  inflection  point  instability  (Etling  and 
Brown  1993).  Comparisons  between  BE  profiles  from  NOMOD,  RESDEP,  and  other 
aircraft  soundings  yield  similar  results  to  those  shown  in  Eigures  5.8  and  5.9. 

2.  Spectral  Analysis  of  Resolvable  Scale  Mean  Fields 

The  horizontal  variations  of  the  mean  quantities  along  the  near-surface  flight  track 
R02,  near  the  south  end  of  the  domain  (see  Eigure  4.7)  are  compared  in  Eigures  5.9a  to 
5.13a.  Both  NOMOD  and  RESDEP  compare  fairly  well  in  magnitude  with  the 
observations  for  the  wind  components  and  water  vapor,  but  the  potential  temperature  is, 
again,  simulated  about  2  degrees  too  low  in  comparison  to  the  near-surface  observations, 
consistent  with  the  profile  of  Eigure  5.8.  The  simulated  vertical  velocity  at  the  flight 
level  is  significantly  weak  in  magnitude  for  both  NOMOD  and  RESDEP,  with  more 
variation  in  RESDEP  that  is  closer  to  that  observed. 
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While  the  inerease  in  horizontal  variability  is  evident  from  the  RESDEP 
simulation  along  the  near-surface  flight  leg,  the  scales  at  which  the  variability  increases  is 
more  clearly  evident  from  the  analysis  of  the  energy  density  spectra,  shown  in  Figures 
5.9b  to  5.13b  for  each  of  the  mean  quantities.  The  spectra  can  only  illustrate  the  resolved 
scale  variations  above  the  resolvable  limit  of  the  model.  In  chapter  IV,  this  resolvable 
limit  was  determined  to  be  6 Ax,  or  3  km  for  our  500  m  resolution  domain.  The  spectra 
are  plotted  versus  wave  number,  so  the  3  km  limit  corresponds  to  a  wave  number  of 
around  2  km'\  At  the  largest  scales  (smallest  wave  numbers),  the  statistical  significance 
of  the  spectra  for  both  the  model  output  and  observations  is  limited  by  the  length  of  the 
data  for  spectral  analysis,  which  is  about  60  km.  Therefore,  we  focus  our  spectral 
comparison  between  the  model  and  observations  to  within  the  scales  of  3  to  12  km  (wave 
numbers  0.5  to  2  km‘'),  where  we  consider  the  spectra  to  present  the  most  accurate 
representation  of  the  perturbation  energy.  Within  this  range,  we  observe  that  the 
NOMOD  spectra  energy  for  all  three  wind  components,  heat,  and  moisture,  is  one  to 
three  orders  of  magnitude  smaller  than  that  for  RESDEP,  while  the  RESDEP  energy 
spectra  compare  to  the  observed  spectra  within  one  order  of  magnitude,  with  vertical 
velocity  being  the  most  underestimated  (but  most  improved  over  NOMOD).  All  of  the 
observed  spectral  quantities  show  a  peak  in  energy  at  around  4  to  6  km,  indicative  of  the 
presence  of  the  convectively-driven  roll  vortices  at  that  scale  range  (LeMone  1976).  In 
general,  the  horizontal  variations  on  the  scale  of  the  rolls  in  RESDEP  are  much  more 
comparable  with  those  from  observations,  although  the  comparisons  in  vertical  velocity 
and  water  vapor  are  less  ideal.  This  is  a  significant  improvement  from  the  control  run 
with  the  original  COAMPS™  length  scale  formulation. 

We  also  calculated  the  co-spectra  between  vertical  velocity  and  scalar  variables 
such  a  potential  temperature  and  water  vapor  to  examine  the  resolvable  scale  flux 
transport  (Figure  5.14).  The  aircraft  observation  shows  sizable  co-spectra  between  wave 
numbers  0.8  and  2  km'\  suggesting  flux  contribution  from  the  roll  scale  perturbation  (~3 
to  8  km  wavelength).  The  RESDEP  simulation  reveals  similar  flux  contribution  by  the 
rolls  with  some  discrepancy  in  sensible  heat  flux  (Figure  5.14a).  In  contrast,  there  is  very 
little  co-spectral  energy  at  all  resolvable  scales  in  NOMOD.  Overall,  the  resolvable 
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Figure  5.9.  Comparison  of  the  horizontal  variability  of  the  U  wind  eomponent  among 
aireraft  observations  (dotted),  NOMOD  (solid),  and  RESDEP  (dash).  All  data  were  taken 
from  the  same  horizontal  path  along  R02  at  40  m  above  sea  level,  (a)  Spatial  series  plot; 
(b)  Energy  density  speetra.  The  speetra  are  cutoff  at  2  km"'  in  wave  number  (k) 
corresponding  to  a  wavelength  of  ~3  km.  The  two  short  vertical  bars  denote  the  range 
over  which  the  spectral  analysis  are  statistically  significant  (see  discussions  in  text). 
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Figure  5,10.  Same  as  in  Figure  5.9,  exeept  for  V  wind  component. 
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Figure  5,11,  Same  as  in  Figure  5.9,  exeept  for  vertieal  veloeity  (W). 
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Figure  5,12,  Same  as  Figure  5.9,  except  for  potential  temperature  (0). 


119 


^  ^  q(8k8 


(k)*k  C  (k)*k 


JES0216 


10  '  10° 
k  (km"^) 


Figure  5,15.  Co-spectra  between  vertical  velocity  (w,  in  ms'^)  and  (a)  potential 
temperature  (0,  in  K);  and  (b)  water  vapor  (q  in  g  kg''). 
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scale  co-spectra,  particularly  at  the  small  scales  close  to  the  resolvable  limit  of  the  model 
and  near  the  scale  of  the  convective  roll  vortices,  appear  to  be  more  realistically 
simulated  in  RESDEP,  based  on  the  comparisons  to  the  available  observations. 

3.  SGS  TKE  and  Fluxes 

The  improved  resolvability  of  the  BL  roll  structure  as  the  result  of  changes  in 
RESDEP  mixing  length  profile  also  result  in  changes  in  the  SGS  TKE  and  flux  fields. 
The  horizontal  cross-sections  for  the  simulated  SGS  sensible  heat  and  latent  heat  fluxes 
in  Figures  5.15  and  5.16  a  and  c,  are  shown  at  5  m  above  sea  level,  where  the  total  (SGS 
and  resolved)  heat  fluxes  are  large  (since  the  forcing  is  from  the  surface),  and  the  SGS 
fluxes  make  up  the  majority  of  the  total  heat  flux  within  the  roll  circulation  (EeMone 
1976).  For  NOMOD  (Figures  5.15a  and  5.6a),  the  sensible  and  latent  heat  fluxes  are 
relatively  invariant  in  the  horizontal.  The  vertical  cross-sections  for  the  NOMOD  heat 
fluxes  (Figures  5.15b  and  5.16b)  have  a  fairly  uniform  vertical  profile  across  the  domain, 
and  illustrate  a  BE  in  the  CAA  regime,  undergoing  both  warming  and  moistening  with 
time,  since  the  heat  and  moisture  fluxes  are  generally  decreasing  with  height  in  the  BE. 
As  expected  for  a  well-mixed  BE,  the  fluxes  are  largest  near  the  surface,  due  to  the  strong 
transfer  of  heat  and  moisture  from  the  relatively  warm  ocean  surface  to  the  cooler,  drier 
air  above.  We  also  see  pockets  of  negative  sensible  heat  fluxes  and  positive  latent  heat 
fluxes  near  the  BE  inversion,  where  warmer  and  drier  air  from  the  free  atmosphere  is 
being  mixed  downward  into  the  cooler,  moister  BE  air  in  the  regions  of  downdraft  eddies 
within  the  BE.  The  NOMOD  SGS  TKE  distribution  in  the  BE  (Figure  5.17  a  and  b) 
shows  no  longitudinal  alignment  in  the  plane  view  of  the  middle  of  the  BE,  and  has  a 
more  uniform  profile  along  the  vertical  cross-section,  showing  no  clear  signal  of 
enhancement  with  roll  convection.  Thus,  we  see  no  sign  of  the  roll  structure  in  the 
NOMOD  SGS  fields,  reflective  of  their  absence  in  the  mean  fields. 

In  contrast  to  NOMOD,  the  horizontal  cross-sections  for  RESDEP  SGS  sensible 
and  latent  heat  fluxes  near  the  surface  (Figures  5.15c  and  5.16c)  show  significantly  more 
spatial  variation  and  slightly  smaller  magnitudes.  This  structure  appears  to  be  realistic  for 
the  fluxes,  based  on  the  description  of  EeMone  (1976)  of  high-frequency  (i.e.,  SGS) 
turbulence  and  fluxes  of  varying  intensity  superimposed  on  the  mean  roll  circulation. 
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Figure  5,15,  Same  as  in  Figure  5.3,  exeept  for  SGS  Sensible  Heat  Flux  (W  m'  ). 
Horizontal  cross-seetion  is  at  a  height  of  5  m.  Contour  interval  is  25  for  all  plots. 
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Figure  5,16,  Same  as  in  Figure  5.15,  except  for  SGS  Latent  Heat  Flux  (W  m'  ). 
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The  cross-sections  of  RESDEP  TKE  (Figure  5.17  c  and  d)  present  a  TKE 
structure  more  indicative  of  a  BE  containing  convective  roll  vortices.  As  was  simulated 
in  the  mean  fields,  the  BE  TKE  structure  in  the  horizontal  at  390  m  above  sea  level  shows 
longitudinal  alignment  in  the  direction  of  the  mean  BE  wind.  The  locations  of  high  TKE 
coincide  with  the  locations  of  vertical  velocity  updraft  cores  (Figure  5.6d)  at  this  height 
in  the  BE.  The  implication  of  the  horizontal  structure  is  seen  in  the  vertical  cross-section 
of  the  RESDEP  SGS  TKE  (Figure  5.17d)  showing  TKE  transport  upwards  in  the  updraft 
regions,  while  the  TKE  in  the  downdraft  regions  remains  concentrated  near  the  surface. 
This  same  structure  is  seen  in  the  vertical  cross-sections  of  the  RESDEP  SGS  latent  and 
sensible  heat  fluxes  (Figures  5.15d  and  5.16d). 

The  structure  shown  in  the  simulated  flux  and  TKE  cross-sections  is  supported  by 
past  observational  and  numerical  studies.  The  observational  study  of  EeMone  (1976), 
states  that  the  circulation  of  roll  vortices  in  the  convective  BE  influences  the  distribution 
of  smaller  scale  turbulence  (analogous  to  the  SGS  turbulence  simulated  here),  with 
turbulent  variances  and  fluxes  systematically  concentrated  in  regions  of  positive  roll 
vertical  velocity.  EeMone  (1976)  refers  to  the  systematic  variation  of  turbulence  across 
the  rolls  as  modulated.  The  turbulence  modulation  was  found  to  be  the  result  of  the 
redistribution  of  turbulence-producing  elements  by  the  mean  roll  circulation.  Near  the 
surface,  the  mean  wind  in  the  along-roll  direction  transports  turbulent  elements 
horizontally  into  regions  of  upward  vertical  motion,  where  the  elements  are  then 
transported  upward.  If  the  turbulent  elements  are  smaller  scale  eddies,  the  eddies  will 
grow  upward  at  an  accelerated  rate.  EeMone  (1976)  also  states  that  the  modulation  of  the 
small-scale  turbulence  by  the  quasi-two  dimensional  rolls  is  maintained  near  the  surface 
as  well  as  at  higher  levels  in  the  BE.  Near  the  surface,  the  maintenance  is  due  to 
horizontal  variation  in  the  shear  production  from  the  mean  flow  and  buoyancy 
production.  At  higher  BE  levels,  upward  vertical  transport  of  the  TKE  and  fluxes  within 
the  roll  updraft  cores  maintain  the  modulation.  Chlond  (1992)  performed  an  idealized 
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COAMPS  JES0216  SGS  TKE 


2  2 

Figure  5,17,  Same  as  in  Figure  5.15,  except  for  TKE  (m  s').  Horizontal  cross-section 
is  390  m  above  sea  level.  Contour  interval  is  0.1  for  all  plots. 
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three-dimensional  large-eddy  simulation  based  on  an  observed  ease  of  cloud  street 
development  during  a  cold  air  outbreak  over  the  Greenland  Sea.  Chlond’s  numerical 
results  of  the  TKE  field  at  2  hours  into  the  simulation  are  illustrated  in  Figure  5.18,  and 
show  a  BL  TKE  structure  similar  to  that  described  in  LeMone  (1976)  and  RESDEP. 
Similar  TKE  structure  is  also  discussed  in  the  two-dimensional  simulation  study  of  roll 
vortices  developed  in  CAA  conditions  in  Raasch  (1990). 

In  validating  the  magnitudes  of  the  COAMPS^'^  SGS  parameterized  fluxes  and 
TKE  near  the  surface  with  that  derived  from  aircraft  observations,  it  is  important  to  note 
that  the  observed  fluxes  and  TKE  derived  from  the  moving  flux  averaging  represent  the 
TOTAE  turbulent  energy  from  all  scales,  while  the  model  parameterized  quantities  only 
represent  the  SGS  portion  of  the  energy.  In  order  to  compare  the  observations  with  the 
SGS  model  output,  it  is  necessary  to  calculate  only  the  observed  turbulent  energy  from 
scales  at  and  below  the  resolvable  limit  of  2Ax,  since  this  is  the  energy  that  should  be 
represented  by  the  model  SGS  parameterization  scheme.  The  observed  fluxes  at  and 
below  the  2Ax  scale  were  calculated  by  first  filtering  out  variances  due  to  scales  above  1 
km,  using  a  high-pass  filter  as  was  done  in  chapter  III.  The  filtered  observations  were 
used  to  calculate  a  moving  average  variance  and  co-variance,  where  the  length  of  the 
averaging  leg  was  10  km.  These  variances  and  co-variances  were  then  used  to  calculate 
the  TKE  and  fluxes.  Figure  5.17  shows  that  the  NOMOD  SGS  TKE  and  fluxes  near  the 
surface  are  slightly  larger  in  magnitude  than  those  from  RESDEP.  This  is  expected,  since 
the  smaller  mixing  length  in  RESDEP  directly  results  in  smaller  eddy  exchange 
coefficients  used  in  the  TKE  and  flux  parameterizations  (see  equations  2.8  and  2.11- 
2.14).  Except  for  the  latent  heat  flux,  which  is  simulated  too  large  in  magnitude  in  both 
simulations,  the  calculated  TKE  and  fluxes  from  observations  better  compare  to  those 
from  RESDEP.  This  close  comparison  is  most  clearly  seen  by  inspection  of  the  mean 
fluxes  and  TKE  for  the  observations  and  RESDEP  along  the  near-surface  leg  shown  in 
Figure  5.19  a-d.  However,  the  biggest  differences  are,  again,  in  the  spatial  variability. 
The  RESDEP  fluxes  and  TKE  have  variability  more  closely  resembling  that  of  the 
observations,  and  the  observed  structure  in  EeMone  (1976). 
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Figure  5.18.  From  Chlond  (1992):  vertical  (a)  and  horizontal  (b)  cross-sections  of 
simulated  TKE  from  a  roll  convection  case  using  a  LES  model. 
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c. 


SUMMARY 


Thorough  comparisons  of  the  NOMOD  and  RESDEP  simulations  have  clearly 
shown  that  modification  to  the  mixing  length  in  the  SGS  turbulenee  parameterization 
within  COAMPS™  yields  a  more  realistic  simulation  of  the  BE  roll  structure  observed  in 
the  JES0216  ease  study.  The  mixing  length  adjustment  resulted  in  a  reduetion  of  the 
mixing  of  BE  properties,  and  perpetuated  the  simulation  of  the  roll  vortices  and 
assoeiated  BE  structure,  reflected  in  both  the  mean  and  SGS  model  fields.  The  RESDEP 
near-surface  spectra  for  the  wind,  temperature,  and  moisture  fields  showed  good 
agreement  with  observations  in  resolvable  scale  energy,  especially  near  the  scales  of  the 
rolls,  at  4  to  6  km.  Validation  with  aireraft  near-surfaee  horizontal  leg  and  BE  sounding 
observations  of  the  mean  quantities  showed  good  agreement  with  both  simulations,  but 
RESDEP,  unlike  NOMOD,  simulated  the  BE  wind  shear  associated  with  the  inflection 
point  instability  that  contributes  to  the  roll  eireulation.  RESDEP  SGS  sensible  and  latent 
heat  fluxes  and  TKE  magnitudes  and  structure  in  the  BE  more  realistieally  illustrate  the 
transfer  of  near-surface  heat  and  moisture  in  the  along-roll  direetion  near  the  surface  and 
upwards  in  the  regions  of  the  updraft  eddies,  as  described  by  EeMone  (1976)  and 
simulated  by  Chlond  (1992)  and  Raaseh  (1990)  for  the  roll  convection  case.  The 
magnitudes  of  SGS  turbulent  fluxes  and  TKE  for  the  RESDEP  simulation  were  slightly 
lower  than  the  NOMOD  simulation  and  closer  to  the  observed  fluxes  and  TKE  calculated 
for  scales  at  1  km  and  below,  indicating  a  more  correet  simulation  of  the  conveetive 
forcing  and  turbulence  near  the  surface. 
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Figure  5.19.  SGS  Fluxes  and  TKE,  and  their  means,  derived  from  along-leg  simulated 
data  and  aircraft  observations  from  the  near  the  surface  (40  m  above  sea  level)  near  the 
center  of  JES0216  domain  4.  Shown  in  the  figure  is  (a)  sensible  heat  flux,  (b)  latent  heat 
flux,  (c)  TKE,  and  (d)  momentum  flux  (units  shown)  for  NOMOD  (NM)  (solid), 
RESDEP  (RD)  (dashed)  and  aircraft  observations  (obs)  (dotted). 
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VI.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


A.  SUMMARY  AND  CONCLUSIONS 

This  study  focuses  on  evaluating  the  use  of  ensemble  turbulence 
parameterization  when  the  resolution  of  the  model  reaehes  that  of  the  turbulent  energy 
scale.  The  Navy  COAMPS™  mesoscale  model,  at  500  m  resolution,  is  used  to  simulate 
an  over-ocean  wintertime  ease  (referred  to  as  JES  0216  in  this  study)  involving  a  VCSAO 
over  the  JES  in  whieh  eonveetive  roll  vortices  were  observed.  The  simulation  was  set  up 
over  the  western  JES  region  using  a  four-nest  domain,  where  the  inner  nest,  at  500  m 
resolution  and  100x100  km^  in  size,  eovered  regions  of  BE  aireraft  observations  from  the 
JES  experiment.  Our  modeling  study  focused  on  the  simulated  output  from  the  500  m 
resolution  nest  at  12  hours  into  the  simulation.  Analysis  of  aireraft  observations  from  this 
ease  revealed  that  the  scale  of  the  convective  rolls  was  around  4  to  7  km.  Testing  of  the 
true  spatial  resolvability  of  the  COAMPS™  model  revealed  that  the  resolvable  limit  of 
the  model  was  near  the  seale  of  6 Ax,  or  3  km  for  the  500  m  resolution  nest,  and  thus  we 
expect  the  COAMPS™  high  resolution  nest  to  resolve  these  small  mesoscale  features. 

The  mixing  length,  whieh  is  a  measure  of  turbulence  mixing,  is  used  within  the 
parameterization  for  the  calculation  of  the  SGS  eddy  viseosity  coefficients.  In  the 
COAMPS™  ensemble  turbulence  parameterization,  the  mixing  length  profile  is  based  on 
the  MY82  formulation,  with  the  assumption  that  all  turbulent  motions  are  below  the 
theoretical  resolvable  limit  of  the  model  (2Ax).  The  resolution  of  domain  4  in  the 
JES0216  simulation,  at  500  m,  falls  within  the  scale  of  turbulence,  sinee  the  model  eould 
theoretically  resolve  a  1  km  turbulent  eddy.  Thus,  the  mixing  length  profile,  using  the 
ensemble  assumptions,  is  not  valid  for  the  500  m  resolution  model,  since  part  of  the 
turbulenee  is  teehnically  resolvable.  Using  a  mixing  length  profile  based  on  the  size  of 
the  largest  turbulent  eddy  would  be  inappropriate  for  the  500  m  resolution  nest,  and  result 
in  overestimated  turbulenee  mixing.  It  is  likely  that  this  overestimation  of  the  mixing 
would  result  in  a  reduction  of  the  resolvable  scale  spatial  resolution  of  the  simulated 
features. 
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The  results  of  the  control  run  simulation  (NOMOD)  revealed  that  this  was  indeed 
the  case.  NOMOD,  using  the  MY82  mixing  length  formulation,  adequately  simulated  the 
mean  structure  of  an  over-the-ocean  boundary  layer  during  a  cold  air  outbreak,  but  failed 
to  produce  a  clear  signature  of  the  observed  convective  rolls.  Too  much  mixing  by  the 
model  could  prevent  the  simulation  of  the  distinct  alternating  updraft/downdraft  structure 
and  resulting  longitudinal  alignment  of  the  roll  eddies  by  the  mean  BL  wind.  NOMOD 
simulated  very  little  spatial  variation  in  the  mean  fields,  with  no  clear  signature  of  roll¬ 
like  convective  eddy  structure  present  in  both  horizontal  and  vertical  cross-sections 
within  the  BL,  the  result  of  a  very  weak  and  unstructured  simulated  vertical  velocity 
field.  The  small  variability  present  at  mid-BL  levels  shows  no  roll-like  structure,  as 
viewed  in  the  horizontal  plane.  Visual  inspection  of  BL  cross-sections  of  the  mean  fields, 
supported  by  spectral  analysis  of  the  simulated  fields  in  comparison  with  surface  layer 
aircraft  observations,  revealed  that  the  NOMOD  simulation  was  not  adequately 
simulating  energies  less  than  about  10  km,  and  not  the  expected  3  km  resolvability  or 
even  the  4  to  7  km  scale  range  of  the  observed  rolls.  The  simulated  SGS  turbulent  flux 
and  energy  structure  was  fairly  uniform  across  the  domain  throughout  the  depth  of  the 
BL,  signifying  the  over-mixing  by  the  turbulence  parameterization  scheme. 

Modifications  are  made  to  the  ensemble  turbulence  parameterization  in 
COAMPS™  through  a  modification  to  the  turbulence  mixing  length  scale.  It  is  hoped 
that  the  new  mixing  length  would  yield  the  appropriate  spatial  resolution  of  the  mean 
fields  in  the  resolvable  scale,  and  an  appropriate  representation  of  the  turbulent  structure 
in  the  SGS.  The  modification  is  based  on  observational  analysis  of  the  integral  length 
scale  (analogous  to  the  mixing  length)  of  over-the-ocean  convective  boundary  layers, 
which  illustrated  that  the  integral  scale  profile  in  the  BL  expectedly  decreases  in 
magnitude  as  the  maximum  scale  of  turbulent  eddies  being  considered  decreases.  From 
the  integral  scale  analysis,  we  derived  a  resolution-dependent  relationship  that  calculates 
the  appropriate  integral  scale  profile  in  the  BL  for  the  maximum  scale  of  turbulent  eddies 
considered.  The  equation  for  the  integral  length  scale,  as  applied  to  the  results  from 
observations,  is  based  on  taking  a  weighted  average  between  the  integral  length  scale 
calculated  for  the  smallest  scales  of  turbulence  considered  (200  m)  and  the  integral  length 
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scale  appropriate  for  the  largest  turbulent  eddies  (6  km),  where  the  weighting  is 
dependent  on  the  scale  being  considered. 

The  resolution-dependent  integral  length  seale  formulation  developed  from  the 
results  of  our  observational  study  is  adapted  for  use  in  a  numerical  model,  and  applied 
within  COAMPS™  for  the  JES0216  case.  This  resolution-dependent  mixing  length 
formulation  in  COAMPS^'^  uses  a  weighted  average  between  a  mixing  length  appropriate 
for  turbulent  eddies  within  the  inertial  sub-range  and  the  original  COAMPS™  MY82 
mixing  length  formulation,  where  again  the  weighting  is  based  on  the  domain  resolution. 
The  resolution-dependent  mixing  length  equation  is  then  applied  only  to  the  500  m 
resolution  inner  nest  in  the  RESDEP  simulation,  where  the  appropriate  mixing  length 
profile  is  calculated  by  considering  mixing  by  only  the  turbulent  eddies  at  and  below  1 
km. 

The  RESDEP  simulation  resulted  in  much  better  simulation  of  the  roll  strueture  in 
all  respects.  RESDEP  simulated  turbulent  roll  eddies  in  both  the  momentum  and  scalar 
mean  fields  that  are  consistent  with  the  aireraft  observations  and  with  past  observational 
and  numerical  studies.  The  model  simulated  a  mueh  stronger  and  well-defined  eddy 
structure  in  the  vertical  velocity  field  compared  to  NOMOD,  with  warm,  moist  updrafts 
and  cooler,  drier,  downdrafts,  aligned  with  the  mean  BE  wind  direction  into  longitudinal 
rolls,  which  scaled  eomparatively  to  the  size  of  the  observed  rolls.  RESDEP  mean  wind 
profiles  illustrated  significant  wind  shear  within  the  BE,  with  an  inflection-point  in  the 
wind  components  near  the  middle  of  the  BE.  This  simulated  shear  signature,  that  was 
present  in  the  observation  profiles,  but  not  present  in  the  NOMOD  wind  profiles, 
indicates  inflection-point  instability  as  a  contributor  to  production  of  the  rolls.  Energy 
speetra  of  the  RESDEP  surfaee  layer  mean  fields  showed  improved  spatial  resolution  in 
the  convective  roll  seale  range,  comparing  to  the  observations  within  one  order  of 
magnitude.  SGS  flux  and  TKE  fields  were  also  more  realistic  in  comparison  to  NOMOD 
for  the  roll  structure,  where  the  flux  and  TKE  modulation  depicted  the  redistribution  of 
turbulence-producing  elements  by  the  mean  roll  eirculation.  This  structure  was  not 
validated  with  observations,  due  to  the  laek  of  spatial  eoverage,  but  was  supported  by  the 
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observational  study  of  LeMone  (1976)  and  numerical  studies  by  Chlond  (1992)  and 
Raasch  (1990). 

The  convective  roll  case  of  JES0216  appropriately  illustrated  the  issue  that 
appropriate  mixing  is  required  to  achieve  the  correct  simulation  of  the  convective  roll 
structure  in  both  the  simulated  mean  fields  and  SGS  fields.  The  implications  for  the  use 
of  a  modified  mixing  length  formulation  within  COAMPS^'^  is  encouraging  for  the 
simulation  of  other  mesoscale  features  whose  development  and  maintenance  via  the 
interaction  with  eddies  in  the  turbulent  scale  is  important.  However,  further  study  to 
completely  understand  this  issue  and  to  improve  the  parameterization  is  obviously 
desired. 

B.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  conclusion  that  the  change  in  mixing  length  profile  within  COAMPS™ 
resulted  in  improvements  to  the  simulated  BL  structure  is  thus  far  tested  with  the  single 
case  of  JES  0216  for  a  resolution  of  500  m.  Eurther  study  is  required  to  determine  if  the 
formulation  for  mixing  length  presented  here  is  adequate  for  general  use  in  COAMPS^'^. 
Pirst,  the  modified  parameterization  should  also  be  tested  for  this  case  at  even  higher 
resolutions,  to  determine  if  the  formulation  is  valid  in  the  limit  as  Ax  approaches  the  100 
m  range.  This  test  was  not  possible  here,  as  the  lack  of  computer  resources  prohibited  a 
simulation  of  an  adequately  sized  domain,  with  resolution  below  500  m,  for  proper 
analysis  of  the  issue.  Model  analysis  becomes  more  difficult  when  the  resolution  of  the 
model  falls  below  200  or  300  m,  where  large  turbulent  eddies  are  above  the  6 Ax 
resolvable  limit  of  the  model.  In  this  case,  both  SGS  and  resolved  scale  turbulence  and 
fluxes  must  be  considered  in  the  model  analysis,  as  in  the  case  of  EES  studies.  In 
addition  to  testing  of  different  resolutions,  sensitivity  testing  of  the  formulation  using 
different  values  of  minimum  mixing  length  limits  and  different  values  of  the  constant 
parameter  R  is  also  suggested.  Different  ensemble  turbulence  parameterization  schemes 
could  also  be  tested,  where  a  formulations  for  the  mixing  length  and/or  the  eddy  mixing 
coefficients  other  than  the  MY82  is  used.  In  COAMPS™,  the  MY82  scheme  is  most 
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commonly  used  for  mixing  length,  but  other  built  in  options  are  available,  sueh  as  the 
Therry  and  Laearrere  (1983)  formulation. 

In  addition,  we  reeommend  testing  of  the  validity  of  this  mixing  length 
modifieation  for  mesoseale  features  of  varying  type  and  seale,  for  oases  over  ooean  and 
over  land.  In  JES02I6,  the  feature  being  simulated  was  very  olose  to  the  6Ax  resolvable 
limit  of  the  model  at  500  m  resolution.  Results  may  differ  as  the  seale  of  the  feature 
oompared  to  the  resolvable  limit  of  the  model  inoreases.  Cases  where  olouds  form  may 
also  yield  different  results,  sinoe  the  parameterization  soheme  in  COAMPS^'^  handles 
fluxes  differently  when  olouds  are  present  in  the  BL.  For  more  oomplex  mesoseale 
phenomena  that  are  more  intense,  sueh  as  squall  lines,  oonveotive  oell  thunderstorms,  or 
sea-breeze  oonvergenoe  for  instanoe,  the  real  data  simulation  may  not  allow 
improvements  to  be  easily  seen.  Testing  of  the  formulation  using  idealized,  vioe  real  data 
simulations,  may  present  olearer  differenoes  between  oontrol  runs  and  modified  runs  for 
these  more  oomplex  oases. 

The  JES  0216  ease  model-observation  validation  was  limited  by  the  laok  of 
spatial  ooverage  of  the  airoraft  observations,  where  only  BE  soundings  and  surfaoe  layer 
flux  legs  were  available.  Cases  where  more  oomplete  airoraft  data  sets  are  available, 
oontaining  long  level  flight  legs  at  different  levels  within  the  BE  and  soundings  to  greater 
vertioal  levels  would  allow  for  a  more  oomplete  inter-oomparison  of  the  simulated  data 
and  observations.  High  resolution  satellite  imagery,  allowing  for  seale  oomparisons  of 
simulated  mesoseale  oloud  features,  would  also  prove  useful. 
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